Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7190TA45.5% 36.0 / 33.1 44intergenic (‑348/‑45)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (18/6);  new base A (0/20);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.57e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGC  >  NC_000913/4287570‑4287861
                                                                                                                                                       |                                                                                                                                               
aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAtc                                                                                                                                                <  2:123535/151‑2 (MQ=255)
aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAt                                                                                                                                                 >  2:713605/1‑150 (MQ=255)
 aaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATaa                                                                                                                                               <  1:291448/151‑1 (MQ=255)
  aaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATaa                                                                                                                                               >  1:879867/1‑150 (MQ=255)
    gagGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATaaaa                                                                                                                                             >  1:888291/1‑150 (MQ=255)
                    aacttaaaggtgtTATTTTACATGCACTTACAATTGATAAAAGACAACATTAACGGTGTGGTTATTTTTTACACAAAGGGGCGAAAAAGGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAAt                                                                                                                           <  1:161620/142‑1 (MQ=255)
                    aaCTATAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAAt                                                                                                                           >  1:128953/1‑151 (MQ=255)
                      ctctAAAGTGTTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACGGTGTGGTTATTTGTTACACATGGGGGCGAGCAATGTCATGACAGTGTGGGGGCGT‑TTACTG‑TCGTAAAAAACAAAAAAAA‑CCGCCATTGCAACAATg                                                                                                                          <  2:394747/150‑1 (MQ=255)
                        ctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGgcg                                                                                                                       >  1:439139/1‑151 (MQ=255)
                             gTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAAtt                                                                                                                  >  1:803802/1‑151 (MQ=255)
                             gTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAAtt                                                                                                                  >  1:543911/1‑151 (MQ=255)
                               gttATTTTACATGCATTTACAATTGATTAAAGAAAACATTAACAGTGTGGTTTTTTGTTCCCCATAGGGGCAAGCAATGTAATGACAGGGTGGGTGCGG‑TTACTG‑TCGTAAAAAAAAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTc                                                                                                                 <  1:135019/148‑1 (MQ=255)
                                  aTTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGAt                                                                                                             >  1:478580/1‑151 (MQ=255)
                                    tttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATg                                                                                                            >  1:646436/1‑150 (MQ=255)
                                    tttACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTTTTTTTTACACATGGGGGCGAGCAATGTCATGACAGTGGAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATg                                                                                                            <  1:543863/150‑1 (MQ=255)
                                          tgaccTTAAAATTGATTAAAGAAAACATGCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTTTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGccc                                                                                                      <  1:101047/146‑1 (MQ=255)
                                           gCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGT‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGcccc                                                                                                     <  2:882840/150‑1 (MQ=255)
                                             aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                  >  2:120177/1‑151 (MQ=255)
                                                 aCAATTGATAAAAGACAACATTCACAGTGTGGTTATTTTTTACACATAGGGGAGAGCAATGTCATGACAGTGTAGGTGCGT‑TTACTG‑TCGTGAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc                                                                                               <  2:80705/150‑1 (MQ=255)
                                                    aTTAATAAAAGACAACATTCACAGGGTGGTTATTTGTTACACATAGGGGCGAGAATTGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAgg                                                                                           <  2:202167/151‑1 (MQ=255)
                                                        aTTAAAGACAACATTCACAGTGTGGTTATTTTTTACAAAAAGGGGCGAGAAATGTCATGACAGTGAAGGTGGGG‑TTACTG‑TCGTAAAAAAAAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaa                                                                                       <  1:506690/151‑1 (MQ=255)
                                                            aaaacaacaTTCAAAGTGTGGTTATTTGTTACACATGGGGGCGACCATTGTCATGACAGTGTAGGTGCGT‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAAt                                                                                   <  2:491560/148‑1 (MQ=255)
                                                               acaacaTTCACAGTGTGGTTATTTTTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATaaa                                                                                <  2:846760/151‑1 (MQ=255)
                                                                 aacaTTCACAGTGTGTTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac                                                                               <  2:886179/150‑1 (MQ=255)
                                                                               gttttttttttACAAATGGGGGAGAGCATTTTAATGACAGT‑TAGGGGGGG‑TTTATGTTCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATCCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTa                                                                     <  1:225834/141‑1 (MQ=255)
                                                                                 ttttttGTAACAAAAGGGGGGGAGAAATGTCATGAAAGTGGAGGTGGGG‑TTACTG‑TCGTAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCtt                                                                 <  2:639802/145‑1 (MQ=255)
                                                                                 ttttttGTAACAAAAGGGGGCGAGAAATGTCAGGACGGTGTAGGGGGGT‑TTACTG‑TCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa                                                               <  1:146825/147‑1 (MQ=255)
                                                                                     tttttACAAATGGGGGCGAACAATGTCATGACGGTGTGGGTGGGT‑TTACTG‑TCGTAAAAAACAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                           <  2:414692/147‑1 (MQ=255)
                                                                                        ttACACATGGGGGCGAAAAAGGTCATGACGGTGTAGGTGCGT‑TTACTG‑TCGGAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt                                                         <  1:110386/149‑1 (MQ=255)
                                                                                           cacaggggggCGAGAAATGTAATGAGGGGGGGGGGGGGT‑TTACGG‑TCGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                                    <  2:872830/145‑1 (MQ=255)
                                                                                              aTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt                                                 >  1:256326/1‑151 (MQ=255)
                                                                                               tAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtt                                                >  2:4096/1‑151 (MQ=255)
                                                                                                 gggggaagaattGTCATG‑CAGTGGTGGGGCGGTTTACTG‑T‑GTAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt                                              <  1:1467/140‑1 (MQ=255)
                                                                                                 ggggCGAAAAATTTAAGGACGGTGTGGGGGGGT‑TTCCTG‑TGGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttt                                               <  1:128085/150‑1 (MQ=255)
                                                                                                 ggggCGAAAAATGTCATGACAGTGTGGGTGCGT‑TTACTG‑TCGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                                    <  1:346519/145‑1 (MQ=255)
                                                                                                 ggggCGAAAAATGTAATGACGGTGTGGGGGCGG‑TAACTG‑TGGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt                                              <  1:197988/151‑1 (MQ=255)
                                                                                                          aaTGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                                     >  1:19295/1‑151 (MQ=255)
                                                                                                                aTGACAGTGT‑GGTGGGGTTTACTG‑CGGAAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                               <  1:869496/151‑1 (MQ=255)
                                                                                                                        gtAGGTGCGT‑TTCCGG‑TGGAAAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACg                        <  2:91265/150‑1 (MQ=255)
                                                                                                                                  gtttcgtgtaataaaaaaCAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAAc                  <  2:519317/135‑1 (MQ=255)
                                                                                                                                     aCTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc           >  1:294098/1‑151 (MQ=255)
                                                                                                                                          tggaaaaaaaaaaaaaaaa‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct          <  1:513511/143‑1 (MQ=255)
                                                                                                                                           cGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTc      >  1:771293/1‑151 (MQ=255)
                                                                                                                                            gTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTcc     >  2:819441/1‑151 (MQ=255)
                                                                                                                                              gAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCg   >  1:765515/1‑151 (MQ=255)
                                                                                                                                              gAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCg   >  1:693807/1‑151 (MQ=255)
                                                                                                                                               aaaaaCAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGc  >  2:545210/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
AAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGC  >  NC_000913/4287570‑4287861

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