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

TTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111612‑4111885
                                                                                                                                         |                                                                                                                                          
ttGTCTTTTATTCAACCCTTCAAGCGATTCAGATAGCCCCAGCTTAATGGGTCCAACAGCGAAGGCCACCCCTTTTC‑CGCCGTTGCCCGCA‑ACAACATACCGCAGTGCCCCTTCTGTCTCGGTGTAATACCTtttttttt                                                                                                                                        <  2:1071070/140‑1 (MQ=255)
ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                                                                                        >  1:1267661/1‑140 (MQ=255)
ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                                                                                        >  1:186065/1‑140 (MQ=255)
             aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTCTTCCCCGCCGTtag                                                                                                                           >  1:801083/1‑140 (MQ=255)
               aCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagta                                                                                                                         >  1:532331/1‑140 (MQ=255)
                    cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                    >  2:260475/1‑140 (MQ=255)
                    cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCATTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                    >  1:67/1‑140 (MQ=255)
                          tttCAGATAGCCCCAGTTTAATCGGTTCAACAGCGAAGGTCA‑CCCCTTTT‑GCGCGTTGGCCCCGAAAAACATACCTAAGTGCCCCTTCTTCCTCGGTATAAGAGGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGt                                                                                                              <  2:820227/139‑1 (MQ=255)
                           ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                             >  1:242659/1‑140 (MQ=255)
                                aTAGCGCCACTTTAATCGTTAAAAAGGCGAGGGCCACCCCTTTTC‑CGCGGTTGTCCGCA‑ACAAAATACCGCAGTCCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                                                                                        <  2:928805/140‑1 (MQ=255)
                                     gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   >  2:30469/1‑140 (MQ=255)
                                     gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   >  1:820227/1‑140 (MQ=255)
                                                 gtttCAACAGCGAAGGCCCCCCCTTTTT‑CGCCGTTGCCCGCG‑ACAACATAACGCGGTCCCCCTTCTGTCTGGGTGTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:293626/138‑1 (MQ=255)
                                                 gttacAACAGCGAGGGCCCGCCCTTTTC‑CGCGGTTGCCCGAA‑AAAACATAGCGCAGTCCCCCTTGTTTCTGGGTTTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:1267661/136‑1 (MQ=255)
                                                 ggTTCAACAGCGAGGGCCCCCCCTTTTT‑CGCCGTTGCCCGCA‑ACAACATACCGCAGTGCCCCTTCTGTCTCGGTTTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:522760/140‑1 (MQ=255)
                                                   tccAACGGCGAGGGCCACCCCTTTTT‑CCCCTTTGTCCGCT‑ACAACAAAACGCAGTGCCCCTTCTGTCTCGGTGTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     <  2:128164/138‑1 (MQ=255)
                                                             agggTCACCCCCTTTG‑CGCTTTTGTCCGCG‑AAAAAATAACGGAGTGCCCCTCCTCTCTGGTTTAAATAATGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:1207990/138‑1 (MQ=255)
                                                             aaGGTCCCCCCTTTTT‑CGCTTTTGTCCGCG‑ACAACATACCGCAGTCCACCTTCTGCCGGGGTTTAATATCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:1021412/140‑1 (MQ=255)
                                                             aaGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           >  1:788431/1‑140 (MQ=255)
                                                             aaGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           >  1:44168/1‑140 (MQ=255)
                                                             aaGGCCCGCCCTTTTT‑CGCCTTTGTCCGCG‑ACAACATAACGCGGTGCCCTTTCTGTCTGGGTTTATTACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:67/140‑1 (MQ=255)
                                                                gTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAg                                                                        >  2:585098/1‑140 (MQ=255)
                                                                   ccccccTTTTCCGCGGTTGCCCGCG‑ACAACATACCGCGGTGCCCCTTCTGTCTGGGTTTAATACCGTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtg                                                                      <  2:532331/138‑1 (MQ=255)
                                                                    cccctTTTT‑CCCCTTTGCCCGCA‑ACAACATACCGCTTTGCCCCTCTTGTCTGGGTTTAATAGTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGtgtg                                                                    <  2:788431/139‑1 (MQ=255)
                                                                               ccctttGCCCGCA‑AAAACATACCGCAGTCCCCCTTCTGTCTTGTTTAAATAATTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                          <  2:242659/136‑1 (MQ=255)
                                                                                 cGTTGTCCGCG‑ACAACATAACGCAGTCCCCTTTCTGTCTGGGTGTAATAGTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                        <  2:1038768/140‑1 (MQ=255)
                                                                                 cGTTGCCCGAA‑AAAAAATACCGCAGTCCCCTTTTTTTCTGGTTGTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                        <  2:919023/140‑1 (MQ=255)
                                                                                      cccGCA‑AAAACATACCGCAGTGCCCTTTCTGTCTGGGTGTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCAt                                                   <  2:1242601/139‑1 (MQ=255)
                                                                                           g‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTAtt                                              >  2:1025509/1‑140 (MQ=255)
                                                                                                      cGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgcttt                                    >  1:504741/1‑140 (MQ=255)
                                                                                                        cAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcg                                  >  2:987760/1‑140 (MQ=255)
                                                                                                         aGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgc                                 >  1:818445/1‑140 (MQ=255)
                                                                                                         aGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgc                                 >  1:818627/1‑140 (MQ=255)
                                                                                                           tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:185045/1‑140 (MQ=255)
                                                                                                           tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:188350/1‑140 (MQ=255)
                                                                                                                  tCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                        >  2:333427/1‑140 (MQ=255)
                                                                                                                     gTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                     <  1:30469/140‑1 (MQ=255)
                                                                                                                     gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                     >  2:22538/1‑140 (MQ=255)
                                                                                                                               aataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc           <  1:585098/140‑1 (MQ=255)
                                                                                                                                atatCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          <  1:1079273/140‑1 (MQ=255)
                                                                                                                                ataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          >  2:157571/1‑140 (MQ=255)
                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:999562/1‑140 (MQ=255)
                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:1004688/1‑140 (MQ=255)
                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:211457/1‑140 (MQ=255)
                                                                                                                                        tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          >  2:43809/1‑132 (MQ=255)
                                                                                                                                        tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          <  1:43809/132‑1 (MQ=255)
                                                                                                                                        tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:1090544/1‑140 (MQ=255)
                                                                                                                                        tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:109007/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                          
TTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111612‑4111885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑