Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7190TA28.1% 58.8 / 12.8 32intergenic (+52/‑45)ytcB/nrfAprotein YtcB/cytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base T (18/5);  new base A (0/9);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 7.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCC  >  NC_000913/4287495‑4287967
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
gCAGGGTAATTTTTGAAGGTCAACAAAAAAAGTTAATAAATTCTTTGAGGAAAATGAAGTTATGCATGCTGTTGAAAAAGAGAAAGATACGGACTAACTCAAAAGTGGTATTTTACATGCATTTACAATTGATTAAAGAAAACATTCACGGTGTGGTTATTTGTTACACATGGGGGCGAGAAATGTCATGC‑CAGTGTGGGTGCGGTTACTGTCGTAAAAAAAAATAAAAACCGCCATTGCAACAATGgcgc                                                                                                                                                                                                                                <  2:255295/251‑1 (MQ=255)
    ggTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATGGATAAAAGACAACATTCACAGGGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAtt                                                                                                                                                                                                                            <  1:197959/251‑1 (MQ=255)
          ttttGAAGGTAAGCAAAAAAAGTTAATAAATTCTTTGAGGAACATGCATTTAGGAATGTTGTAGAAAAAGAGGAAGATACTCATTAATTCTAAGGTGGTATTTAACGTGCTTTTAAAATTGATAAAAGACAAAATTCACAGTGTGGTTAGTTGTTACCAATGGGGGAGTACAGTGTCATTA‑CGTTGGGGGGGCGTTTAGTGTCGTAAAAAACAATGAAAACCGCCATTGCAACAATGGCGCAATTCGGAt                                                                                                                                                                                                                       <  2:144761/250‑1 (MQ=255)
                    cagaaaaaaaaGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTAAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGGGGGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGA‑CAGTGTAGGGGCGGTTACTGTCGTGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                                                                                                                            <  2:75083/244‑1 (MQ=255)
                     aGCAACAAAAGTTGATTAATTTTTTGAGGAACATGCAGTTATGCATGCTGTTAAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                                                                                                                                           <  1:141847/251‑1 (MQ=255)
                       caacaaAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc                                                                                                                                                                                                         >  1:190942/1‑251 (MQ=255)
                        aacaaAAGTTGATTAATTCTTGGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGGGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa                                                                                                                                                                                                        <  1:225743/251‑1 (MQ=255)
                             aaGTTGATTAATTCTTTGAGGACCATGCAGTTATGCATGCTGTTGAAAAAGTGGAAGATACTGCATAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTGAGAGTGTGGTTATGTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGGAGGTGCGGTTACTGTCGTGAAAAACAATACAAACCGCCGTTGACACAATGGCGCAATTCGGATGACGCCCCGATGACCAGGa                                                                                                                                                                                                    >  2:119231/1‑250 (MQ=255)
                                        ttCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTCCTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGAt                                                                                                                                                                                                   <  1:26057/240‑1 (MQ=255)
                                        ttCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGAt                                                                                                                                                                                                   >  2:26057/1‑240 (MQ=255)
                                                        aatttaggcaggctttaaaaaaaGAGAAAGATACTAACTACCTAAAAGGTGGTTTTTTACTTGCACTTACAATTAATAAAAGACAACATCCACAGTGTGTTTTTTTGTTACACATGGGGGGGAGAAATTTCAGGA‑CGGGGAAGGTGGGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAg                                                                                                                                                                        <  2:114395/233‑1 (MQ=255)
                                                               gcatgcTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                                                                                                                                 >  1:212323/1‑251 (MQ=255)
                                                                        tGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt                                                                                                                                                        >  1:287899/1‑251 (MQ=255)
                                                                          aaaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                                                                                                                                      >  2:175159/1‑251 (MQ=255)
                                                                                    aGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACt                                                                                                                                            >  1:366615/1‑251 (MQ=255)
                                                                                     gATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACtt                                                                                                                                           >  1:297702/1‑251 (MQ=255)
                                                                                        attgaaaaactaaaagtttgTATTTTAAATGCATTTAAAATTGTTAAAAAAAAACATTCAGAGTGTGTTTATTTGTTACAAAGAGGGGAGAACAATGTAATGA‑CGGTGTGGGTGCGTTTACTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAAAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                                                                                                                                         <  2:68869/231‑1 (MQ=255)
                                                                                                                  aCATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTccc                                                                                                              >  1:295169/1‑251 (MQ=255)
                                                                                                                                ttGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATGGGGGCGAGAAATGTCAGGA‑CAGTGTAGGGGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                                                                                                                                          <  2:191202/129‑1 (MQ=255)
                                                                                                                                ttGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                                                                                                                                          >  1:191202/1‑129 (MQ=255)
                                                                                                                                   aTTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACtgt                                                                                             >  1:43870/1‑251 (MQ=255)
                                                                                                                                       aaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGa                                                                                          >  1:3560/1‑250 (MQ=255)
                                                                                                                                                    caGTGTGGTTTTTTGTAACGCATAGGGGCGAAAAATTGAATGA‑CGTTGTAGGTGCGTTTCCTTTCTTAAAAAACAAAAAAAACCGCCATTGCAACAAGGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGaaa                                                                            <  1:239459/251‑1 (MQ=255)
                                                                                                                                                        gtgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCtt                                                                        >  1:252765/1‑251 (MQ=255)
                                                                                                                                                                         aatgggggaaagaaTGTCAAGAACAGGGTAGGTGCGTTTACTGTGGAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATccc                                                        <  1:72078/239‑1 (MQ=255)
                                                                                                                                                                                                  gtgtGGGTGGGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCtt                                                                        <  1:184163/210‑1 (MQ=255)
                                                                                                                                                                                                  gtgtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCtt                                                                        >  2:184163/1‑210 (MQ=255)
                                                                                                                                                                                                     taggggcgtttactttctaaaaaacCAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCAcc                             <  1:399097/230‑1 (MQ=255)
                                                                                                                                                                                                            gtttactgtggaaaaaaaaaaaaaaACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACa                     <  2:43870/238‑1 (MQ=255)
                                                                                                                                                                                                                       tGAAAAACAATAAAAACCGCCATTGCAACAATGGAGCAATTCCGATGAAGCCCCTATGACAAGGATACACATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCgtg          >  1:374226/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaaCAATAAAAACCGCCATTGCAACAATGGAGCCATTCGGATGAAGCCCCTATGACACGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACCAACGGCTGCTCCCGCACAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTg       >  1:183798/1‑250 (MQ=255)
                                                                                                                                                                                                                               aaTAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGcc  >  1:27403/1‑251 (MQ=255)
                                                                                                                                                                                                                               aaTAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGcc  >  1:244751/1‑251 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
GCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGA‑CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCC  >  NC_000913/4287495‑4287967

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