Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT49.3% ‑2.8 / 112.8 73intergenic (+266/+376)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (24/13);  new base T (18/18);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA  >  NC_000913/4295996‑4296123
                                                                |                                                               
cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAc                                                                 >  1:817471/1‑65 (MQ=35)
 aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCGCCGCATCCGACATCAACg                                                                <  1:443696/65‑1 (MQ=25)
  aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                               <  1:3473610/65‑1 (MQ=25)
  aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                               >  1:2157945/1‑65 (MQ=35)
  aTAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                               <  1:1722377/65‑1 (MQ=255)
    aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                             >  1:3470913/1‑65 (MQ=25)
    aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                             >  1:3456833/1‑65 (MQ=25)
    aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                             >  1:2414844/1‑65 (MQ=25)
    aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                             >  1:2244382/1‑65 (MQ=25)
    aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                             >  1:334901/1‑65 (MQ=35)
    aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                             >  1:2176498/1‑65 (MQ=35)
      ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                           >  1:1341486/1‑65 (MQ=25)
      ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTgt                                                           >  1:126303/1‑64 (MQ=25)
      ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                           >  1:1266925/1‑65 (MQ=35)
      ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                           >  1:114823/1‑65 (MQ=35)
      ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                           >  1:3051480/1‑65 (MQ=35)
      ccTGATTTAGCGTGATTTTGTAGGTCCGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                           >  1:1068603/1‑65 (MQ=255)
       cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                          <  1:3476046/65‑1 (MQ=25)
       cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                          <  1:373057/65‑1 (MQ=35)
       cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                          >  1:913045/1‑65 (MQ=35)
       cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACACCGCATCCGACATCAACGCCTGAt                                                          <  1:663506/65‑1 (MQ=25)
       cTGATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                          <  1:2363009/65‑1 (MQ=255)
         gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                        >  1:616901/1‑65 (MQ=25)
           tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                      >  1:151821/1‑65 (MQ=25)
            ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                     >  1:3059198/1‑65 (MQ=25)
            ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                     >  1:1372118/1‑65 (MQ=25)
             tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACg                                                    <  1:391735/65‑1 (MQ=35)
                cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                 >  1:1644060/1‑65 (MQ=255)
                cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                 >  1:2905496/1‑65 (MQ=255)
                cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                 >  1:2340718/1‑65 (MQ=255)
                cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                 >  1:2734985/1‑65 (MQ=255)
                cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCtt                                                 >  1:3616357/1‑65 (MQ=35)
                 gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                <  1:3063487/65‑1 (MQ=37)
                 gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                <  1:78555/65‑1 (MQ=37)
                      tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                           <  1:2448149/65‑1 (MQ=255)
                      tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                           <  1:3457359/65‑1 (MQ=255)
                      tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                           <  1:2834983/65‑1 (MQ=255)
                      tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                           <  1:2382318/65‑1 (MQ=255)
                       ttGTAGGTGGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                          <  1:1086734/65‑1 (MQ=255)
                         gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt                                        <  1:158580/65‑1 (MQ=255)
                         gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt                                        <  1:1585024/65‑1 (MQ=255)
                         gTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                        <  1:3054923/65‑1 (MQ=255)
                           aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                      >  1:1848322/1‑65 (MQ=255)
                           aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                      >  1:665235/1‑65 (MQ=255)
                           aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                      <  1:841233/65‑1 (MQ=255)
                           aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                      >  1:2469855/1‑65 (MQ=255)
                                gggtAAGGCGTTTGCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  1:597647/62‑1 (MQ=14)
                                ggATAAGTCGTTTCCGCCGCATCCGACATCAATGCCTAATGCGACGCTTGCCGCGTCTTATCAgg                                 <  1:709920/65‑1 (MQ=14)
                                ggATAAGGCGTTTCCCCCTCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  1:329021/65‑1 (MQ=11)
                                ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  1:619147/65‑1 (MQ=255)
                                ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  1:722898/65‑1 (MQ=255)
                                ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  1:3545885/65‑1 (MQ=255)
                                ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 >  1:1653406/1‑65 (MQ=255)
                                ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  1:1771012/65‑1 (MQ=255)
                                ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  1:2065052/65‑1 (MQ=255)
                                 gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                <  1:1844172/65‑1 (MQ=255)
                                   tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGCCTTATCAGGCCt                              >  1:2294720/1‑65 (MQ=14)
                                    aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                             <  1:793782/65‑1 (MQ=255)
                                     gggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                            <  1:3260788/64‑2 (MQ=255)
                                     gggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:1405539/64‑1 (MQ=255)
                                     aGGCGTTTCCGCCGCATCCGCCATCAACGCCTGATGCGGCGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:1630990/65‑1 (MQ=14)
                                     aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:3048262/65‑1 (MQ=255)
                                      ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                           >  1:2927970/1‑63 (MQ=255)
                                      ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                           <  1:712276/65‑3 (MQ=255)
                                      ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                           <  1:1184290/65‑3 (MQ=255)
                                      ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                           <  1:1309054/65‑3 (MQ=255)
                                      ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                           <  1:336283/65‑3 (MQ=255)
                                      ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGGCTTTTCAGGCCTatc                           >  1:418236/1‑63 (MQ=14)
                                        cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                         >  1:3400232/1‑65 (MQ=255)
                                        cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                         >  1:289819/1‑65 (MQ=255)
                                        cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                         >  1:2544750/1‑65 (MQ=255)
                                        cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                         >  1:1355740/1‑65 (MQ=255)
                                             aCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                    >  1:199690/1‑65 (MQ=255)
                                             aCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGCCAGACa                    >  1:188579/1‑65 (MQ=255)
                                             aCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGCCCTACGCCAGACa                    >  1:208394/1‑65 (MQ=255)
                                              cgccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                   >  1:688070/1‑65 (MQ=255)
                                               gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                  >  1:664717/1‑65 (MQ=40)
                                               gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                  >  1:724205/1‑65 (MQ=40)
                                               gccccACCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                  <  1:1532981/65‑1 (MQ=255)
                                                 cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcgc                <  1:2274849/65‑1 (MQ=40)
                                                 cccATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcgc                <  1:2818385/63‑1 (MQ=255)
                                                       cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc          <  1:1008134/65‑1 (MQ=40)
                                                          cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTg       >  1:2397480/1‑65 (MQ=40)
                                                           aTCAGCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGa      <  1:2843067/65‑1 (MQ=25)
                                                           aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGa      >  1:2171905/1‑65 (MQ=40)
                                                             cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt    <  1:2386759/65‑1 (MQ=37)
                                                              aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt   >  1:1604735/1‑65 (MQ=35)
                                                               aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa  >  1:3196651/1‑65 (MQ=33)
                                                                |                                                               
CAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA  >  NC_000913/4295996‑4296123

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