Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT42.1% 39.2 / 95.3 76intergenic (+266/+376)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (29/15);  new base T (10/22);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 4.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC  >  NC_000913/4295994‑4296125
                                                                  |                                                                 
cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                  <  1:737703/68‑1 (MQ=35)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                 >  1:1628998/1‑68 (MQ=25)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                 <  1:1148099/68‑1 (MQ=25)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                 <  1:3838362/68‑1 (MQ=25)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                 >  1:2102398/1‑68 (MQ=35)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                 >  1:1032900/1‑68 (MQ=35)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                 >  1:3880662/1‑68 (MQ=35)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                 >  1:1248328/1‑68 (MQ=35)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            >  1:3182624/1‑68 (MQ=255)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACGTCAACGCCTGAt                                                            >  1:3591649/1‑68 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            <  1:991252/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            <  1:3450314/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                            <  1:3510723/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                            >  1:1385103/1‑68 (MQ=35)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                            >  1:2693718/1‑68 (MQ=35)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                            <  1:2241031/68‑1 (MQ=35)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                            >  1:1854740/1‑68 (MQ=35)
        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACCTCAACGCCTGATGc                                                          >  1:3685066/1‑68 (MQ=25)
        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                          >  1:1154797/1‑68 (MQ=25)
        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                          >  1:1256483/1‑68 (MQ=35)
        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                          >  1:4197414/1‑68 (MQ=35)
        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                          >  1:1752663/1‑68 (MQ=35)
         cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                         >  1:4010640/1‑68 (MQ=25)
         cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                         >  1:1139257/1‑68 (MQ=25)
         cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                         >  1:2127554/1‑68 (MQ=35)
           gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                       <  1:2356558/68‑1 (MQ=25)
            aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                      <  1:477977/68‑1 (MQ=25)
            aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCAGACATCAACGCCTGATGCGACg                                                      <  1:2784183/68‑1 (MQ=25)
              ttAGCGTGATTTTGTAGGTCGGATAANGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                    <  1:3256041/68‑1 (MQ=11)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                    >  1:936859/1‑68 (MQ=25)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                    >  1:1736264/1‑68 (MQ=35)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                    >  1:3853954/1‑68 (MQ=35)
                  cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                >  1:2190116/1‑68 (MQ=255)
                  cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                >  1:37320/1‑68 (MQ=40)
                     gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                             <  1:3881000/68‑1 (MQ=255)
                      aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                            >  1:1412240/1‑68 (MQ=255)
                        tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt                                          >  1:177275/1‑68 (MQ=255)
                         ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                         <  1:3779051/68‑1 (MQ=255)
                          tggaGGTCGGATAAGGCGTTTACGCCGCGTCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                        <  1:3105645/65‑1 (MQ=14)
                          tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                        <  1:1352081/68‑1 (MQ=255)
                          tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                        <  1:3985597/68‑1 (MQ=255)
                          tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                        <  1:1040338/68‑1 (MQ=255)
                          tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                        <  1:1524219/68‑1 (MQ=255)
                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                     >  1:3076521/1‑68 (MQ=255)
                               ggcGGATAAGGGGTTTACGCCGCGTCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:965542/66‑1 (MQ=14)
                               gTCGTATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:1821093/68‑1 (MQ=255)
                               gTCGGGTAAGGCGTTTCCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:3106701/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTGCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:534109/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:2238617/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   >  1:4154710/1‑68 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:1897733/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:1922756/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:1369054/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   >  1:2961365/1‑68 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTCATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:1723090/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:3300078/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:922660/68‑1 (MQ=255)
                               acaggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                   <  1:595467/65‑1 (MQ=255)
                                tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                  <  1:727626/68‑1 (MQ=255)
                                gcGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                  <  1:221146/67‑1 (MQ=255)
                                   gATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                               <  1:2364534/68‑1 (MQ=255)
                                    gcaaGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                              <  1:1946039/66‑1 (MQ=255)
                                    aTAATGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                              <  1:3832207/68‑1 (MQ=255)
                                    aTAAGGCGTTTCCGTCGCAGCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:2900065/68‑2 (MQ=11)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:818447/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                              <  1:2170820/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                              <  1:3875702/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                              <  1:183271/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                              <  1:3809250/68‑1 (MQ=255)
                                     tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                             <  1:3515844/68‑3 (MQ=255)
                                     tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                             <  1:132654/68‑3 (MQ=255)
                                     tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                             >  1:1865450/1‑68 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAgc                        >  1:3050045/1‑67 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAgc                        >  1:1150505/1‑67 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                        >  1:1931393/1‑68 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                        >  1:3671419/1‑68 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                        >  1:1848335/1‑68 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                        >  1:1207061/1‑68 (MQ=255)
                                                 gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCa                 >  1:643570/1‑68 (MQ=40)
                                                 gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCa                 >  1:1374219/1‑68 (MQ=40)
                                                 gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTACGCCAGACGGCGCa                 >  1:1997597/1‑68 (MQ=255)
                                                 gccgcACCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCa                 <  1:1107941/68‑1 (MQ=255)
                                                     cATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg             <  1:1795756/68‑1 (MQ=40)
                                                     cATCCGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg             <  1:1417550/68‑1 (MQ=255)
                                                      ctCCGACGTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:1211495/67‑1 (MQ=255)
                                                      aTCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:2903703/68‑1 (MQ=255)
                                                      aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:1660186/68‑1 (MQ=40)
                                                      aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc            <  1:3368001/68‑1 (MQ=40)
                                                          gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGa        <  1:2290044/68‑1 (MQ=40)
                                                           aCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTCAt       >  1:3744892/1‑68 (MQ=255)
                                                            cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt      >  1:276462/1‑68 (MQ=40)
                                                             aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt     >  1:371493/1‑68 (MQ=40)
                                                              tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa    >  1:1673681/1‑68 (MQ=39)
                                                              tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa    >  1:2324660/1‑68 (MQ=39)
                                                                agcGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc  <  1:1362626/66‑1 (MQ=14)
                                                                  |                                                                 
CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC  >  NC_000913/4295994‑4296125

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