Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913508,7450TG20.0% 61.7 / 10.7 30intergenic (+186/+130)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base T (12/12);  new base G (4/2);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGG  >  NC_000913/508604‑508886
                                                                                                                                             |                                                                                                                                             
gccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCGACAACCTTGTGCCTGAACCGTAGGCCGGATAAGGCGGTTACGCCGCCTCCCGCAGCCGAGCACCATTGCCTGAAGCGACGCGTAcgcga                                                                                                                                        >  8:93755/1‑148 (MQ=17)
                 aGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACaaa                                                                                                                       >  3:16220/1‑149 (MQ=255)
                    aCCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCt                                                                                                                    <  4:16220/149‑1 (MQ=255)
                     ccATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCtg                                                                                                                   >  6:18777/1‑149 (MQ=255)
                         tGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGcc                                                                                                               >  3:108543/1‑149 (MQ=255)
                          gCCTGATAAGACGGGTAAGCGTCGCATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCt                                                                                                              >  6:19566/1‑149 (MQ=18)
                             tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAAc                                                                                                                      >  7:127595/1‑138 (MQ=255)
                             tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAAc                                                                                                                      <  8:127595/138‑1 (MQ=255)
                                       cTTACCCGGCTTTTCAGGCCTAGATTCTTGTGCCTGAACCGTAGGCCGGATAAGGCGATTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCg                                                                                                 <  7:93755/149‑1 (MQ=17)
                                          aCCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATAAGACGCGTAAGCGTcgcatc                                                                                                                                  <  5:173026/113‑6 (MQ=18)
                                          aCCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATAAGACGCGTAAGCGTcgcatc                                                                                                                                  >  6:173026/1‑108 (MQ=18)
                                            ccGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCtg                                                                                                                   >  8:95268/1‑126 (MQ=255)
                                            ccGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCtg                                                                                                                   <  7:95268/126‑1 (MQ=255)
                                               tcgcaTCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGc                                                                                         <  5:19566/145‑1 (MQ=255)
                                                          cTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGTCGGATAAGGCGTTTa                                                                                    <  1:178817/143‑1 (MQ=255)
                                                          cTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGTCGGATAAGGCGTTTa                                                                                    >  2:178817/1‑143 (MQ=255)
                                                                aCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                        <  4:108543/149‑1 (MQ=255)
                                                                     tgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGt                                                                    >  5:100433/1‑148 (MQ=255)
                                                                         tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGc                                                               >  8:18028/1‑149 (MQ=255)
                                                                                  ggCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGc                                                       >  4:107629/1‑148 (MQ=255)
                                                                                                      ccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAcc                                          >  8:27860/1‑141 (MQ=255)
                                                                                                      ccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAcc                                          <  7:27860/141‑1 (MQ=255)
                                                                                                        gcatccggcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAagcagc                                <  6:100433/149‑1 (MQ=255)
                                                                                                          atccggcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAgcgc                              <  7:18028/149‑1 (MQ=255)
                                                                                                                      aGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAatgat                  <  5:18777/149‑1 (MQ=255)
                                                                                                                          ccATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTAc              >  5:39010/1‑149 (MQ=255)
                                                                                                                           gtaggCCTGATAAGACGCGTAAGCGTCGCATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACt             >  5:151369/5‑149 (MQ=25)
                                                                                                                               gCCTGATAAGACGCGTAAGCGTCGCATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACt             <  7:28447/145‑1 (MQ=25)
                                                                                                                               gCCTGATAAGACGCGTAAGCGTCGCATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACt             >  8:28447/1‑145 (MQ=25)
                                                                                                                                       cGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCgg  <  3:107629/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
GCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGG  >  NC_000913/508604‑508886

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑