Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT29.4% 50.3 / 24.6 35intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (9/15);  new base T (4/6);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA  >  NC_000913/4295995‑4296123
                                                                 |                                                               
gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                               >  2:355694/1‑68 (MQ=25)
 cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGcc                                                              <  1:447099/68‑1 (MQ=35)
     aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAATGCCTGAt                                                          >  1:794974/1‑68 (MQ=255)
          gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                     <  2:549081/68‑1 (MQ=25)
             ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCt                                                  >  1:298486/1‑68 (MQ=25)
                 cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                              >  1:65052/1‑68 (MQ=40)
                  gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcg                                             >  2:257175/1‑68 (MQ=255)
                    gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACAGCAACGCCTGATGCGACGCTTGCcgcg                                           <  2:700627/68‑1 (MQ=255)
                       tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt                                        <  2:65052/68‑1 (MQ=255)
                        ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                       <  2:314647/68‑1 (MQ=255)
                         tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                      <  2:378935/68‑1 (MQ=255)
                              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 >  2:174531/1‑68 (MQ=255)
                              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  2:51810/68‑1 (MQ=255)
                               tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                >  2:578276/1‑68 (MQ=255)
                                   ttAAGGCGTTTACGCCGCTTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:732386/67‑1 (MQ=255)
                                   aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                            <  2:600641/68‑2 (MQ=255)
                                   aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                            <  1:46193/68‑2 (MQ=255)
                                   aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                            <  2:428690/68‑2 (MQ=255)
                                   aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  2:561521/68‑1 (MQ=255)
                                   aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  2:140306/68‑1 (MQ=255)
                                   aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  2:278762/68‑1 (MQ=255)
                                   aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:754378/68‑1 (MQ=255)
                                   aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  2:487511/68‑1 (MQ=255)
                                   aTAAGGCGTTTACGCCGCATCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                            <  1:719598/68‑2 (MQ=14)
                                        gCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  2:403470/58‑1 (MQ=255)
                                        gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 >  1:403470/1‑58 (MQ=255)
                                        gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAg                       <  1:734002/68‑1 (MQ=255)
                                         cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                      >  2:31546/1‑68 (MQ=255)
                                         cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                      >  2:183772/1‑68 (MQ=255)
                                             tACGCCGCATCCGACATCACCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                  >  1:521598/1‑68 (MQ=255)
                                                gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCa               <  2:706577/68‑1 (MQ=40)
                                                 ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                       >  2:183874/1‑43 (MQ=38)
                                                     ctCCGACATCAACTTCTGATGCGACGCTTTCCTCTGCTTATCAGTCCTACGCCAGACAGCGCAATAGc          <  2:790959/67‑1 (MQ=11)
                                                           cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGAGCAATAGCCAGAtt    >  2:2890/1‑68 (MQ=255)
                                                             tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa  >  2:163743/1‑68 (MQ=39)
                                                                 |                                                               
GCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA  >  NC_000913/4295995‑4296123

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