Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT44.0% 35.1 / 99.3 75intergenic (+266/+376)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base T (12/21);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC  >  NC_000913/4295994‑4296125
                                                                  |                                                                 
cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                  <  1:668997/68‑1 (MQ=35)
cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                  <  1:340800/68‑1 (MQ=35)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                 >  1:668942/1‑68 (MQ=25)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                 <  1:1536261/68‑1 (MQ=35)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                 >  1:1683074/1‑68 (MQ=35)
    aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                              >  1:2303618/1‑68 (MQ=25)
    aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                              >  1:2191049/1‑68 (MQ=25)
    aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                              >  1:2601803/1‑68 (MQ=35)
    aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                              >  1:1884705/1‑68 (MQ=35)
     tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                             <  1:726366/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATGAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                            <  1:2286598/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            >  1:3029701/1‑68 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            >  1:2676891/1‑68 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            >  1:2499649/1‑68 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            <  1:3399017/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            >  1:1786457/1‑68 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            >  1:645599/1‑68 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            >  1:1342366/1‑68 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                            >  1:1922456/1‑68 (MQ=35)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                            >  1:3330386/1‑68 (MQ=35)
      aGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            >  1:2128715/1‑68 (MQ=255)
        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                          >  1:2487770/1‑68 (MQ=35)
         cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                         >  1:1215689/1‑68 (MQ=35)
          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                        <  1:3320500/68‑1 (MQ=25)
           gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAc                                                       <  1:1403667/68‑1 (MQ=35)
            aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                      <  1:3125562/68‑1 (MQ=25)
            aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACg                                                      <  1:2184412/68‑1 (MQ=35)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                    <  1:2672109/68‑1 (MQ=25)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                    >  1:435734/1‑68 (MQ=35)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                    >  1:2852414/1‑68 (MQ=35)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                    >  1:2429603/1‑68 (MQ=35)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                    >  1:1525440/1‑68 (MQ=35)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCAGATGCGACGCt                                                    >  1:2694579/1‑68 (MQ=25)
                aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTg                                                  <  1:3082095/68‑1 (MQ=255)
                aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                  <  1:3247712/68‑1 (MQ=37)
                  cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                >  1:1591498/1‑68 (MQ=255)
                     gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                             <  1:661868/68‑1 (MQ=255)
                     gATTTTGTAGGGCGGATAAGGCGTTTGCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                             <  1:652252/68‑1 (MQ=14)
                          tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                        <  1:1113843/68‑1 (MQ=255)
                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                     >  1:505078/1‑68 (MQ=255)
                               gtgggATAGGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:1911071/65‑1 (MQ=14)
                               gTCGTATAGGGCATTTCCGCCCCATAAGACGTAAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:2418971/68‑1 (MQ=9)
                               gTCGGATAAGGCTTTTACGCCGCGTCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:2473315/68‑1 (MQ=14)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:4209/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:3274267/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:408857/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:434148/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:635097/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:2918239/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:1883966/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   >  1:3124689/1‑68 (MQ=255)
                               acaggATAAGCCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:845125/65‑1 (MQ=255)
                                    aTAAGGGATTTCCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                              <  1:1950087/68‑1 (MQ=14)
                                    aTAAGGCGTTTGCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:629178/68‑2 (MQ=14)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:1009583/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:821034/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:1174596/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:1706529/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:1963377/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                              <  1:2728361/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                              <  1:217136/68‑1 (MQ=255)
                                    aTAAGGCGATTACGCGGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:237028/68‑2 (MQ=14)
                                    aTAAGGCGATTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:2521195/68‑2 (MQ=14)
                                     tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCNCGTCTTATCAGGCCTatc                             <  1:2572111/68‑3 (MQ=14)
                                      aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                            >  1:1817892/1‑65 (MQ=255)
                                        ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                          >  1:741836/1‑68 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACTCCAGa                        >  1:144718/1‑68 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                        >  1:3272181/1‑68 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                        >  1:265992/1‑68 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCGAccccggc                        >  1:726006/1‑65 (MQ=14)
                                            tttCCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                      <  1:1770050/68‑1 (MQ=255)
                                            tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                      <  1:2832746/68‑1 (MQ=255)
                                              tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGActgt                    >  1:3271271/1‑65 (MQ=255)
                                              gccgccgcATCCGACATAAACGCGTGATGCGACGATTGCCTCGTCTTTTCAGGCCTACGCCAGACAgc                    <  1:3412787/66‑1 (MQ=14)
                                               aCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                   >  1:701122/1‑68 (MQ=255)
                                               aCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                   >  1:1340680/1‑68 (MQ=255)
                                               aCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                   >  1:334803/1‑68 (MQ=255)
                                                cgccgcATCCGACATCGCCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcgc                  <  1:2942073/68‑1 (MQ=255)
                                                      aTCGGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:1783423/68‑1 (MQ=255)
                                                      aTCCGACATCACCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:1785392/68‑1 (MQ=255)
                                                      aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:14270/68‑1 (MQ=40)
                                                      aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:3379545/68‑1 (MQ=40)
                                                      aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:3045000/68‑1 (MQ=40)
                                                      aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:2215542/68‑1 (MQ=40)
                                                            cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt      >  1:3143094/1‑68 (MQ=40)
                                                              tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa    >  1:1524285/1‑68 (MQ=39)
                                                                aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc  <  1:1921615/68‑1 (MQ=35)
                                                                  |                                                                 
CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC  >  NC_000913/4295994‑4296125

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