Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913550,1880GT25.0% 87.9 / 21.8 36A188S (GCG→TCG) ylbFDUF2877 domain‑containing protein YlbF
Reads supporting (aligned to +/- strand):  ref base G (11/16);  new base T (5/4);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACTCCTGCACTTTGTTCATGCT  >  NC_000913/550049‑550319
                                                                                                                                           |                                                                                                                                   
tCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATg                                                                                                                                     <  1:290800/140‑1 (MQ=255)
tCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATg                                                                                                                                     <  1:407168/140‑1 (MQ=255)
    tGTTTCCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATgcgcg                                                                                                                                 <  1:491474/140‑1 (MQ=255)
               cgcgCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGgtcgtc                                                                                                                      <  1:574976/140‑1 (MQ=255)
                gcgcTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTcc                                                                                                                     >  1:739983/1‑140 (MQ=255)
                         ggCGTTAAGACCGACTGGCGGCACCGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGc                                                                                                            <  2:739983/140‑1 (MQ=255)
                            gTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTg                                                                                                         <  2:9469/140‑1 (MQ=255)
                                aGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTcc                                                                                                     >  1:36959/1‑140 (MQ=255)
                                      aCTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAAt                                                                                               >  2:40640/1‑140 (MQ=255)
                                      aCTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAAt                                                                                               >  2:679536/1‑140 (MQ=255)
                                           cggcACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATTCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACa                                                                                                 <  2:300676/133‑1 (MQ=255)
                                           cggcACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATTCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACa                                                                                                 >  1:300676/1‑133 (MQ=255)
                                            ggcACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                                                                         <  1:234812/140‑1 (MQ=255)
                                            ggcACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                                                                         <  1:679536/140‑1 (MQ=255)
                                            ggcACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                                                                         <  2:468328/140‑1 (MQ=255)
                                                   gctggGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATTCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTa                                                                                  >  1:70279/1‑140 (MQ=255)
                                                           aaGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATTCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCg                                                                          >  1:312651/1‑140 (MQ=255)
                                                            aGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATTCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCgg                                                                         <  1:319480/140‑1 (MQ=255)
                                                                  cGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCg                                                                   >  2:467638/1‑140 (MQ=255)
                                                                    gATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATTCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTc                                                                 <  2:312651/140‑1 (MQ=255)
                                                                              gAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTAc                                                       >  2:353431/1‑140 (MQ=255)
                                                                               aGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCa                                                                <  2:776601/130‑1 (MQ=255)
                                                                               aGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCa                                                                >  1:776601/1‑130 (MQ=255)
                                                                               aGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACg                                                      >  2:430114/1‑140 (MQ=255)
                                                                                   atgatgACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTAt                                                  >  2:146009/1‑140 (MQ=255)
                                                                                      atgaCACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGcc                                               >  2:749458/1‑140 (MQ=255)
                                                                                       tgaCACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCcg                                              <  1:749458/140‑1 (MQ=255)
                                                                                           acGCTGAGCGGAACGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATTCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACg                                                      >  1:628741/1‑128 (MQ=255)
                                                                                           acGCTGAGCGGAACGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATTCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACg                                                      <  2:628741/128‑1 (MQ=255)
                                                                                              cTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAgg                                       <  1:146009/140‑1 (MQ=255)
                                                                                                            gctgGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTtcgcctcgc                         <  1:40640/140‑1 (MQ=255)
                                                                                                                gcggcCTGGTATTATGGCGCTTTAGATTCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACt                     >  1:567254/1‑140 (MQ=255)
                                                                                                                 cggcCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACTc                    <  2:260378/140‑1 (MQ=255)
                                                                                                                       ggTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACTCCTGCAc              <  1:353431/140‑1 (MQ=255)
                                                                                                                              tGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACTCCTGCACTTTGTTc       <  2:36959/140‑1 (MQ=255)
                                                                                                                              tGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACTCCTGCACTTTGTTc       <  1:618328/140‑1 (MQ=255)
                                                                                                                              tGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACTCCTGCACTTTGTTc       <  2:98321/140‑1 (MQ=255)
                                                                                                                                   cTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACTCCTGCACTTTGTTCATGct  >  2:183061/1‑140 (MQ=255)
                                                                                                                                   cTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACTCCTGCACTTTGTTCATGct  <  1:467638/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
TCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTATTTACGTTATGCCGCGCAAGGATATTTCGCCTCGCCACTCCTGCACTTTGTTCATGCT  >  NC_000913/550049‑550319

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