Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA4,296,0900CT22.5% 54.1 / 14.5 41intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (20/11);  new base T (5/4);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTA  >  NC_000913_3_vch_tpiA/4295945‑4296230
                                                                                                                                                 |                                                                                                                                            
ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                           >  2:311660/1‑149 (MQ=17)
ggAGAAGGCGTGTAGGCCGCAGGCGACATGAACGCCTGGTGCGACGCTTAACGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGcc                                                                                                                                           <  2:375990/149‑1 (MQ=11)
             aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCTTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCtt                                                                                                                              <  1:223018/149‑1 (MQ=255)
               gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGCCGCTTGc                                                                                                                            >  1:363675/1‑149 (MQ=18)
                    aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                       <  1:311660/149‑1 (MQ=21)
                     tCCGACATCAACGCCTGATGTGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                      <  1:55067/149‑1 (MQ=255)
                     tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                      >  1:489233/1‑149 (MQ=37)
                          cATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                 >  1:410470/1‑149 (MQ=25)
                             cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                              >  2:373113/1‑149 (MQ=38)
                             cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                              <  1:520213/149‑1 (MQ=39)
                               aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                            >  2:136073/1‑149 (MQ=21)
                                   cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                        <  1:52497/149‑3 (MQ=16)
                                   cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCATatc                                                                                                        >  1:466976/1‑147 (MQ=12)
                                   cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                                                                                                        >  1:231397/1‑149 (MQ=37)
                                    tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCATGATTTTGTAGGTCAGAAAAGACGTTTACGCCGCATCCGACATCAACGCGTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTAcac                                                                                                       >  1:140715/1‑147 (MQ=11)
                                        gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                   <  1:323900/149‑1 (MQ=38)
                                               ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGTCGAATAAGGAGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCCGCGCa                                                                                            >  1:267063/1‑149 (MQ=11)
                                               ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGACAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAAACAGTGCa                                                                                            >  2:352811/1‑149 (MQ=255)
                                                         ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTCAGCCAGACAGCGCAATAGCCTAAt                                                                                  >  1:117178/1‑149 (MQ=11)
                                                         ttATCAGGCCTACGCCAGACAGCGCAATACCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAGGCGACGCTAGCCGCGTCTTATCAGGCCTCCGCACGCCAGCCAAACAGCCtaca                                                                                  >  2:221851/1‑146 (MQ=11)
                                                              aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTAATGCGACGCTTGCCGCGTCTTATCAGGCCTACTCCAGACAGCGCACTAGCGTGATTTatc                                                                             >  1:111743/1‑147 (MQ=255)
                                                                         aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                  <  2:489233/149‑1 (MQ=38)
                                                                          gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCAGCCTCCCACATCAACGCCTGATGCGACGCTTGCCGCGTATTATCAGGCCTACGCCATCCAGTCCAATAACGTGATTTAGCGGGATTTTGTAg                                                                 >  1:47146/1‑149 (MQ=11)
                                                                                  aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATa                                                         >  2:292740/1‑149 (MQ=38)
                                                                                      gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                                     <  1:192745/149‑1 (MQ=37)
                                                                                        cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGGCTTATCAGGCCTACGCCAGAAAGCGCAATAGCCTGATTTCACGCAATTCTGTCGATCGGATAAGCCgc                                                   >  2:73043/1‑148 (MQ=11)
                                                                                        cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCCCATCCGTCTTAAGCGGATCAGTCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGt                                                   <  1:542724/149‑1 (MQ=11)
                                                                                         tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAACCTGATTTAGCGTAACTTTACAGGTCGAACAAGGCGtg                                                  >  1:266628/1‑148 (MQ=11)
                                                                                         tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATACGACATCAACGCCTGATGCTACGCTTGCCGCGTCTTAGAAGGCCTACTCCAGCCAGCGCAATAGCCTGATTTAGCGCGATTTTGTAGGTCGTGTCAGGCGtt                                                  >  1:500191/1‑149 (MQ=11)
                                                                                                          gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCACGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCCGATAAGGCGCTTATGCCGCATCCGACac                                 >  2:298150/1‑148 (MQ=255)
                                                                                                              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGc                                          >  2:343472/1‑136 (MQ=40)
                                                                                                              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGc                                          <  1:343472/136‑1 (MQ=40)
                                                                                                              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCAGATAGGGCGTTTATGCCGCAGCCGACATCAAc                             >  2:3689/1‑149 (MQ=255)
                                                                                                                       ggCGTTTACGCCGCATCCGCTATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTTGTCTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGc                    <  1:106072/149‑1 (MQ=11)
                                                                                                                       ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                                                                                              <  2:102001/43‑1 (MQ=255)
                                                                                                                       ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCtt                                                                                                                              >  1:102001/1‑43 (MQ=40)
                                                                                                                          gTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATGCCAGACAGTGCAATAGTCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAc                 <  1:78894/149‑1 (MQ=255)
                                                                                                                                 ccgcATCCGACATCAACGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCGTACGCCCGACAGCGCACTAGCCTGATTTAGCGTGAGTTTGTAGGCCGGATAAGGCGTTTATGCCGCATCCGAGATCAACGCCTGATGCGACGTGTAAc          >  1:466930/1‑149 (MQ=11)
                                                                                                                                  cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCAGACAGACAGCGCAATAACCTAATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAGCCGACATCAACGCCTGATGCGACGCTTAAcg         >  1:47178/1‑149 (MQ=11)
                                                                                                                                   gcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc        <  1:373113/149‑1 (MQ=35)
                                                                                                                                         gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTa  >  2:102880/1‑149 (MQ=35)
                                                                                                                                                 |                                                                                                                                            
GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTA  >  NC_000913_3_vch_tpiA/4295945‑4296230

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