Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,988,0510AT21.1% 72.1 / 15.8 38intergenic (‑92/‑171)azuC/yecRacid‑inducible small membrane‑associated protein/lipoprotein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (16/14);  new base T (4/4);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTT  >  NC_000913/1987906‑1988198
                                                                                                                                                 |                                                                                                                                                   
gAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAAt                                                                                                                                                   <  8:85976/148‑1 (MQ=255)
   cGTCTGGCAATAGTTATTGAACATACGTTTCAGGATGTTGCGCAGTTTCATCGCGGCACTCCGACCCTTGGTGATACAGGTGTTAGTGTCTTTGCGCTTATAATATGACACCCATCACAAAAATCAATCTTTATGTGATACAAATcaca                                                                                                                                               >  6:5674/1‑149 (MQ=255)
      cTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATa                                                                                                                                             >  6:72866/1‑148 (MQ=255)
         gCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaa                                                                                                                                           <  3:134845/147‑1 (MQ=255)
         gCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATa                                                                                                                                             >  7:105470/1‑145 (MQ=255)
          cAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaaga                                                                                                                                         <  8:105470/148‑3 (MQ=255)
           aaTAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaagatt                                                                                                                                       <  4:22055/149‑5 (MQ=255)
           aaTAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaaga                                                                                                                                         <  3:114542/147‑3 (MQ=255)
           aaTAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaaga                                                                                                                                         >  4:114542/1‑145 (MQ=255)
                      aaCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGt                                                                                                                             >  4:5444/1‑148 (MQ=255)
                       aCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTta                                                                                                                           <  5:72866/149‑1 (MQ=255)
                                 cAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGa                                                                                                                  <  2:7277/148‑1 (MQ=255)
                                    gATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATaa                                                                                                               >  4:45052/1‑148 (MQ=255)
                                    gATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAAt                                                                                                              >  1:13931/1‑149 (MQ=255)
                                        ttGCGCATTTTCATCGCGTCACTCAGACCATTTGTTATACATGTTTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCaa                                                                                                           <  7:116539/148‑1 (MQ=255)
                                              aGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAg                                                                                                    <  2:13931/149‑1 (MQ=255)
                                                tttCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATca                                                                                                                                                 >  8:80854/1‑102 (MQ=255)
                                                tttCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATca                                                                                                                                                 <  7:80854/102‑1 (MQ=255)
                                                                cGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                     >  2:43406/1‑146 (MQ=255)
                                                                cGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                     <  1:43406/146‑1 (MQ=255)
                                                                                  gtgtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGaaa                                                                 >  1:130888/1‑148 (MQ=255)
                                                                                   tgtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGaaa                                                                 >  2:88019/1‑147 (MQ=255)
                                                                                           tCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGc                                                       >  5:29407/1‑149 (MQ=255)
                                                                                           tCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGc                                                       >  5:31733/1‑149 (MQ=255)
                                                                                                gcgcTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATa                                                                                                                                             <  8:132442/58‑1 (MQ=255)
                                                                                                gcgcTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATa                                                                                                                                             >  7:132442/1‑58 (MQ=255)
                                                                                                 cgcTTATAATATGACACCCATCACAAAAATCAATCTTTATGTGCTACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAACCAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGc                                                  <  5:5674/148‑1 (MQ=255)
                                                                                                       taataTGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCttt                                            <  3:5444/148‑1 (MQ=255)
                                                                                                         atatGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTcc                                         <  1:88019/149‑1 (MQ=255)
                                                                                                            tGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGc                                      >  2:171834/1‑149 (MQ=255)
                                                                                                              aCAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCa                                    >  8:128926/1‑149 (MQ=255)
                                                                                                                    aTCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCg                              <  6:31733/149‑1 (MQ=255)
                                                                                                                    aTCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCg                              <  6:29407/149‑1 (MQ=255)
                                                                                                                        caAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgcag                           >  5:21485/1‑148 (MQ=255)
                                                                                                                                ggtatttATGTGATTTGTATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAAc                  >  1:147518/5‑149 (MQ=255)
                                                                                                                                 gtatttATGTGATTTGTATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAAcg                 >  4:98872/4‑149 (MQ=255)
                                                                                                                                  tcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAATAACAGCACCCCTATTTCTGACACCTAAAGATGGCAAGAAATAGAGCCTGCCAGGCGTCTTTTCAGTCCATTGTCGCAGCACTGTAAcgc                >  7:42886/1‑149 (MQ=255)
                                                                                                                                             tACAAATCACATAAATACCACATTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCATCACTGTAACGCGTAAAATAGTg     <  8:42886/149‑1 (MQ=255)
                                                                                                                                                aaaTCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCtt  >  6:99738/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTT  >  NC_000913/1987906‑1988198

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