Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT40.9% 42.4 / 87.6 66intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (16/23);  new base T (9/18);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA  >  NC_000913/4295994‑4296123
                                                                  |                                                               
cgcAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACAGCAACg                                                                <  2:1175116/68‑1 (MQ=255)
cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                <  2:697179/68‑1 (MQ=35)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                               >  1:923638/1‑68 (MQ=25)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                               >  2:1155363/1‑68 (MQ=25)
 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                               <  2:1068102/68‑1 (MQ=35)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                          <  1:521545/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                          <  1:494941/68‑1 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                          >  2:521543/1‑68 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                          >  2:494939/1‑68 (MQ=25)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                          >  2:528463/1‑68 (MQ=35)
      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                          >  2:483586/1‑68 (MQ=35)
      aGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                          <  1:250551/68‑1 (MQ=255)
        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                        <  1:156236/68‑1 (MQ=25)
        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                        >  2:1045965/1‑68 (MQ=35)
         cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                       >  1:218859/1‑68 (MQ=25)
             tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGc                                                   >  1:1256086/1‑68 (MQ=35)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                  >  2:259390/1‑68 (MQ=25)
              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                  >  2:89298/1‑68 (MQ=35)
                aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTg                                                <  1:367886/68‑1 (MQ=255)
                   gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                             >  1:49323/1‑68 (MQ=255)
                     gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                           <  1:1304995/68‑1 (MQ=255)
                        tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                        <  1:1052281/68‑1 (MQ=255)
                        tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCATGCCGCGTCt                                        >  2:951630/1‑68 (MQ=255)
                         ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt                                       <  2:1276567/68‑1 (MQ=255)
                         ttGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt                                       <  1:594168/68‑1 (MQ=255)
                          tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                      <  2:1099315/68‑1 (MQ=255)
                             gggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                   >  2:1192743/2‑68 (MQ=255)
                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                        >  2:1053750/1‑63 (MQ=255)
                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                        <  1:1053751/63‑1 (MQ=255)
                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                   >  1:309495/1‑68 (MQ=255)
                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                   >  1:127547/1‑68 (MQ=255)
                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                   >  2:1227979/1‑68 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  2:742968/68‑1 (MQ=255)
                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                 <  1:1211663/68‑1 (MQ=255)
                                    aTAAGTCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:153786/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                            <  2:969038/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                            <  1:856985/68‑2 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:483588/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  2:146423/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  2:1157939/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  2:70677/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:973/68‑1 (MQ=255)
                                    aTAAGGCGTTTACGCCGCAGCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:521975/68‑1 (MQ=255)
                                    aTAAGGCGCTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  2:976381/68‑1 (MQ=255)
                                     tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                           <  1:1155364/68‑3 (MQ=255)
                                     tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            >  2:965362/1‑67 (MQ=255)
                                     tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                            <  1:965363/67‑1 (MQ=255)
                                     tAAGGCGTTTACGCCGCATCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                           <  2:654152/68‑3 (MQ=14)
                                     tAAGGCGTTTACGCCGCATCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                           <  2:1099273/68‑3 (MQ=14)
                                     tAAGGCGGTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                           <  2:141776/68‑3 (MQ=14)
                                      aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                          <  1:394267/68‑4 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                          <  1:1332294/48‑1 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                          >  2:1332293/1‑48 (MQ=255)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACg                                                    <  2:1144488/38‑1 (MQ=40)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACg                                                    >  1:1144491/1‑38 (MQ=40)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACg                                                    >  1:1144489/1‑38 (MQ=40)
                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGa                      >  2:1242654/1‑68 (MQ=25)
                                            tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                    <  1:179427/68‑1 (MQ=255)
                                              tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                  <  2:1016903/68‑1 (MQ=255)
                                               aCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                 >  2:410769/1‑68 (MQ=255)
                                                cgccgcATCCGACTTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcgc                <  2:1058243/68‑1 (MQ=255)
                                                cgccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                          <  1:283875/58‑1 (MQ=255)
                                                cgccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCGACGCCAGACAgcgc                <  2:996330/68‑1 (MQ=255)
                                                      aTCCGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc          <  1:621797/68‑1 (MQ=255)
                                                          taCATCAACGCCTGATGCGACTCTTGCCGCGTCTTATCTGGCCTACTCCAGCCAGCGCAATAGCCtta      <  2:757819/67‑3 (MQ=14)
                                                              tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa  >  1:920891/1‑68 (MQ=39)
                                                                  |                                                               
CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA  >  NC_000913/4295994‑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: