Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7100GA29.0% 50.0 / 10.1 31intergenic (+43/‑54)ytcB/nrfAprotein YtcB/cytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base G (18/4);  new base A (0/9);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 3.55e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCG  >  NC_000913/4287462‑4287957
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        
gATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGa                                                                                                                                                                                                                                                         >  2:81509/1‑250 (MQ=255)
      gCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGCAAAA‑Ca                                                                                                                                                                                                                                                  >  2:58077/1‑250 (MQ=255)
      gCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑c                                                                                                                                                                                                                                                   <  1:156211/249‑1 (MQ=255)
      gCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑c                                                                                                                                                                                                                                                   >  2:156211/1‑249 (MQ=255)
                   tgtgCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCAt                                                                                                                                                                                                                                    >  1:23778/1‑251 (MQ=255)
                                      gTAATTTTTAAAGGTCAGAAAAAAAATTTGATTAATTCTTGGAGAAACATGCGGTAATGAAGGCTGTTGAAAAAGAGGAAGTTCCTGACTAAC‑CTAAAAGTGGTATTTTACATGCACTTAAAATTGATTAAAGACACCATCCACAGTGTGGTT‑GTTT‑TTTCCGCAT‑AGGGGC‑GAGAAATTTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                                                                                                                                                 <  1:289277/251‑1 (MQ=255)
                                              tGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGa                                                                                                                                                                                                           >  2:139642/1‑249 (MQ=255)
                                                                             ttGAGGAACATACAGTTATGCATGCTTTTGAAAAAGAGGAAGATCCTGACTAACTCTAAA‑GTGGTATTTTATATGCACTTCCAATTGATTAAAGACAAACTTCACAGTGTGGTT‑TTTT‑GTTACACAT‑AGGGGC‑GATCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac                                                                                                                                                                           <  1:267687/250‑1 (MQ=255)
                                                                                         cAGTTATGCATGCTGTTGAACAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GAG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCAATTGCAACAATGGCGCAATTCGGATGAAGCCCCTCTGACAAGGCTAAAAATAAACGAACGCCGTAGCTTc                                                                                                                                                                >  2:229510/1‑249 (MQ=255)
                                                                                              atgcatgcTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt                                                                                                                                                         >  2:206701/1‑251 (MQ=255)
                                                                                                            aaaaGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc                                                                                                                                           >  1:117721/1‑251 (MQ=255)
                                                                                                                            aaaaaactaaaaa‑GTGGTTTTTTAAATCCATTAAAAATTAATAAAAAACAACATTCACAGT‑TGGTTTATTT‑GTTA‑AAATGAGGGGA‑GAGAAATGTAAGG‑ACAGTGGAG‑GTG‑GGGTTACTGTCGTAAAAAA‑AAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc                                                                                                                           <  2:21989/241‑1 (MQ=255)
                                                                                                                             aCTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc                                                                                                                           >  2:57837/1‑250 (MQ=255)
                                                                                                                                     tAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAAc                                                                                                                  >  1:321145/1‑251 (MQ=255)
                                                                                                                                        a‑GTGTTATTTTACATGCATTAAAAATTGATAAAAGACAACTTTCAAATTGTGGTT‑ATTT‑GTTCCACAT‑GGGGGG‑GAGAAATGTCATG‑AAGGGGTGG‑GTG‑GGGTTCCTGT‑GTAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaaa                                                                                                               <  1:93285/251‑1 (MQ=255)
                                                                                                                                        a‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGA‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa                                                                                                                >  1:380094/1‑250 (MQ=255)
                                                                                                                                                                    gataaaaGACAAAATTCACAGTGTGTTT‑ATTT‑GTTAAAAAG‑GGGGGC‑GAGAAATGTCAGA‑ACAGTGGGG‑GTG‑GGGTTATTGTCTTAAAAAA‑AAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACt                                                                                     <  1:81509/246‑1 (MQ=255)
                                                                                                                                                                      ttAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGgc                                                                                                                                                                                                                        >  1:316857/1‑117 (MQ=255)
                                                                                                                                                                          aagaaaacttccacggttgGTT‑ATTTTGTAACACAG‑GGGGGC‑GA‑AAATGGCATGGACAGTGGGG‑GTG‑CGTTTACTGACGAAAAAAA‑AAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGaa                                                                               <  2:321145/234‑1 (MQ=255)
                                                                                                                                                                                accatccacattttgttttttttttACACAA‑AGGGGG‑GACAAATGCCATG‑ACAGTGTGG‑TTG‑CGGTTACTGTCTTAAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCgaa                                                                          <  1:139642/228‑1 (MQ=255)
                                                                                                                                                                                  gtcatagattgGTT‑ATTT‑GTGACACAG‑GGGGGC‑GAACAATGTCAAG‑AAAGTGTTG‑GTT‑AGGTTGCTATAGTAAAAAA‑AAATAATACCCGCCATTGCAACAATGGCGCAATGCGGATGAAGCCCCAATGACAAGGATAAAAATAACCGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCgaattat                                                                      <  2:446789/242‑5 (MQ=255)
                                                                                                                                                                                   tcgggtggtgGTT‑TTTT‑TTTACACAT‑AGGGGC‑GACTAATGTAATG‑AATGTG‑AGTGTG‑CGTTTACTTTGGAAAAAAA‑AAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATTAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCGCCCGCAAAACCTGTAACTGTGGAAGCGAAGAAt                                                                      <  2:172378/243‑1 (MQ=255)
                                                                                                                                                                                    cacaGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATg                                                                     >  1:396822/1‑250 (MQ=255)
                                                                                                                                                                                     acaGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGa                                                                    >  1:380548/1‑250 (MQ=255)
                                                                                                                                                                                      aaGTGTGGTT‑TTTT‑GTTACAAAT‑GGGGGG‑GAGAAATGTCATG‑ACGGTGTGG‑GTG‑CGGTTACTGTCGTAAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGaaa                                                                  <  1:307069/250‑1 (MQ=255)
                                                                                                                                                                                       ggtgtGGTT‑ATTT‑GTTACACA‑‑GGGGGCTGAGAAATGTAAGC‑ACGGTGGGG‑GTG‑CGGTTACTTCCGAGAAAAA‑CAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGaaa                                                                  <  1:108290/249‑1 (MQ=255)
                                                                                                                                                                                             ggtgttttgtgtcaaaaatgggggA‑GAGAAATGTCATG‑A‑TGTGTTG‑GTGGCGTTTACCGTCGAAAAAAA‑CAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCtt                                                              <  2:178549/231‑1 (MQ=255)
                                                                                                                                                                                                        tACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCgcagc                                                   >  2:180623/1‑250 (MQ=255)
                                                                                                                                                                                                                                             g‑gTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCg                >  1:320692/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       tACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTc        >  1:290176/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             cGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCg  >  1:22260/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        
GATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAA‑GTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTT‑ATTT‑GTTACACAT‑AGGGGC‑GAGCAATGTCATG‑ACAGTGTAG‑GTG‑CGGTTACTGTCGTGAAAAA‑CAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCG  >  NC_000913/4287462‑4287957

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