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

CATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111616‑4111885
                                                                                                                                     |                                                                                                                                          
tttttATTAAACCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCCCCCCTTTTC‑CGCCGTTGCCCGCGACAACATAA‑CGCAGTGCCCCTTCTGTCTCGGTGTAATAAC‑TTTTTTTTTcccc                                                                                                                                   <  2:111605/138‑1 (MQ=255)
              ttCAAGCTATTCAGATAGCCCCAGCTTAATCGTTTCAACAGCGAAGGCCAGCCCCTTTC‑CGCCGTTGCCCGCGAAAACAAAA‑CGCAGTGCCCCTTCTGTCTGGGTGTAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCg                                                                                                                      <  2:88410/140‑1 (MQ=255)
                aaaGCGATTCAGATCGCCCCACTTTAATCGGTTCAACAGCGAAGGTCAGCCCTTTTC‑CGCCGTTGCCCGCGAAAACATAC‑CGCAGTGCCCCTTCTTTCTGGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGAc                                                                                                                    <  2:250766/139‑1 (MQ=255)
                       ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                             >  2:470378/1‑140 (MQ=255)
                            atggcgcCACCTAAATCGGTCAAACAGAGAGGGTCACCCCTTTTC‑CGCTGTTGTCCGCGAAAAAAAAC‑CGCAGTGCCCTTTCTGTCTGGGTTTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                                                                                        <  2:202915/137‑1 (MQ=255)
                            aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTCTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGACGTTAGTAGCGACTGCAGTTTTTTc                                                                                                        <  1:921744/140‑1 (MQ=255)
                              ggcgcCAGTTTAATCGTTACAACAGCAAAGGTCCCCCCCTTTT‑CGCCGTTGCCCGCAACAACATAC‑CGCAGTGCCCCTTCTGTCTGGGTTAAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTg                                                                                                      <  2:316929/139‑1 (MQ=255)
                               gcgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                                                                                     >  2:345356/1‑140 (MQ=255)
                                         aattggTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                                                                                           >  1:92718/5‑140 (MQ=255)
                                             ggTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       >  1:204646/1‑140 (MQ=255)
                                               ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     >  1:844322/1‑140 (MQ=255)
                                               ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     >  2:666796/1‑140 (MQ=255)
                                                   aCAGCGAAGGTCACCCCTTTTC‑CGCGTTTGCCCGCGACAAAATAA‑CGCAGTCCCCCTTCTGTCTGGGTGTAATATTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                                 <  2:1095894/140‑1 (MQ=255)
                                                         aaGGTCACCCCTTTTT‑CGCCGTTGCCCGCGACAACATAA‑CGCATTCCCCCTTCTGTCTGGGTTTAATACCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:741333/140‑1 (MQ=255)
                                                         aaGGCCACCCCTTTTT‑CCCGTTTGCCCGCAACAACATAA‑CGCAGTCCCCTTTCTGTCTGGGTTTAATATTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:857836/140‑1 (MQ=255)
                                                          acggcccccccTTTTCCCCCTTTGTCCGC‑AAAACAAAAGCGCGGTCCCCTTTCTTTGTGGGTGTAAAACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:1122495/133‑1 (MQ=255)
                                                                 ccccTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc                                                                   >  1:1122495/1‑140 (MQ=255)
                                                                           gCCGTTGTCCGCAACAACATAA‑CGCAGTGCCCTTTCTGTCTCGGTGTAATAGGGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                          <  1:470378/140‑1 (MQ=255)
                                                                             gtttGCCCGCGACAACATAA‑CGCGGTGCCCTTTCTGTCTGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                        <  2:741486/138‑1 (MQ=255)
                                                                              gTTGTCCGCGACAACATAA‑CGCAGTGCCCCTTCTGTCTGGGTGTAATATTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                       <  2:3637/140‑1 (MQ=255)
                                                                              gTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                       >  1:3637/1‑140 (MQ=255)
                                                                                tGTCCGCGAAAACATAA‑CGCAGTCCCCTTTCTGTCTCGGTATAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                     <  1:143576/140‑1 (MQ=255)
                                                                                         caacaTAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTcc                                            >  1:1163693/1‑140 (MQ=255)
                                                                                                  cGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCATTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgcttt                                    >  1:522963/1‑140 (MQ=255)
                                                                                                   gCAGTCCCCTTTCTGTCTGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttc                                   <  1:345356/140‑1 (MQ=255)
                                                                                                    cAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcg                                  >  2:372020/1‑140 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:917462/1‑140 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  2:685988/1‑140 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:850633/1‑140 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  2:1131214/1‑140 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTACCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:601744/1‑140 (MQ=255)
                                                                                                            cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                          >  1:372013/1‑140 (MQ=255)
                                                                                                             ttCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCgg                         >  1:238592/1‑140 (MQ=255)
                                                                                                             ttCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCgg                         >  1:1147724/1‑140 (MQ=255)
                                                                                                              tCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                        >  1:590222/1‑140 (MQ=255)
                                                                                                                 gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                     >  1:760201/1‑140 (MQ=255)
                                                                                                                  tctcGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCtttt                    >  1:345640/1‑140 (MQ=255)
                                                                                                                      ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                >  2:1135486/1‑140 (MQ=255)
                                                                                                                      ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                >  2:1057631/1‑140 (MQ=255)
                                                                                                                        tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              >  1:601275/1‑140 (MQ=255)
                                                                                                                                ttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  2:601275/138‑1 (MQ=255)
                                                                                                                                ttttttttttCCCCCGCCATTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  2:522963/138‑1 (MQ=255)
                                                                                                                                  tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:215348/1‑140 (MQ=255)
                                                                                                                                   ttgtttTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGAt   <  2:352421/140‑1 (MQ=255)
                                                                                                                                    tgtttTTCCACGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  1:880927/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                          
CATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111616‑4111885

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