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

TGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  NC_000913/4111613‑4111883
                                                                                                                                        |                                                                                                                                        
tGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                                                                                     >  1:960258/1‑140 (MQ=255)
          tCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGtt                                                                                                                           >  2:717975/1‑140 (MQ=255)
                  tCAAGCGATTCAGATACCCCCACTTTAATCGGTTCAACAGCGAAGGCCCCCCCTTTTC‑C‑GCCGTTGCCCGCAACAACATAACGCGGTCCCCCTTCTGTCGGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGa                                                                                                                   <  2:943789/140‑1 (MQ=255)
                   cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                  >  2:883984/1‑140 (MQ=255)
                   cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCCTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                  >  1:40662/1‑140 (MQ=255)
                              gATAGCGCCAGTTTAATCGTTTAAACAGCGAGGGTCACCCCTTTTC‑C‑GCCTTTGCCCGCAAAAACATACCGCAGTCCCCCTTCTGTCTGGGTGTAATAAGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                                                                       <  2:286167/140‑1 (MQ=255)
                                    gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 >  1:286167/1‑140 (MQ=255)
                                    gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 >  1:732589/1‑140 (MQ=255)
                                    cccattttAATCGGTTAAACAGCGAAGGCCCCCCCCTTTT‑C‑GCCGTTGCCCGCGACAACATAACGCAGTCCCCTTTCTGTCTCGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 <  2:807944/134‑1 (MQ=255)
                                        gCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGCGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTTTGGCTTtgt                                                                                             >  1:374634/1‑140 (MQ=255)
                                                      aCAGCGAAGGTCAGCCCTTTTT‑C‑GCCTTTGCCCGCGACAACATAACGGAGTGCCCCTTCTGTCTGGTTGTAATAGCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                               <  2:616472/140‑1 (MQ=255)
                                                            aaGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         >  1:385258/1‑140 (MQ=255)
                                                            aaGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         >  1:447912/1‑140 (MQ=255)
                                                            aaGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         >  1:559513/1‑140 (MQ=255)
                                                            aaGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         >  2:541710/1‑140 (MQ=255)
                                                            aaGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         >  1:830200/1‑140 (MQ=255)
                                                             ggtgtcgCCTCTTTT‑CGGGCGTCGTCCGCAAAACCATAACGTGGTCCTCCTTCTCCCGGGGTGAAATGCGGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:960258/134‑1 (MQ=255)
                                                             aaggccgCCCCTTTTCC‑GCGGTTTCCCGCAACAAAATACCGCGGTCCCCCTTCTGTCTGGGTGTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:447912/134‑1 (MQ=255)
                                                                    ccccTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc                                                                 >  1:300950/1‑140 (MQ=255)
                                                                            cc‑GCTGTGGTCCGCAAAAACATACCGCTGTCCCCTTTCTGTCTGGGTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAc                                                          <  2:72044/139‑1 (MQ=255)
                                                                                 cGTTGTCCGCGACAACATAACGCAGTCCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                      <  1:338221/140‑1 (MQ=255)
                                                                                        cgcgACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGt                                               >  1:268190/1‑140 (MQ=255)
                                                                                                     cGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgcttt                                  >  1:641251/1‑140 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  2:607989/1‑140 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  2:997736/1‑140 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:956012/1‑140 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:97843/1‑140 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:86164/1‑140 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:942546/1‑140 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTCGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:225079/1‑140 (MQ=255)
                                                                                                               tttctttctGGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                        <  2:300950/139‑1 (MQ=255)
                                                                                                               cTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                        <  1:201038/140‑1 (MQ=255)
                                                                                                               cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                        >  1:181723/1‑140 (MQ=255)
                                                                                                               cTTCGGTTTCGGTATAATTTTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                        <  2:268190/140‑1 (MQ=255)
                                                                                                                 tCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                      >  2:404752/1‑140 (MQ=255)
                                                                                                                  cTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCt                     >  2:82204/1‑140 (MQ=255)
                                                                                                                         ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct              >  1:562421/1‑140 (MQ=255)
                                                                                                                         ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct              >  1:129268/1‑140 (MQ=255)
                                                                                                                           tgtAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct            <  1:82204/140‑1 (MQ=255)
                                                                                                                           tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct            <  1:607989/140‑1 (MQ=255)
                                                                                                                              aataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc         <  1:883984/140‑1 (MQ=255)
                                                                                                                                   ttttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    <  2:942546/138‑1 (MQ=255)
                                                                                                                                   ttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    <  2:562421/138‑1 (MQ=255)
                                                                                                                                   cGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    >  2:924024/1‑140 (MQ=255)
                                                                                                                                     tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  2:968771/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑C‑GCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  NC_000913/4111613‑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: ‑