Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,892,1950AT41.2% 33.4 / 14.2 34intergenic (‑297/‑19)cysJ/queDsulfite reductase, flavoprotein subunit/6‑carboxy‑5,6,7,8‑tetrahydropterin synthase
Reads supporting (aligned to +/- strand):  ref base A (7/13);  new base T (14/0);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 1.60e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATG  >  NC_000913/2891953‑2892439
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
aaaCAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTg                                                                                                                                                                                                                                                <  1:85408/251‑1 (MQ=255)
     gaagagaagGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTagag                                                                                                                                                                                                                                           <  2:233591/251‑1 (MQ=255)
                         ggggCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTGTTGAGAAAtttttattgtcaccacg                                                                                                                                                                                                                       >  1:168683/1‑236 (MQ=255)
                            gCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTGTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTa                                                                                                                                                                                                                    <  2:74569/251‑1 (MQ=255)
                                 cGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTTGTCCTGTAGAGAAATTATGATGTCCACCCCGttttttt                                                                                                                                                                                                                >  1:29324/1‑249 (MQ=255)
                                 cGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTGTAGAGAAATTATGATTTCCACCCCGGTATTTa                                                                                                                                                                                                                >  1:275779/1‑250 (MQ=255)
                                          ttAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTTGTCCTGTAGAGAAATTATGATTTTCACCCCGTTATTTAAAGATTTCAc                                                                                                                                                                                                      >  1:4697/1‑251 (MQ=255)
                                           tAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTGTAGAGAAATTATGATGTCCACCCCGGTATTTAAAGATTTTAcc                                                                                                                                                                                                     >  1:359894/1‑251 (MQ=255)
                                                         cGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCGACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCa                                                                                                                                                                                      <  2:256135/251‑1 (MQ=255)
                                                                             gACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATTGTCCTGGAGAGCAATTATGGTGTCCACCCCGTTTTTTTAAGATTTTACCTTT‑CAAGACGCTCCACGCTTACCCCACGTccacg                                                                                                                                                                  >  1:196901/1‑248 (MQ=255)
                                                                                tAACCTGTTTTAGCACTAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAACCAACTGAATTTAGCCAATCATTAAGATACATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTATTTTTGTCCTGTTGAGGAAATATTATGTACTCCCCGTTACTTTAAGATTTTATCTTCAGAA‑CCGCTCACCGCTTCCcacccccccccgg                                                                                                                                                                 >  2:340925/1‑239 (MQ=255)
                                                                                             cAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTTGTCCTGTAGAGAAATTATGATGTCCACCCCGTTATTTAAAGATTTCACCCTC‑GAAGCCGCTCACCGCTTACACCACGTCCCGGAAGGGCATAAATGTgg                                                                                                                                                  >  1:125014/1‑251 (MQ=255)
                                                                                                       tAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATTGTCCTGGAGAGAAAATTTGATGTTCACCACGTTATTTTAAGATTTCACCCTC‑GAAACCGCCCACCCCCTACCCCACCTCCCCGGAGGGCATAAAAGTGGTCGCCTGccc                                                                                                                                        >  2:137961/1‑249 (MQ=255)
                                                                                                               aaaTAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTTGGCCTGTAGAGAAATTATGATGTTCCCCACCTTTTTTAAAGATTTTACCCTC‑GAAGCCGCCCACCCCTTTCCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGccacggccttc                                                                                                                                >  2:316951/1‑241 (MQ=255)
                                                                                                                        aaccaaccTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCATCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATgg                                                                                                                       <  2:252204/250‑1 (MQ=255)
                                                                                                                                       aaaCAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCCACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGGTATTTAAAGATTTCACCTTC‑GAAGACGCTCACCGCTTACCCCACGTCCCGGGAAGGCATAAATGTGGTCGCCTGCAAGGGCA‑TTCCTTTTTGGTGGGACTGGACATTAc                                                                                                       >  1:319738/1‑251 (MQ=255)
                                                                                                                                         aCAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTTTAGAGAAATTATTATGTCCACCAAGGTTTTTAAAGAATTTACCTTC‑GCAGCCGGCCACCGCTTACCACACCTCCCCGAAGGGGATAAAAGTGGTCGCCTGCCCCGGCA‑TTCCTTTaatgtggggatggaaaataacg                                                                                                     >  2:298226/1‑230 (MQ=255)
                                                                                                                                                    tAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCCCGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCCCTTACCACACGTCCCGGAAGGGGATAAATGTGGTCGCCTTCCCGGGCA‑TTCCTTTATGGTGCGACTGGAAATTAACGGGGAAGTCGAt                                                                                          >  1:155446/1‑251 (MQ=255)
                                                                                                                                                     aGCCAATCATTAAGGTAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTTGGTCAGTTGAGAAATTTTGATGTCCACCCCGGTTTTTTAAGATTTTACCTTC‑GAAGGCGCCCACCCCCTACCACAACTACACGAAAGGCATAAATATGGTCGGCTGCACGGGCAATTCTTTTTTGGT‑CGGCTGGGAAATACCgggggaggggaa                                                                                          >  2:172403/1‑242 (MQ=255)
                                                                                                                                                       ccAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCg                                                                                       <  1:2628/251‑1 (MQ=255)
                                                                                                                                                                       aTCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTTTAGAGAAATTATGATTTCCACCACGTTATTTTAAGATTTCACCTTC‑CAAAACGCCCACCCCTTACCAAACCTCCCGGAAGGGCATAAATGTGGTCGCCCGCACGGGCC‑TTCCTTTTTGGTGCGACTGGAAATTACCCGGGGAGGCGATccccctactgggtgggttt                                                                       >  2:95374/1‑234 (MQ=255)
                                                                                                                                                                          aGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTAtag                                                                    <  2:111519/251‑3 (MQ=255)
                                                                                                                                                                                        acaaAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCCCCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTTGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTa                                                      <  2:309527/251‑1 (MQ=255)
                                                                                                                                                                                          aaaGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTGTAGAGAAATTATGATGTCCCCCACGTTATTTTAAAATTTTACCTTT‑GAAGCCGCTCACCGCTTTCCACACGTCCCCGAAAGGGATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGGGCGACTGGGAATTACCGGGGAAGTCGATCCGCATACGGGCTGGGTTATCGATTTTGCTGAACTaa                                                     >  1:355025/1‑250 (MQ=255)
                                                                                                                                                                                          aaaGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAACCGCCCCACGCTTACCACACGGCCCGGGAGGGGATAAATGTGGGCGCCCGCCCGGGCA‑TTCCTTTTTTGGGCGGCTGGAAATTACCGGGGAAGGCGATCCGCATACGGGCTGGGTTTTCGATTTTGCTGGACTaa                                                     >  1:139935/1‑250 (MQ=255)
                                                                                                                                                                                               cGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTTGGCCTGTAGAGAAATTATGGTGTCCACCACGTTATTTAAAAATTTTACCTTC‑GAAGCCGCTCACCGCCTACCACACGTTCCGGAAGGGGATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATTGGGCGACTGGAAATTTACGGGGAAGTCGATCCGCAAACGGGCTGGATTTTCGGTTTTACTTAACTAAAAgcgg                                               >  1:284832/1‑251 (MQ=255)
                                                                                                                                                                                                   ttAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGATTTATCGTCATATAGAGAAATTATGATGTCCACCACGGTATTTTAAGATTTCACCCTA‑CAAGCCCCCCACCGCCTACCACACGGCCCGGAAGGGCGTAAATGGGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTGGGAATTACCGGGGGAGTCGATCCGCATACGGGGTGGAGTATCGATTTCGCTGAATTAAAAgcggcggt                                           >  1:249409/1‑249 (MQ=255)
                                                                                                                                                                                                    tAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGtt                                           <  2:34109/250‑1 (MQ=255)
                                                                                                                                                                                                      gATACTGCGTAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATTGTCCTGTAGAGAAATTATGATTTCCACCACGTTATTTTAAGATTTCACCTTC‑CAAGCCGCCCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGGCGCCTGCACGGGCA‑TTCCttt                                                                                                                           >  1:270678/1‑168 (MQ=255)
                                                                                                                                                                                                      gATACTGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCttt                                                                                                                           <  2:270678/168‑1 (MQ=255)
                                                                                                                                                                                                      gATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTaa                                        <  1:74815/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaCAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCCCCCCGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCCCACGTCCCGGGAGGGCATAAATGTGGTCGCCTGCACGGGCA‑CTCCTTTATGGGGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAAATAAAAACGGCGTTTAAACCAACCtccg                             >  1:46253/1‑246 (MQ=255)
                                                                                                                                                                                                                                          gcggTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCGCGATCACCATTATCTCaa    <  2:33014/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ggTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATg  <  2:55305/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
AAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTC‑GAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCA‑TTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATG  >  NC_000913/2891953‑2892439

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