Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA601,4130TA27.3% 37.7 / 12.4 22I900I (ATT→ATAcusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (10/6);  new base A (3/3);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTATCGCCCTCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGC  >  NC_000913_3_pae_tpiA/601288‑601546
                                                                                                                             |                                                                                                                                     
tcCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCtg                                                                                                               >  4:62581/1‑150 (MQ=255)
          tCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCa                                                                                                     <  3:62581/150‑1 (MQ=255)
                      aGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATAATcagcaac                                                                                                                              >  2:1098/1‑111 (MQ=255)
                      aGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATAATcagcaac                                                                                                                              <  1:1098/113‑3 (MQ=255)
                            ccAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGa                                                                                   >  4:77167/1‑150 (MQ=255)
                             cAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGAt                                                                                  >  1:48427/1‑150 (MQ=255)
                                  aTAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATAATcagca                                                                                                                                <  1:77966/99‑1 (MQ=255)
                                  aTAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATAATcagca                                                                                                                                >  2:77966/1‑99 (MQ=255)
                                          aaaCTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCtt                                                                     >  4:90859/1‑150 (MQ=255)
                                              tCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTa                                                                                                                           >  2:98757/1‑92 (MQ=255)
                                              tCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTa                                                                                                                           <  1:98757/92‑1 (MQ=255)
                                              tCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCt                                                                                                  >  3:36245/1‑117 (MQ=255)
                                              tCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCt                                                                                                  <  4:36245/117‑1 (MQ=255)
                                               cATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGt                                                                <  2:48427/150‑1 (MQ=255)
                                                     gCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGCTGCTGATAATcagcaac                                                                                                                              <  2:126414/82‑3 (MQ=255)
                                                     gCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGCTGCTGATAATcagcaac                                                                                                                              >  1:126414/1‑80 (MQ=255)
                                                          tGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCa                                                     >  2:113396/1‑150 (MQ=255)
                                                                                   gttgtATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTATCGCCCTCGCCGGGGt                            <  3:90859/150‑1 (MQ=255)
                                                                                    ttgtATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTATCGCCCTCGCCGGGGTc                           >  1:78353/1‑150 (MQ=255)
                                                                                         tCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTATCGCCCTCGCCGGGGTcgccgc                      >  4:113857/1‑150 (MQ=255)
                                                                                                     tCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTATCGCCCTCGCCGGGGTCGCCGCCGAATTTGGCgt          <  3:113857/150‑1 (MQ=255)
                                                                                                             gCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTATCGCCCTCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGc  >  3:75363/1‑150 (MQ=255)
                                                                                                                             |                                                                                                                                     
TCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTATCGCCCTCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGC  >  NC_000913_3_pae_tpiA/601288‑601546

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