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

ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGG  >  NC_000913/4111617‑4111886
                                                                                                                                  |                                                                                                                                           
ttttATTCAACCCTTCAAGCTATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCACCCCTTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCcccc                                                                                                                                    <  2:22967/139‑1 (MQ=255)
ttttATTAAAACCTCCAAGCTATTAAGATGGCGCCATCTTAATGGGTTAAACAGCGACGGTCCCCCCTTTTTCGCCGTTGCCCGCAACAACATACCGCAGTGCACCTTCTGTCTGGTTATAAAAACGTTTTTTTTTcccc                                                                                                                                    <  2:512225/139‑1 (MQ=255)
                aaGCGTTTAAGATAGCCCCAGCTTAACCGTTTCAACAGGGAAGCTCCCCCCTTTTTCGCCGTTGCCCGCGACAACATAACGCGGTCCCCCTTCTGTCTCGGTGTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACt                                                                                                                    <  2:42436/140‑1 (MQ=255)
                aaGCAATACAGGTACCGCGAGCTTAATTGTTTCAACAGCGAGGGTCACCCCCTTTCCGCCGTGCTCCGCCACACCATAGCGCCGTCTTCCTCCTGTCTGGGTGTAATAGTGTTTGTTTTTCCCCGCCGTTAGTAGCGACt                                                                                                                    <  2:1033754/140‑1 (MQ=255)
                      ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                              >  1:21691/1‑140 (MQ=255)
                      ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                              >  1:618444/1‑140 (MQ=255)
                             aGCGCTAGCTTAAGCGTTTCAACAGGAAAGCCCACCCTTTTTTCGCCTTTTCCCGCGAAAACATACGGCGGTCCCCCTTCTGGCTCGGTATAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTg                                                                                                       <  2:813038/140‑1 (MQ=255)
                                gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                    >  2:75522/1‑140 (MQ=255)
                                            ggtaaaaCAGCGAAGCCCCCCCCTTTTTCGCCTTTGCCCGCAAAAAAATACCGCTGTCCCCCTTCTGTCTGGGTATAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                        <  2:885650/135‑1 (MQ=255)
                                              ttCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                      >  1:409701/1‑140 (MQ=255)
                                              tccAACAGGGAGGGCCACCCCTTTTCCGCCGTTGTCCGCAAAAAAAAAACGCAGTCCCCTTTCTTTCTGGTTGTATTATTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                      <  2:1199313/138‑1 (MQ=255)
                                              tccAACAGCGAGGGCCAGCCCCTTTTCGCCGTTGCCCGCGACAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATATTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                      <  2:1013985/138‑1 (MQ=255)
                                                aaaCAGCGAGGGTCAGCCCTTTTCCGCCGTTGCCCGCAACAACATAACGCAGTGCCCCTTCTGTCTGGTTTAAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                    <  2:727803/139‑1 (MQ=255)
                                                  aCAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                                  >  2:328856/1‑140 (MQ=255)
                                                        aaGGTCCCCCCTTTTTCGCCGTTGCCCGCAACAAAATAACACAGTCCCCCTTCTTTCTCGGTGTAATAATTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            <  2:982016/140‑1 (MQ=255)
                                                        aaGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggggg                                                                            >  2:839494/1‑139 (MQ=255)
                                                        aaGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            >  1:651020/1‑140 (MQ=255)
                                                        aaGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            >  1:484057/1‑140 (MQ=255)
                                                        aaGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTG‑TTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTg                                                                           >  1:470090/1‑140 (MQ=255)
                                                        aaGGTCACCCCTTTTCCGCCGTTGCCCGCGAAAACATAACGCGTTCCCCTTTCTGTCTGGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            <  2:1170550/140‑1 (MQ=255)
                                                          ggTCCGCCCTTTTTCGCCGTTGCCCGCAACAACATAACGCAGTCCCTTTTCTGTCTTGGTGTAATACTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa                                                                          <  2:473650/140‑1 (MQ=255)
                                                                ccccTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc                                                                    >  1:142383/1‑140 (MQ=255)
                                                                            gTTGCCCGCGACAACAAACCGCGGTCCCCTTTCTGTCTGGGTTTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                        <  2:1040747/140‑1 (MQ=255)
                                                                                                  aGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgc                                  >  1:237582/1‑140 (MQ=255)
                                                                                                    tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTATTGGGGGGAGTGTGCGGATCACGGCTGCGATTTATTCCTGAgctttcgctt                                >  2:222527/1‑140 (MQ=255)
                                                                                                    tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                >  1:129080/1‑140 (MQ=255)
                                                                                                              gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                      >  1:480489/1‑140 (MQ=255)
                                                                                                              gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                      >  2:38416/1‑140 (MQ=255)
                                                                                                                   ggTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                 <  1:1014098/140‑1 (MQ=255)
                                                                                                                    gtataatatCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc                <  2:484057/137‑1 (MQ=255)
                                                                                                                     tgtAATATCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct               <  1:747235/140‑1 (MQ=255)
                                                                                                                       taataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg             >  1:126742/1‑140 (MQ=255)
                                                                                                                        atgatcGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATGAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc            <  2:977509/135‑1 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  2:273600/1‑140 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  2:926294/1‑140 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCCGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  2:1210781/1‑140 (MQ=255)
                                                                                                                                 tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg   >  1:309749/1‑140 (MQ=255)
                                                                                                                                 tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg   >  2:272355/1‑140 (MQ=255)
                                                                                                                                  gtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg  <  1:926294/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                           
ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGG  >  NC_000913/4111617‑4111886

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