Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT47.8% 7.7 / 64.2 47intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (6/18);  new base T (7/15);  total (14/33)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGT  >  NC_000913/4295994‑4296127
                                                                  |                                                                   
cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                    <  2:553885/68‑1 (MQ=35)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                   >  2:40557/1‑68 (MQ=25)
    aTAGCCTGATTAAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                >  1:899652/1‑68 (MQ=255)
      agacTGATTTAGCGTGATTTTGTAGGTCGGAAAAGGCGTTTACGCCGCAACCGACATCAACGACTGAt                                                              >  2:779483/4‑68 (MQ=11)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                              <  1:809285/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                              <  1:809284/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                              >  2:789329/1‑68 (MQ=35)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                              <  1:420287/68‑1 (MQ=35)
             tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGc                                                       <  2:202820/68‑1 (MQ=35)
                  cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                  >  1:929268/1‑68 (MQ=255)
                     gATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                               <  2:348058/68‑1 (MQ=25)
                     gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                               <  1:654678/68‑1 (MQ=255)
                     gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                               <  1:383370/68‑1 (MQ=255)
                     gATTTTGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACAGCAATGCCTGATGCGACGCTTGCcgcg                                               <  1:682067/68‑1 (MQ=14)
                     gATTTTGAAGTTCGGATAAGGCGTTTACGCCGCATCCTACATCAACGCCTGATGCGACGCTTGCcgcg                                               <  2:183964/68‑1 (MQ=14)
                      tttttGTAGGTCGGATAATGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGt                                              <  2:817580/67‑1 (MQ=14)
                      aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                              <  2:818397/68‑1 (MQ=255)
                       ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAAGCTTGACGCGTc                                             >  2:51409/1‑68 (MQ=255)
                        tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                            <  2:910856/68‑1 (MQ=255)
                        tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTAGACGCGTCt                                            >  2:737717/1‑68 (MQ=11)
                          tGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                          <  1:544527/68‑1 (MQ=255)
                           gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                         >  1:883871/1‑68 (MQ=255)
                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                       >  2:800482/1‑68 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                     <  2:883875/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                     <  2:591247/68‑1 (MQ=255)
                               gTCGCATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATTCGACGCTTGCCGCGTCTTATCAgg                                     <  2:576628/68‑1 (MQ=14)
                                gcGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                    <  2:307821/67‑1 (MQ=255)
                                  ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                  >  2:746359/1‑68 (MQ=255)
                                    ttAAGGCGTTTACGCCGCATCCGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                <  1:13420/67‑1 (MQ=255)
                                    ctAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                <  2:705016/67‑2 (MQ=255)
                                    ctAAGGCGTTTACGCCGCATCCGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                <  1:876428/67‑1 (MQ=255)
                                    aTAAGTCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                <  2:156346/68‑2 (MQ=14)
                                    aTAAGGCGTTTACGCCGCATCCGACATCATTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                <  2:790650/68‑2 (MQ=14)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                <  1:579101/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                <  1:19158/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                <  1:19157/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                <  2:232613/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                <  2:89493/68‑1 (MQ=255)
                                         gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCaa                           <  1:856770/68‑2 (MQ=255)
                                          cGTATACGCCGCATCAGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGACAGa                          >  2:883515/1‑68 (MQ=14)
                                              tACGCCGCATCCGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                      <  2:342434/68‑1 (MQ=255)
                                                cgccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcgc                    <  1:472651/68‑1 (MQ=255)
                                                  ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGACAGACAGCGCaa                  >  1:419626/1‑68 (MQ=255)
                                                      aTCCGAGAGCAACGCCTGATGCGACGCTTTCCGCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGc              <  2:399923/68‑1 (MQ=11)
                                                      aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc              <  1:574912/68‑1 (MQ=40)
                                                            cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt        >  2:249270/1‑68 (MQ=40)
                                                            cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt        >  1:797857/1‑68 (MQ=40)
                                                                  cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGt  >  2:856510/1‑68 (MQ=31)
                                                                  |                                                                   
CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGT  >  NC_000913/4295994‑4296127

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