Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT31.0% 61.9 / 10.6 29N156K (AAC→AAA) yiiQDUF1454 domain‑containing protein YiiQ
Reads supporting (aligned to +/- strand):  ref base G (19/1);  new base T (0/9);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 9.99e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  NC_000913/4111615‑4111883
                                                                                                                                     |                                                                                                                                        
ttttttATTCAACCCTTCAAGCGATTCAGATAGCCCCATCTTAATCGGTTCAACAGCGAAGGCCACCCCTTTTCCGC‑CTTTGCCCGAAAAAAAATAACGCAGTCCCCCTTCTGTCTGGGTGTAATATC‑TTTTTTTTTccc                                                                                                                                  <  2:1008266/137‑1 (MQ=255)
  ttttATCCAACCCTTAAAGCGATTCAGATAGCGCCAGTTTAATCGGTTCAACGGCGAGGGTCACCCCTTTTCCGC‑CGTTGCCCGCGACAACAAAACGCGTTGCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTCcccc                                                                                                                                 <  2:624851/139‑1 (MQ=255)
      aTTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTTAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgccg                                                                                                                             >  2:211898/1‑140 (MQ=255)
                 cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                  >  1:788112/1‑140 (MQ=255)
                 cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                  >  1:554721/1‑140 (MQ=255)
                      gATTCAGATCGCCCCAGCTTATTCGTTACAACAGCGAGGGTCCCCCCTTTTCCGC‑CGTTGTCCGCGAAAACATACCGCAGTGCCCCTTCTGTCTGGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAg                                                                                                             <  2:554721/140‑1 (MQ=255)
                        ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGACAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                           >  2:649307/1‑140 (MQ=255)
                               aGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTg                                                                                                    >  2:240224/1‑140 (MQ=255)
                                  gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 >  2:285661/1‑140 (MQ=255)
                                  cccAGTTTAATCGGTACAACAGGGAGGGTCCCCCCTTTTTCGC‑CGTTGCCCGCGACAACATACCGCAGTGCCCCTTCTGTCTGGGTGTAATACTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 <  2:108149/139‑1 (MQ=255)
                                  cccAGCTAAATCGGTCAAACAGCG‑AGGGCCCCCCCTTTTCCCGCGTTGCCCGCGACAACAAACCGCAGTGCCCTTTCTTTCTGGGTGTAAAACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 <  2:1122224/139‑1 (MQ=255)
                                      gCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                             >  2:46738/1‑140 (MQ=255)
                                      gCTTAATCGGTTCAACAGCGAAGGTAAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                             >  2:197883/1‑140 (MQ=255)
                                      gCTTAATCGGTTCAACAGCGAAGGTAAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                             >  2:197882/1‑140 (MQ=255)
                                       cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                            >  1:108149/1‑140 (MQ=255)
                                                ttCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                   >  2:117788/1‑140 (MQ=255)
                                                 tCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                                                                                  >  1:964576/1‑140 (MQ=255)
                                                    aCAGCGAAGGTTAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                               <  1:211898/140‑1 (MQ=255)
                                                          agggTCACCCCTTTTCCGC‑GGTTGTCCGCAACAACATAACGCAGTGCCCCTTCTGTCTGGGTTTAATAGTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:1123490/138‑1 (MQ=255)
                                                          aaGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         >  2:584499/1‑140 (MQ=255)
                                                          aaGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         >  1:999310/1‑140 (MQ=255)
                                                          aaGGTCACCCCTTTTCCGC‑CGTTGCCCGCAACAACAAACCGCAGTGCCCCTTCTGTCTCGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:999310/140‑1 (MQ=255)
                                                          aaGGCCACCCCTTTTCCGC‑CGTTGCCCGCGAAAACATACCGCAGTGCCCCTTCTGTCTCGGTGTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:208918/140‑1 (MQ=255)
                                                              ccacccctTTTTCGC‑CGTTGCCCGCGACAACATACCGCAGTGCCCTTTCTGTCTGGGTGTAATACTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt                                                                     <  2:1031146/136‑1 (MQ=255)
                                                                            c‑cGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTCATAACGTTTGTTTTTCCCAGCCGTTAGTAGCTACTGTAGTTTTTTCTGGCTTTGTGTTTTGGTCTTTAATGGGGTGCGTGTGCGGATCACCgt                                                       >  2:811058/1‑139 (MQ=255)
                                                                                ttGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGc                                                    <  1:670534/140‑1 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:1008794/1‑140 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:213216/1‑140 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:104827/1‑140 (MQ=255)
                                                                                                            cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                        >  1:335040/1‑140 (MQ=255)
                                                                                                                       gtgtAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc             <  2:947393/140‑1 (MQ=255)
                                                                                                                        tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTCTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct            >  2:849451/1‑140 (MQ=255)
                                                                                                                        tataatatCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct            <  2:788112/138‑1 (MQ=255)
                                                                                                                           aatagttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc         <  2:335040/133‑1 (MQ=255)
                                                                                                                           aagatggtttgtttTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc         <  1:811058/134‑1 (MQ=255)
                                                                                                                                ttttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    <  1:881932/138‑1 (MQ=255)
                                                                                                                                gttttgtttTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    <  1:240224/138‑1 (MQ=255)
                                                                                                                                  tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  1:179696/1‑140 (MQ=255)
                                                                                                                                  tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  2:871470/1‑140 (MQ=255)
                                                                                                                                  tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  2:314850/1‑140 (MQ=255)
                                                                                                                                  tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  1:897368/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                        
TCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGC‑CGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  NC_000913/4111615‑4111883

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