Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,097,1180CA38.1% 15.0 / 16.8 21intergenic (+289/+447)ycdU/serXuncharacterized protein YcdU/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base A (5/3);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTAATAAGTGAGGAAGTGATAGGAAGTGACCAGATAATACATATATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAC‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTC  >  NC_000913/1096819‑1097412
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                       
ttAATAAGTGAGGAAGTGATAGGAAGTGACCAGATAATACATATATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATCTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAc                                                                                                                                                                                                                                                                                                         <  1:140624/300‑1 (MQ=255)
       gTGAGGAAGTGATAGGAAGTGACCAGATAATACATATATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCCGTTGTCGTc‑cccccgt                                                                                                                                                                                                                                                                                                 >  1:172597/1‑297 (MQ=255)
               gTGATAGGAAGTGACCAGATAATACATATATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAA‑CCCCGGTCGGGGCtt                                                                                                                                                                                                                                                                                         >  1:699419/1‑300 (MQ=255)
                          ggACCAGATAATACATATATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACNGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGTGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGNCTGTCCAACTGGCTGTCGCCAGTTGTCGAC‑CCCCGGTCGGGGCTTCTCACcccccc                                                                                                                                                                                                                                                                              <  2:733520/299‑1 (MQ=255)
                              cAGATAATACATATATTTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGNCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAA‑CCCCGGTCGGGGCTTCTCATCCCCCCGgtg                                                                                                                                                                                                                                                                          <  1:546881/300‑1 (MQ=255)
                                   aataCATATATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGTGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGTGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAC‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGc                                                                                                                                                                                                                                                                     <  1:519718/300‑1 (MQ=255)
                                   aataCATATATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAC‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGc                                                                                                                                                                                                                                                                     <  2:447160/300‑1 (MQ=255)
                                   aagatagacgtGTTCTGTACTCCCTTGCGAAGGTTGATTGGTGACGGAATAACCAGACAGGGGAGGTGCACGATTTCCCTTACCCCTAACAGGCGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAA‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGc                                                                                                                                                                                                                                                                     <  2:686640/290‑1 (MQ=14)
                                      aCATATATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGACAGTTGTCGAC‑CCCCGATCGGGGCTTCTCATCCCCCCGGTGTGTGCAAt                                                                                                                                                                                                                                                                  >  1:527/1‑300 (MQ=255)
                                                                              aCTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAA‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTcttctt                                                                                                                                                                                                                          >  1:461183/1‑300 (MQ=255)
                                                                                 gAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCAAGTTGTCGAA‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCACGTACTTTCGTACGAGCTTTTCTTTaa                                                                                                                                                                                                                       >  2:786075/1‑300 (MQ=255)
                                                                                              gTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCACCCAACCCGCCCTACTATCAGCCCTCTCACTGCGTTCACGATCCGCCCCACCGACTGTCACAAGTTG‑CCACACCCCGGTCGGGGCTTCTCATCCCCCAGGTGTGTGCAATACACGAAAAAAGAGCACGTACTTTAGTACGAGCTCTTCTTTAAATATGGCGGTGAg                                                                                                                                                                                                          >  1:758972/1‑300 (MQ=17)
                                                                                              gTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTAACTCGCCCCACTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGCCGAC‑CCCCGGTCGGGGCTTCTCATCCACCCGGTGTGTGCAATACACGAAAAAAAAGCACGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAg                                                                                                                                                                                                          >  2:462166/1‑300 (MQ=255)
                                                                                                           aTTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAC‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACt                                                                                                                                                                                             <  1:698786/300‑1 (MQ=255)
                                                                                                                         tAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCCGTCCCCCTAGCTACCTCCAATTGCCCAC‑CACAAATCGGGGCGTCTCATACACCCGGTGTGTGCAATACACGAAAAAAAAGCAAGTAGTGTCGTAGGAGCTCTTATGTGAATGTGGAGGTGAGGGGGGGATTGACTCGCTTAGCTCGccc                                                                                                                                                                               >  1:693332/1‑300 (MQ=14)
                                                                                                                              gACGTGGGCGGGGGCACCCGGGGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAC‑CCCCGTTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCgg                                                                                                                                                                          <  1:116522/300‑1 (MQ=255)
                                                                                                                                         gggCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCGGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAA‑CCCCGTTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTc                                                                                                                                                               <  2:455540/300‑1 (MQ=2)
                                                                                                                                            cACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTATGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAC‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTATGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                                                                                                                            >  1:201059/1‑300 (MQ=21)
                                                                                                                                                                                                        atatGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCCACTGGCTATCTCCAATTGTCAAC‑CCAAGGCCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGNGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCCACTGGCTATCGCCAATTGTCAAAccccagccagggattctca                                                                                                >  1:541479/1‑285 (MQ=11)
                                                                                                                                                                                                              cGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTC‑AACTGGCTGTCGCCAGTTGTCGAC‑CCCCGTTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATcccccc                                                                                         <  1:403402/300‑1 (MQ=17)
                                                                                                                                                                                                                                                                                                     cttCGAC‑CCCCGTGCGGTGTTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGTCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGGGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgctt   <  1:190866/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      gTCGAA‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGTGGTGAGGGGGGGATTGACtcgcttc  >  1:283871/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                                                      gTCGAA‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgctta  >  1:465864/1‑299 (MQ=11)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                       
TTAATAAGTGAGGAAGTGATAGGAAGTGACCAGATAATACATATATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAC‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTC  >  NC_000913/1096819‑1097412

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