Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,295,9760CT26.9% 31.8 / 11.5 26intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base T (5/2);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCG  >  NC_000913_3_pae_tpiA/4295846‑4296111
                                                                                                                                  |                                                                                                                                       
gCCGCATCCGACATCAACGCCAGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCCACGCCTg                                                                                                                                    >  1:213285/1‑136 (MQ=12)
     gtCCGACATCAACGCCTGATGCGACGGTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                                                                                               <  2:256665/135‑1 (MQ=17)
     aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGTCGCTTTTGACATCAACGCCTGATGCg                                                                                                                               <  1:142272/136‑1 (MQ=255)
      tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                                                                                              >  2:209964/1‑136 (MQ=17)
        cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAACGCAATAGCCTGATAAAAACAACTGTGGTGTGTCGGATAAGGCGATAACGCCGCATCCGACATAAATCCCTGATGCGACg                                                                                                                            >  1:52941/1‑136 (MQ=2)
                   cttGATGCGACGCTTGCTGCGCCTTATAAGGCCTTCGCCAGATAGCGCAAGAGCCTGATATAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                 <  2:128966/134‑1 (MQ=11)
                    cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                >  2:136857/1‑136 (MQ=21)
                      gatgcgatgcTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTTTGCCGCATCTGTCATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                              <  1:32365/136‑1 (MQ=11)
                      gATGCGACGCTGAACGCGTCGTATCAGGACTACCCCTGACGCACCAATAGCATAATTGAGGGGGATTTTGTTGGTGGTAGGAGGGGTGTTGGCCGGATGCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                              <  2:217030/136‑1 (MQ=2)
                       aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                             >  1:167766/1‑136 (MQ=21)
                            aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGTCGTTTGCTGCGTCTTATCAGGc                                                                                                        <  1:199041/136‑1 (MQ=255)
                             cGCTTAACGCGTCTTATCAGGCCTACGCCAAAAAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGACTGATGCGCCGCTTGCCGCGTCTTATCAGGcc                                                                                                       >  1:152867/1‑136 (MQ=11)
                                ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATATAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACACCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTAc                                                                                                    >  1:175013/1‑136 (MQ=12)
                                   aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctt                                                                                                 >  2:181526/1‑132 (MQ=9)
                                      cgTCATATCAGGCCCACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                              <  1:149668/136‑1 (MQ=255)
                                         cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGCACTACGCCAGAccg                                                                                           >  2:256501/1‑134 (MQ=255)
                                          ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                                                                                          <  1:256501/136‑1 (MQ=38)
                                          ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                                                                                          >  1:251685/1‑136 (MQ=38)
                                                     aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCATCGCAATAGTCTg                                                                               >  1:97103/1‑136 (MQ=255)
                                                                              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                                                                                                    <  2:250574/58‑1 (MQ=35)
                                                                                   gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                          >  2:56445/1‑111 (MQ=255)
                                                                                   gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                          <  1:56445/111‑1 (MQ=255)
                                                                                                tctgATCAGGCTTTTACGCCGCATCCGACATCAACGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGTCAGACAGCGTAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCa                                    <  1:264150/133‑1 (MQ=11)
                                                                                                 cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGGAATAGCCTGATTTAGCGTGATTTTCTAGGTCGGATAAGGCGTTTATGcaccat                                   >  2:125682/1‑131 (MQ=255)
                                                                                                                          gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAAACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGAATAAGGCGTTTCTGCCGCATCCGACATCAACGCCCGACGCGAccc          >  2:18831/1‑134 (MQ=11)
                                                                                                                           gcATCAACCCCTGATGCGACGCTTGCCGCGTCCTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGTCGCATCCGACATCAACGCCTGATGCGACGCt         <  1:254860/135‑1 (MQ=255)
                                                                                                                                  cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgcg  >  2:19994/1‑136 (MQ=11)
                                                                                                                                  |                                                                                                                                       
GCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCG  >  NC_000913_3_pae_tpiA/4295846‑4296111

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