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

ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  NC_000913/4111617‑4111883
                                                                                                                                   |                                                                                                                                        
aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTNANCGNTTCAACAGCGAGGGTCANCCNTTTTCCGCCG‑TTGCCCGCGAAAACATACCGCAGTGCCCCTTCTGTCTGGGTGAAATAACTTTTTTTTTTcccc                                                                                                                                 <  2:1192457/140‑1 (MQ=255)
aTTTATTAAAACCTACAAGCGATTAAGACAGCGCCAGTTTAATGGGTCCAACAGCGAAGGCCACCCCTTTTCCCCCG‑TTGTCCGCAACAACATACCGCAGTGCCCCTTCTGTCTCGGTGTAATAGCGTTTTTTTTCcccc                                                                                                                                 <  2:304413/140‑1 (MQ=255)
               cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                  >  2:617210/1‑140 (MQ=255)
                  gCGATTAAGATGGCCCCAGCTTAATCGGTTAAACGGCGAAGGCCACCCCTTTTCCGCCG‑TTGCCCGCAAAAACATACCGCAGTGCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGc                                                                                                               <  2:74545/140‑1 (MQ=255)
                     aTTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGt                                                                                                            >  1:1358734/1‑140 (MQ=255)
                      ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                           >  1:565153/1‑140 (MQ=255)
                          gacaGCGCCGGCTAATTGGGTCAAACAGCG‑AGGTCCCCCCCTTTTCCCCGTTTGCCCGCAACAACATACCGCAGTCCACTTTCTGTCTGGGTATAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                                                                       <  2:565153/137‑1 (MQ=255)
                                gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 >  1:1295701/1‑140 (MQ=255)
                                   aGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtg                                                                                              >  1:383779/1‑140 (MQ=255)
                                            gtttCAACAGCGAGGGCCACCCCTTTTCCGCCT‑TTGCCCGCGAAAACATACCGCAGTCCCCTTTCTGTCGGGGTTTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                     <  2:1184866/138‑1 (MQ=255)
                                              ttCAACAGCGAAGGTCAGCCCCTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                   >  1:1016490/1‑140 (MQ=255)
                                                  aCAGCAAAGGCCCCCCCCTTTTCGCCG‑TTGCCCGCAAAAACATAACGCAGTCCCCCTTCTGTCTGGGTTAAATATCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                               <  2:1214285/140‑1 (MQ=255)
                                                      ggagggTCACCCCTTTTCCGCCG‑TTGCCCGCAACAACATACCGCAGTGCCCCTTCTGTCTGGGTGTAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                                                                           <  2:345238/136‑1 (MQ=255)
                                                        agggcccgCCCTTTTCCGCCG‑TTGCCCGCAAAAACAAAACGCAGTCCCCCTTCTGTCTGGGTTTAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:382600/133‑1 (MQ=255)
                                                        aaggcgaGCCCTTTTTCGCCT‑TTGCGCGCGAAAAAATATCCCAGGCCACCTTCTGTTGCGGTTTAATAGCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:1295701/134‑1 (MQ=255)
                                                        aaGGTCACCCCTTTTCCGCCG‑TTGCCCGCAAAAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:632836/140‑1 (MQ=255)
                                                           gTCAGCCCCTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAg                                                                      >  1:566210/1‑140 (MQ=255)
                                                             cccccctTTTTCGCCG‑TTGCCCGCGACAACATAACGCAGTCCCCCTTCTGTCTCGGTTTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtg                                                                    <  2:1358734/137‑1 (MQ=255)
                                                              cgCCCTTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCCCTTTCTGTCTGGGTTTAATAGCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtgt                                                                   <  2:1016490/139‑1 (MQ=255)
                                                                ccccTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc                                                                 >  1:919442/1‑140 (MQ=255)
                                                                     tttCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                            >  2:967941/1‑140 (MQ=255)
                                                                       tCGCCT‑TTGCCCGCAACAACATACCGCATTCCTCCTTCTGTCTGGTTTTAATATTTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAc                                                          <  2:566210/140‑1 (MQ=255)
                                                                          ccG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGc                                                       >  1:572968/1‑140 (MQ=255)
                                                                             tttGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATAACTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                     <  2:919442/139‑1 (MQ=255)
                                                                                  ccGCAACAACATAGCGAAGTGCCCTTTCTTTGTGGTTGTAAAAATGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg                                                <  2:791701/140‑1 (MQ=255)
                                                                                   cgcgACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGt                                               >  2:420825/1‑140 (MQ=255)
                                                                                        caacaTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTcc                                          >  2:1418045/1‑140 (MQ=255)
                                                                                           caTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa                                       >  1:439849/1‑140 (MQ=255)
                                                                                                   aGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgc                               >  1:478485/1‑140 (MQ=255)
                                                                                                     tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:513555/1‑140 (MQ=255)
                                                                                                     tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:479151/1‑140 (MQ=255)
                                                                                                     tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  2:760585/1‑140 (MQ=255)
                                                                                                          cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                        >  1:1097667/1‑140 (MQ=255)
                                                                                                                    ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttctgtct              >  1:357762/1‑134 (MQ=255)
                                                                                                                      tgtAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct            <  1:1029260/140‑1 (MQ=255)
                                                                                                                      tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct            >  1:47219/1‑140 (MQ=255)
                                                                                                                      taaaatatCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct            <  2:383779/137‑1 (MQ=255)
                                                                                                                           taaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGccc       >  2:1365089/1‑140 (MQ=255)
                                                                                                                              ttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    <  2:479151/138‑1 (MQ=255)
                                                                                                                                ttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc        <  1:799302/134‑1 (MQ=255)
                                                                                                                                tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc        >  2:799302/1‑134 (MQ=255)
                                                                                                                                tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  1:1065778/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                        
ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCG‑TTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  NC_000913/4111617‑4111883

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: ‑