Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT30.0% 66.5 / 15.2 30N156K (AAC→AAA) yiiQDUF1454 domain‑containing protein YiiQ
Reads supporting (aligned to +/- strand):  ref base G (19/2);  new base T (0/9);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 3.84e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111610‑4111885
                                                                                                                                           |                                                                                                                                          
ttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACgtttgtt                                                                                                                                          <  2:1122273/140‑1 (MQ=255)
       ttttATTAAACCCTAAAAGCGATTCAGATAGCGCCAGTTTAATCGTTCAAACAGCGAGGGCCACCCCTTTTT‑CGCCGTT‑GCCCGCGAAAACATACCGCAGTGCCCTTTCTGTCTCGGTGTAATAACTTTTTTTTTTcccc                                                                                                                                   <  2:746943/139‑1 (MQ=255)
                      cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                    >  1:101215/1‑140 (MQ=255)
                           gtttCAGATAGCGCCAGTTTACTTGTTTCAACAGCGAGGGTCCCCCTTTTTT‑CGCGTTT‑GTCCGCGACAAAATAAGGCAGTGCCCTTTTTTTTTGGTTTTAATACTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAg                                                                                                               <  2:817851/138‑1 (MQ=255)
                           gATTAAGATGGCGCCAGTTTAATCGGTCAAACAGCGAGGGTCACCCCTTTTC‑CGCCGTT‑GCCCGCGAAAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAg                                                                                                               <  2:122523/140‑1 (MQ=255)
                             ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                             >  1:738377/1‑140 (MQ=255)
                                  aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                                                                                        >  1:110415/1‑140 (MQ=255)
                                       gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGcc                                                                                                   >  1:817851/1‑139 (MQ=255)
                                       gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   >  1:1072393/1‑140 (MQ=255)
                                       cccAGTTTAATCGGTTAAACAGCGAAGGCCACCCCTTTTT‑CGCCGTT‑GCCCGCAACAACATAACGCGGTGCCCCTTCTGTCTGGGTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   <  2:1072393/139‑1 (MQ=255)
                                         cAGTTTAATGGTTACAACAGCAAGGGCCCCCCCTTTTT‑CGCTGTT‑GCCCGCGACAACATACCGCAGTGCCCCTTCTGTCTGGGTTTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                                                                                 <  2:801283/140‑1 (MQ=255)
                                                 tCGGTCCAACAGCGAGGGCCCCCCCTTTTT‑CGCCGTT‑GCCCGCAACAACAAAACGCGGTGCCCTTTCTGTCTGGGTGTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                         <  2:712830/140‑1 (MQ=255)
                                                   ggTTAAACAGCGAGGGTCCGCCCTTTTT‑CGCTGTT‑GCCCGCAACAACATAACGCGGTCCCCTTTCTGTCTGGGTGTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:738377/140‑1 (MQ=255)
                                                   aaaaaaaCAGCGAAGGCCCCCCCTTTTT‑CGCGGTG‑GCCCGCAAAAACAAACCGCGTTCCCCCTTTTTTCTGGGTTTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:178966/135‑1 (MQ=255)
                                                             cGAAGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                                                                             >  1:750527/1‑140 (MQ=255)
                                                               aaGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           >  1:952138/1‑140 (MQ=255)
                                                               aaGGCCACCCCTTTTC‑CGCCGTT‑GCCCGCAACAACATAACGCAGTGCCCCTTCTGTCTGGTTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:520779/140‑1 (MQ=255)
                                                                acggccgCCCCTTTTCCGCCGTG‑CCCCGCCACAACATAACCCCGTCCCCCTTCGCTCTGGGTATATTATTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:750527/134‑1 (MQ=255)
                                                                    cAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggggctgtct                                                                      >  1:178966/1‑134 (MQ=255)
                                                                       ccccTTTT‑CCCCGTTGGCCCGCAAAAAAAAAACGCTGTCCCCCTTTTTTCTGGGTTAAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGtgtg                                                                    <  2:101215/140‑1 (MQ=255)
                                                                                   gttt‑GGCCGCGACAACATACCGCAGTGCCCCTTCGGCTTGGTTTTAATACCTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                        <  1:661258/138‑1 (MQ=255)
                                                                                    gTT‑GTCCGCGAAAACATAACGCAGTGCCCTTTCTGTCTGGGTTTAATATCTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                       <  1:583486/140‑1 (MQ=255)
                                                                                       tgcccGCGACAACATACCGCGGTCCCCCTTCTGTCTGGGTATAATATTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTTGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                     <  2:870152/137‑1 (MQ=255)
                                                                                           cgcgACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGt                                                 >  1:15528/1‑140 (MQ=255)
                                                                                                 aacaTACCGCAGTCCCCTTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCt                                           <  1:939977/140‑1 (MQ=255)
                                                                                                        cGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgcttt                                    >  1:533620/1‑140 (MQ=255)
                                                                                                             tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:537476/1‑140 (MQ=255)
                                                                                                             tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  2:292275/1‑140 (MQ=255)
                                                                                                              gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                              >  2:919783/1‑140 (MQ=255)
                                                                                                                cccTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGctttcgctttcg                            <  1:124826/139‑1 (MQ=255)
                                                                                                                  tttCTGTCTGGTTTTATTATTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                          <  2:746458/139‑1 (MQ=255)
                                                                                                                       gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                     >  1:703289/1‑140 (MQ=255)
                                                                                                                       gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                     >  1:746458/1‑140 (MQ=255)
                                                                                                                            ggTGTAATGACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                >  1:917942/1‑140 (MQ=255)
                                                                                                                            ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                >  1:390817/1‑140 (MQ=255)
                                                                                                                             gtgtAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc               <  2:952138/140‑1 (MQ=255)
                                                                                                                              tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              >  1:908631/1‑140 (MQ=255)
                                                                                                                              tataataaCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  1:1080262/138‑1 (MQ=255)
                                                                                                                                  atacctttttttttTCGCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          <  2:533620/134‑1 (MQ=255)
                                                                                                                                      ttttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  1:671909/138‑1 (MQ=255)
                                                                                                                                      cGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      >  2:961571/1‑140 (MQ=255)
                                                                                                                                        ttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          <  1:1096584/134‑1 (MQ=255)
                                                                                                                                        tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          >  2:1096584/1‑134 (MQ=255)
                                                                                                                                        tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:141302/1‑140 (MQ=255)
                                                                                                                                        tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:827062/1‑140 (MQ=255)
                                                                                                                                        tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:866277/1‑140 (MQ=255)
                                                                                                                                          tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:1047981/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                          
TTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTT‑GTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111610‑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: ‑