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

TTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111618‑4111885
                                                                                                                                  |                                                                                                                                          
tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccg                                                                                                                                  >  1:71878/1‑140 (MQ=255)
   aTTAAACCCTTCAAGCGATTCAGATGCCCCCATTTTAATCGGTTCAACAGCGAGGGCCACCCCTTTTC‑CGCCGTTGCCCGCAACAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTTCCccgccg                                                                                                                               <  2:198247/140‑1 (MQ=255)
             acAAGCAATTAAGATAGCCCCAGTTTAATGGTTTAAACAGCGAAGGCCACCCCCTTTC‑CGCCGTTGCCCGCAAAAACATACCGCAGTCCCCCTTCTGTCGGGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGa                                                                                                                     <  2:211319/139‑1 (MQ=255)
              cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                    >  2:807584/1‑140 (MQ=255)
              cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                    >  2:731337/1‑140 (MQ=255)
              cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                    >  1:577524/1‑140 (MQ=255)
              cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                    >  1:163081/1‑140 (MQ=255)
                aGCGATTCAGATAGCGCCAGTTTAATGGGTTCAACAGCGAGGGTCACCCCTTTTC‑CGCCGTTGTCCGCGACAACAAACCGCGGTGCCCCTTCTGTCTGGGTGAAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTg                                                                                                                  <  2:71878/140‑1 (MQ=255)
                     ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                             >  1:771762/1‑140 (MQ=255)
                            ggcgcCAGCTTAAGCGGTACAACAGCGAAGGTCACCCCTTTTT‑CCCCGTTGTCCGCGAAAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATAAGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTg                                                                                                      <  2:771762/139‑1 (MQ=255)
                            aGCGCAACCTTAACGGGTTCAACAGCGAAGGCCACCCCCTTTT‑CGCCGTTGTCCGCAAAAACAAAACCCAGTGCCCTTTCTGTCTGGGTTTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTg                                                                                                      <  2:399421/140‑1 (MQ=255)
                                  aGCTTAATCGGTCCAACAGGGAAGGCCCCCCCTTTTT‑CGCCGTTGCCCGCGACAACATACCGCGGTGCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtg                                                                                                <  2:753586/140‑1 (MQ=255)
                                             ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     >  1:753414/1‑140 (MQ=255)
                                              tCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                                                                                    >  2:666723/1‑140 (MQ=255)
                                              tCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                                                                                    >  2:173789/1‑140 (MQ=255)
                                                       aaGGTCAGCCCTTTTC‑CGCCGTTGCCCGCGACAACATACCGCGGTGCCCTTTCTGTCTGGGTTAATTATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:18259/140‑1 (MQ=255)
                                                        aaggcccccccTTTTCCGCCGTTGCCCGCGAAAACATAACGCAGTCCCCCTTCTGTCTCGGTTTAATATCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:163081/133‑1 (MQ=255)
                                                            cAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtg                                                                      >  1:130790/1‑140 (MQ=255)
                                                                          ccGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGc                                                         >  2:251022/1‑140 (MQ=255)
                                                                           cGTTGCCCGCAAAAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                        <  2:753414/140‑1 (MQ=255)
                                                                             ttGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCTTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGc                                                      >  2:735636/1‑140 (MQ=255)
                                                                                tCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGcctgtct                                                   >  1:131474/1‑134 (MQ=255)
                                                                                 ccGCAAAAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg                                                  <  1:122359/140‑1 (MQ=255)
                                                                                 ccGCAAAAACAAAACGCATCCCCCTTTCTGTCTCGGTAAAATACCTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg                                                  <  1:540915/140‑1 (MQ=255)
                                                                                                    tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:115690/1‑140 (MQ=255)
                                                                                                     gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                              >  1:924833/1‑140 (MQ=255)
                                                                                                         cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                          >  1:7458/1‑140 (MQ=255)
                                                                                                                     tttaatatCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGTAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  1:807584/138‑1 (MQ=255)
                                                                                                                         ataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          >  2:753393/1‑140 (MQ=255)
                                                                                                                         ataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          >  2:855090/1‑140 (MQ=255)
                                                                                                                             ttttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  2:868606/138‑1 (MQ=255)
                                                                                                                              gtttgtttTTCCCCGCCTTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTgg     <  1:735636/140‑1 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:554226/1‑140 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:507009/1‑140 (MQ=255)
                                                                                                                                 tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc               >  2:230832/1‑127 (MQ=255)
                                                                                                                                 tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc               <  1:230832/127‑1 (MQ=255)
                                                                                                                                 tgtttTTCCCCGCCGTTAGGAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGGTTTCTCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:672868/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                          
TTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111618‑4111885

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