Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,838,044 +G coding (498/1458 nt) ascF → PTS system enzyme IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,838,0431.G99.3% 374.6 / ‑3.6 135coding (497/1458 nt)ascFPTS system enzyme IIBC
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (59/75);  total (59/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 155 total aligned reads displayed.

CAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATC  >  NC_000913/2837968‑2838114
                                                                            |                                                                        
cAAAAGGTGCGCGGACCTTAACCATTCTGAAGGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                           <  2:5411426/76‑1 (MQ=255)
cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                           <  2:1858360/76‑1 (MQ=255)
cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                           <  1:947981/76‑1 (MQ=255)
cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                           <  2:1553418/76‑1 (MQ=255)
cAAAAGGATCGCCGACCTTAACCAGGCTGAATGTGATTGGTGACGGTACTTTATTCTTCCTGCCGCTCATGGTcgc                                                                           >  2:5721843/1‑76 (MQ=255)
 aaaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcg                                                                         <  2:2709471/76‑1 (MQ=255)
 aaaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcg                                                                         <  2:4921886/76‑1 (MQ=255)
  aaaGGATGTCCGGCCTTAACCATTCTGAAGGGGGTGGGGGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcgg                                                                        <  2:5807691/76‑2 (MQ=255)
  aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcgg                                                                        <  1:3764089/76‑2 (MQ=255)
  aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcgg                                                                        <  1:2419636/76‑2 (MQ=255)
  aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATAGTcgcgg                                                                        >  2:1543990/1‑75 (MQ=255)
   aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcggc                                                                       <  1:4249232/76‑3 (MQ=255)
    aGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcggca                                                                      <  1:1153058/76‑4 (MQ=255)
    aGGATCGCCGAACCTAACCAGTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGGcgcggct                                                                      >  1:1585814/1‑73 (MQ=255)
       aTCGCCGACCTTACCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCt                                                                   <  2:2743327/76‑1 (MQ=255)
       aTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCt                                                                   <  1:5763768/76‑1 (MQ=255)
         cGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGc                                                                 <  1:4078889/76‑1 (MQ=255)
          gCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGcc                                                                >  1:4530278/1‑76 (MQ=255)
           ccGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcg                                                               >  2:2220765/1‑76 (MQ=255)
              aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgcca                                                            >  1:5513691/1‑74 (MQ=255)
              aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgcca                                                            >  2:2272862/1‑74 (MQ=255)
               ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgccat                                                           >  1:4516960/1‑73 (MQ=255)
                 ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCa                                                        <  1:1398161/76‑1 (MQ=255)
                  tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCaa                                                       <  1:4271503/76‑1 (MQ=255)
                   aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCaaa                                                      >  1:201312/1‑76 (MQ=255)
                    aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAt                                                     <  2:4279371/76‑1 (MQ=255)
                    aCCATTCTGAAGGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAt                                                     <  2:5687357/76‑1 (MQ=255)
                     ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAtt                                                    <  1:1635582/76‑1 (MQ=255)
                      cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAttt                                                   <  1:3106051/76‑1 (MQ=255)
                       aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa                                                  <  2:5513689/76‑1 (MQ=255)
                       aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa                                                  >  2:1285939/1‑76 (MQ=255)
                       aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa                                                  >  2:5000138/1‑76 (MQ=255)
                       aTTCTGAAGGTGATTGGGGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGTCATCTGCCGCCATCAAATTTa                                                  <  2:2591639/76‑1 (MQ=255)
                        ttCTGAATGTGATTGGTGAGGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGGATCTGCCGCCATCAAATTTaa                                                 <  2:5062261/76‑1 (MQ=255)
                        ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATAGTCGCGGCATCTGCCGCCATCAAATTTaa                                                 <  1:1543990/76‑1 (MQ=255)
                         tCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAtttgaa                                                <  1:3919457/76‑4 (MQ=255)
                         tCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTaaa                                                <  1:3171984/76‑1 (MQ=255)
                           tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAc                                              >  1:3177152/1‑76 (MQ=255)
                           tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAc                                              >  1:3166985/1‑76 (MQ=255)
                           tGAATGTGATTGGTGACGGTGCTTTCTTCTGCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAc                                              <  2:1338462/76‑1 (MQ=255)
                            gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTaaagcc                                             <  1:3960141/76‑4 (MQ=255)
                            gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAcc                                             >  1:1197635/1‑76 (MQ=255)
                             aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCa                                            <  1:3804303/76‑1 (MQ=255)
                             aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAACa                                            >  1:1798492/1‑76 (MQ=255)
                              aTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCaa                                           <  1:852149/76‑1 (MQ=255)
                                gtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAata                                         >  1:1328641/1‑76 (MQ=255)
                                gtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAACAata                                         >  1:73580/1‑76 (MQ=255)
                                 tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAatat                                        <  1:1860781/76‑1 (MQ=255)
                                   atgggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGt                                      <  1:146098/73‑1 (MQ=255)
                                   aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGt                                      >  1:3211076/1‑76 (MQ=255)
                                   aTTGGTGAAGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGt                                      >  2:5468933/1‑76 (MQ=255)
                                    ttGGTGCCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTc                                     <  2:1925065/76‑1 (MQ=255)
                                    ttGGTGACGGTTCTTTCTTCTCCCGGCCGCTGAGGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTc                                     <  2:3079443/76‑1 (MQ=255)
                                     tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCg                                    <  2:4395943/76‑1 (MQ=255)
                                      ggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGc                                   >  1:1067961/1‑76 (MQ=255)
                                       gTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCt                                  <  2:2019847/76‑1 (MQ=255)
                                        tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTg                                 >  2:3371653/1‑76 (MQ=255)
                                        tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTg                                 <  2:4179156/76‑1 (MQ=255)
                                        ggccGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTg                                 <  2:3114954/73‑1 (MQ=255)
                                          aCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGc                               >  2:3985829/1‑76 (MQ=255)
                                           cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCg                              >  2:4449218/1‑76 (MQ=255)
                                            ggTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGATGGCGa                             >  1:2688552/1‑76 (MQ=255)
                                            ggTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCAGCAATCAAATGTAACACCAATATGTCGCTGGCgc                             >  1:2995540/1‑75 (MQ=255)
                                             gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGAt                            <  1:3639583/76‑1 (MQ=255)
                                             gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGAt                            <  1:5298233/76‑1 (MQ=255)
                                              gtgtttCTTCTTAGTGCCGCGGGTGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGAtt                           <  2:1958036/73‑1 (MQ=255)
                                               tcTTTCTCCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTg                          <  1:104342/75‑1 (MQ=255)
                                               gCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTg                          >  1:2162145/1‑76 (MQ=255)
                                                 tttCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGcc                        >  1:5370498/1‑76 (MQ=255)
                                                 tttCTTCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCTCTTGCGATTGcc                        <  1:446907/76‑1 (MQ=255)
                                                  ttcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                       <  1:3002075/76‑1 (MQ=255)
                                                   tcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                      <  2:2999405/76‑1 (MQ=255)
                                                    cttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                     <  1:2919283/76‑1 (MQ=255)
                                                    cttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                     <  1:2912314/76‑1 (MQ=255)
                                                    cttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                     <  1:846448/76‑1 (MQ=255)
                                                     ttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                    >  1:4970273/1‑76 (MQ=255)
                                                     ttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                    <  2:3568574/76‑1 (MQ=255)
                                                        ttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCCTCGCgg                 >  2:1883242/1‑76 (MQ=255)
                                                        ttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCgg                 <  1:2541677/76‑1 (MQ=255)
                                                        ttcCTGCCGCTGATGGTCGCGGCATCTGACGCCATCAAATTTAAAAACTATATGTCGCTGGCGCATGCCCGCGCgg                 >  1:1337948/1‑76 (MQ=255)
                                                          ccTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCgggg               >  1:402754/1‑75 (MQ=255)
                                                          ccTGCCGCTGATGGCCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggtt               <  1:5000139/76‑3 (MQ=255)
                                                           cTGCCGCTGATGTGCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg              <  1:2240174/76‑1 (MQ=255)
                                                           cTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg              <  1:937639/76‑1 (MQ=255)
                                                            tGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCTCTGGCGATTGCCATCGCGGgtgt             >  2:2183899/1‑76 (MQ=255)
                                                            tGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt             <  2:3474377/76‑1 (MQ=255)
                                                             gCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGtgtg            >  1:2085750/1‑76 (MQ=255)
                                                             gCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGtgtg            >  1:651206/1‑76 (MQ=255)
                                                              ccGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc           <  1:480037/76‑1 (MQ=255)
                                                               cGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt          >  2:2987462/1‑76 (MQ=255)
                                                                gCTGAGGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg         <  2:114416/76‑1 (MQ=255)
                                                                  tGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt       >  2:589615/1‑76 (MQ=255)
                                                                   gATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGtc      >  2:4461202/1‑75 (MQ=255)
                                                                   gATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa      >  1:3138044/1‑76 (MQ=255)
                                                                   gATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGAGTGCCATCGCGGGCGTGCTGGTa      >  2:4926238/1‑76 (MQ=255)
                                                                    aTGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc     >  1:5186034/1‑76 (MQ=255)
                                                                     tGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa    <  2:2085750/76‑1 (MQ=255)
                                                                     tGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa    >  2:782797/1‑76 (MQ=255)
                                                                     tGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTCGTACa    <  1:4903198/76‑1 (MQ=255)
                                                                       gTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc  <  2:2422419/76‑1 (MQ=255)
                                                                            |                                                                        
CAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATC  >  NC_000913/2837968‑2838114

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