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

TTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111609‑4111885
                                                                                                                                           |                                                                                                                                          
tttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                                                                           <  2:369136/140‑1 (MQ=255)
   ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                                                                                        >  1:876676/1‑140 (MQ=255)
                       cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                    >  1:803557/1‑140 (MQ=255)
                                        gcCAGCTTAATCGGTTCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   >  2:159683/1‑140 (MQ=255)
                                        gcCAGCTTAATCGGTTCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   >  1:895035/1‑140 (MQ=255)
                                        gcCAGCTTAATCGGTTCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   >  2:665302/1‑140 (MQ=255)
                                         ccAGTTTAATCGGTTCAACAGCGAGGGT‑CCCCCCTTTTCCGCCGTTGCCCGCGAAAACATACCGCAGTCCCCCTTCTGTCTGGGTGTAATATCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCtt                                                                                                  <  2:1235844/140‑1 (MQ=255)
                                          cAGCTTAATCGTTCCAACAGCGAAGGC‑CCCCCCTTTTCCGCCGTTGCCCGCAAAAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATATTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                                                                                 <  2:1027838/140‑1 (MQ=255)
                                               tAATCGTTCCAACAGC‑AAGGTCCACCCCTTTTTCGCCTTTGCCCGCAAAAACATACCGCAGTCCCCTTTCTGTCTGGGTATAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtt                                                                                            <  2:803557/140‑1 (MQ=255)
                                                aaTCGGTCCAACAGGGAAGCT‑CCCCCCTTTTCCGCCTTTGCCCGCAAAAACATACCGACGTCCACCTTCTGGCTGGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                                                                                           <  2:1455371/140‑1 (MQ=255)
                                                    gtttaaaCAGCGAAGCC‑CACCCCTTTTTCGCTTTTGCCCGCGAAAAAATAACGCAGTGCCCTTTCTGTCGCGTTGAAAAATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:1369914/135‑1 (MQ=255)
                                                    gttacAACAGCGAGGGT‑CCCCCCTTTTCCGCCGTTGCCCGCGAAAACATAACGCTGTCCCCCTTCTGTCTGGGTGTAATATTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:1497436/136‑1 (MQ=255)
                                                    ggTTAAACAGCGAAGGT‑CACCCCTTTTTCGCCGTTGCCCGCAAAAACATACCGCAGTGCCCCTTCTGTCGGGGTTTAATATTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:666301/140‑1 (MQ=255)
                                                      ttCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     >  2:1488813/1‑140 (MQ=255)
                                                      ttCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     >  1:1534747/1‑140 (MQ=255)
                                                      ttCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     >  1:1497436/1‑140 (MQ=255)
                                                      taaaaCAGAGAAGGC‑CACCCCTTTTCCGCCGTTGCCCGCAACAACATACCGCAGTCCCCTTTCTGTCTGGGTGAAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     <  2:758484/137‑1 (MQ=255)
                                                       tCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                                                                                    >  1:228160/1‑140 (MQ=255)
                                                           cAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                                                                                >  1:69581/1‑140 (MQ=255)
                                                               gAAGGT‑CCCCCCTTTTTCGCCTTTGCCCGCAAAAACATACCGCAGTCCCCCTTCTGTCTGGGTTTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg                                                                            <  2:895035/140‑1 (MQ=255)
                                                                aaGGT‑CATCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAtctgtc                                                                           >  1:1069901/1‑135 (MQ=255)
                                                                  ggT‑CATCCCTTTTGGGCGGTTCTCCGCAATAATATAACGCCGTGCTCCTCGTGTGTGGGTATAAAATTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATAGGGTGa                                                                         <  2:1114652/140‑1 (MQ=255)
                                                                  ggT‑CACCCCTTTTCCGCGGTTGTCCGCAACAACATACCGCAGTCCCCCTTCTGTCTCGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa                                                                         <  2:15861/140‑1 (MQ=255)
                                                                      cAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtg                                                                      >  2:336576/1‑140 (MQ=255)
                                                                              tttCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                              >  1:15861/1‑140 (MQ=255)
                                                                                  gCCGTTGCCCGCGACAACATACCGCAGTGCCCCTTCTGTCTCGGTGTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                          <  1:923389/140‑1 (MQ=255)
                                                                                     gtggcccGCAAAAAAATAACGCAGTCCCCTTTCTTTCGGGGTTTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                       <  2:1561466/135‑1 (MQ=255)
                                                                                     gTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                       >  1:80470/1‑140 (MQ=255)
                                                                                                 aacaTAACGCAGTCCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCt                                           <  1:555236/140‑1 (MQ=255)
                                                                                                             tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:120718/1‑140 (MQ=255)
                                                                                                             tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:263041/1‑140 (MQ=255)
                                                                                                             tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACAGCTGCCATGTATTCCTGAgctttcgctt                               >  1:670699/1‑140 (MQ=255)
                                                                                                              gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                              >  1:351573/1‑140 (MQ=255)
                                                                                                               cACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctttc                             >  2:154248/1‑140 (MQ=255)
                                                                                                                       gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                     >  1:656250/1‑140 (MQ=255)
                                                                                                                              tgtAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  2:263041/140‑1 (MQ=255)
                                                                                                                              tgtAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  2:1046950/140‑1 (MQ=255)
                                                                                                                              tgtAATAACGTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  1:336576/140‑1 (MQ=255)
                                                                                                                                  ttacttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          <  2:228160/134‑1 (MQ=255)
                                                                                                                                      ttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  2:351573/138‑1 (MQ=255)
                                                                                                                                      ttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  2:80470/138‑1 (MQ=255)
                                                                                                                                        tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:309325/1‑140 (MQ=255)
                                                                                                                                          tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                <  1:896499/126‑1 (MQ=255)
                                                                                                                                          tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                >  2:896499/1‑126 (MQ=255)
                                                                                                                                          tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  <  2:656250/140‑1 (MQ=255)
                                                                                                                                          tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTAGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:426369/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                          
TTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGT‑CAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111609‑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: ‑