Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,022,9430CA36.6% 53.6 / 17.0 41intergenic (+52/‑447)doeA/proP_6Ectoine hydrolase/Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base C (19/7);  new base A (0/15);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 2.75e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACT  >  ACB122/5022811‑5023076
                                                                                                                                    |                                                                                                                                     
aTTACCGAGAGCATCATTATAACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtt                                                                                                                                   >  3:89708/1‑136 (MQ=255)
   aCCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGa                                                                                                                                <  3:52717/137‑1 (MQ=255)
            aTCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCtt                                                                                                                       >  1:220225/1‑137 (MQ=255)
                            cTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTcc                                                                                                       >  4:264945/1‑137 (MQ=255)
                            cTGGGGGCGAGACGTTCTGCTCTACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTcc                                                                                                       <  5:149042/137‑1 (MQ=255)
                                  tCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCAc                                                                                                 >  3:106776/1‑137 (MQ=255)
                                       aCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCtt                                                                                            >  6:75012/1‑137 (MQ=255)
                                              gCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCCGCCACTCCTTTGTTTTa                                                                                     >  6:230247/1‑137 (MQ=255)
                                                      cctcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACaaaaaaa                                                                             >  1:78271/1‑137 (MQ=255)
                                                       ctcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAg                                                                            >  6:36654/1‑137 (MQ=255)
                                                       ctcGAAAGCTTTTCGTCACGTATTTACTAAAAGTTTGAAAGTCTCTGGAGTTATTTTGACACGGCTATCTTCAGCTAAAATGTTATGGTGCCTAAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACaaaaaagc                                                                            <  2:218684/137‑3 (MQ=255)
                                                        tcGCACGCTTTTGGTCATGGACTGTCTAAAATTTTGCATGTGGCTGTTTTTTTTGTGCAACGGCTATATTCAGCTAAATTGTGTTGGTGCCTTAATGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGc                                                                           <  8:109928/137‑1 (MQ=255)
                                                        tcGCAAGCTTTTCGTCAAGGACTTACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGAAAATGTTATGGTGCCTTAAAGCACAATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGc                                                                           <  7:272855/137‑1 (MQ=255)
                                                        tcGAAAGCTTTTCGTCAAGGACTGACTAAAAGTTTGAATGTCGCTGGAGTTATAGGGCCCCGGCTATCTTCAGCGAAAATGTGATGGTGCCTTAAACCACAATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGc                                                                           <  6:7016/137‑1 (MQ=255)
                                                         cGCAAGCTTTTCGTCAAGGACTGACTAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCTACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCa                                                                          <  5:114715/137‑1 (MQ=255)
                                                          gCAAGCTTTTCGTCAAGTACTGACTAAATGTTTGAAAGTCGCTGGATTTATAGTGCAACGGCTATATTCAGCGAAAATGTTATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCAc                                                                         <  8:17677/137‑1 (MQ=255)
                                                           cAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGACTTAAAGCACAATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGAAcc                                                                        >  5:59584/1‑136 (MQ=255)
                                                            aaGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGc                                                                       >  8:44989/1‑137 (MQ=255)
                                                                     gTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGaa                                                              >  6:21190/1‑137 (MQ=255)
                                                                          ggACTTACTAAATGTTTGCAAGTCGCTGTAGTTATAGTGCCACGGCTATCTTCAGCGAAAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGt                                                         <  4:106776/137‑1 (MQ=255)
                                                                                   aaatttttGAATGTCGTTGGAGTTATAGTGCCCCGGATATCTTCAGCGAAAATGTCAGGGGGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTg                                                <  7:20546/132‑1 (MQ=255)
                                                                                   aaaaGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGttt                                                >  1:166443/1‑136 (MQ=255)
                                                                                   aaaaGTTTGCAAGTAGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGAACCATCACTCCAGCCAATCCTTTGTTTTACAAAAAAAGCACTCAGCCTTGAAGAGGTTGAAAGTTg                                                >  6:117645/1‑137 (MQ=255)
                                                                                         ttGCAAGTCGCTGGAGTTATAGTCCCCCGGCTATATTCAGCGAAATTGTGATGGTGCCTTAAACCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGtt                                          <  1:242217/137‑1 (MQ=255)
                                                                                              aGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCa                                     <  5:95915/137‑1 (MQ=255)
                                                                                               gTCGCTGGAGTTATTGTGCCACGGCTATATTCAGCTAAAATTTCATGGTGCCTAAAACAACAATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAg                                    <  3:264945/137‑1 (MQ=255)
                                                                                               gTCGATGGAGTATAAGTCCCACGGCTAAATTCAGCGAAAATGTTATGGTGCCTAAAACCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAg                                    <  1:181415/137‑1 (MQ=255)
                                                                                                    tGGAGATATAGTGCCACGGCTATCTTCAGCTAAAATGTGATGGTCCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGc                               <  1:58283/137‑1 (MQ=255)
                                                                                                                  cACGGCTTTCTTCACCTAAATTGTTATGGTGCATAAAAGAACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGc                 <  3:120069/137‑1 (MQ=255)
                                                                                                                  cACGGCTATCTTCAGCTAAAATGTGATGGTGCCTTAAAGCACAATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGc                 <  5:71896/137‑1 (MQ=255)
                                                                                                                   ccGGCTATCTTCAGCGAAAATTTGATGGTGCCTTAAACCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGcc                <  3:79089/136‑1 (MQ=255)
                                                                                                                   aCGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGcc                >  3:14273/1‑137 (MQ=255)
                                                                                                                    cGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCt               >  5:286252/1‑137 (MQ=255)
                                                                                                                    cGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCt               >  5:190662/1‑137 (MQ=255)
                                                                                                                    cGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCt               >  4:128621/1‑137 (MQ=255)
                                                                                                                      gCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGc                 >  4:172579/1‑133 (MQ=255)
                                                                                                                      gCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGc                 <  3:172579/133‑1 (MQ=255)
                                                                                                                            ttCAGCTAAAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTg       <  3:74581/137‑1 (MQ=255)
                                                                                                                             taaGCGACAATGTGATGGTGCCTTAAAGAACCATCACTCCAGCCAATCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTTATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGa      <  6:59584/135‑1 (MQ=255)
                                                                                                                             tCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGa      <  6:286252/137‑1 (MQ=255)
                                                                                                                                 cGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACt  >  8:47164/1‑137 (MQ=255)
                                                                                                                                    |                                                                                                                                     
ATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACT  >  ACB122/5022811‑5023076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑