Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,838,055 Δ1 bp coding (509/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,0550G.99.2% 593.0 / ‑3.9 133coding (509/1458 nt)ascFPTS system enzyme IIBC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (60/72);  minor base C (1/0);  total (61/72)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 133 total aligned reads displayed.

TTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATT  >  NC_000913/2837985‑2838125
                                                                        |                                                                      
ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCa                                                                   >  1:3568574/1‑76 (MQ=255)
ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCa                                                                   <  1:1398161/76‑1 (MQ=255)
 tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCaa                                                                  <  1:4271503/76‑1 (MQ=255)
  aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCaaa                                                                 >  2:151704/1‑76 (MQ=255)
  aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCaaa                                                                 >  1:201312/1‑76 (MQ=255)
   aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAt                                                                <  2:4279371/76‑1 (MQ=255)
   aCCATTCTGAAGGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAt                                                                <  2:5687357/76‑1 (MQ=255)
     cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAttt                                                              <  1:3106051/76‑1 (MQ=255)
      agtCTGAAGGTGATGTGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa                                                             <  2:1706952/74‑1 (MQ=255)
      aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa                                                             >  2:4186956/1‑76 (MQ=255)
      aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa                                                             <  1:2949752/76‑1 (MQ=255)
      aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa                                                             >  2:5000138/1‑76 (MQ=255)
      aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa                                                             <  2:5513689/76‑1 (MQ=255)
      aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa                                                             >  2:1285939/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:3857791/1‑76 (MQ=255)
          tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAc                                                         >  1:3166985/1‑76 (MQ=255)
          tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAc                                                         >  1:3177152/1‑76 (MQ=255)
          tGAATGCGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCCCGGCATCTGCCGCCCTCCAATTTCaaaa                                                          >  2:692442/1‑76 (MQ=255)
           gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTaaagcc                                                        <  1:3960141/76‑4 (MQ=255)
           gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAcc                                                        >  1:1197635/1‑76 (MQ=255)
            aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCa                                                       >  1:1975411/1‑76 (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)
                ggATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAatat                                                   <  2:1968547/75‑1 (MQ=255)
                  atgggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGt                                                 <  1:146098/73‑1 (MQ=255)
                  aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGt                                                 >  1:3211076/1‑76 (MQ=255)
                   ttGGTGCCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTc                                                <  2:1925065/76‑1 (MQ=255)
                   ttGGTGACGGTTCTTTCTTCTCCCGGCCGCTGAGGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTc                                                <  2:3079443/76‑1 (MQ=255)
                   ttGGTGACGGTGCTTTCTTCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTc                                                <  1:5742860/76‑1 (MQ=255)
                    tGGTGAGGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCg                                               <  2:4528061/76‑1 (MQ=255)
                    tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCg                                               <  2:4395943/76‑1 (MQ=255)
                    tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCg                                               <  1:5328120/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)
                       ggccGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTg                                            <  2:3114954/73‑1 (MQ=255)
                       ggACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTg                                            <  2:3841144/75‑1 (MQ=255)
                        gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTgc                                           >  2:2594692/1‑75 (MQ=255)
                         aCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGc                                          >  2:3985829/1‑76 (MQ=255)
                          cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCg                                         <  1:1352023/76‑1 (MQ=255)
                           ggTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGATGGCGa                                        >  1:2688552/1‑76 (MQ=255)
                           ggTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCAGCAATCAAATGTAACACCAATATGTCGCTGGCgc                                        >  1:2995540/1‑75 (MQ=255)
                            gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGAt                                       <  1:5298233/76‑1 (MQ=255)
                            gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGAt                                       <  1:3639583/76‑1 (MQ=255)
                            gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCCGGCGAt                                       >  2:2240174/1‑76 (MQ=255)
                             gtgtttCTTCTTAGTGCCGCGGGTGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                      <  2:1958036/73‑1 (MQ=255)
                              tcTTTCTCCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                     <  1:104342/75‑1 (MQ=255)
                              gCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                     >  1:2162145/1‑76 (MQ=255)
                                tttCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGcc                                   >  2:560817/1‑76 (MQ=255)
                                tttCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGcc                                   <  2:2303534/76‑1 (MQ=255)
                                 ttcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                  <  1:3002075/76‑1 (MQ=255)
                                  tcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                 <  2:2999405/76‑1 (MQ=255)
                                   cttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                <  1:846448/76‑1 (MQ=255)
                                   cttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                <  1:2912314/76‑1 (MQ=255)
                                    ttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                               <  2:3568574/76‑1 (MQ=255)
                                    ttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                               <  2:2974755/76‑1 (MQ=255)
                                    ttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                               >  1:4970273/1‑76 (MQ=255)
                                     tcttcCTGCCGCTGAAGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgc                              >  1:3173738/1‑76 (MQ=255)
                                      cttcCTGCCGCTGATGGTCGCGGCATATGCCGCCATCAAATTTAAAACCACTATGTCGCTGGCGATTGCCATcgcg                             <  1:641622/76‑1 (MQ=255)
                                       ttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCCTCGCgg                            >  2:1883242/1‑76 (MQ=255)
                                       ttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCgg                            >  1:1705983/1‑76 (MQ=255)
                                        tcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                           >  1:2708595/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)
                                           tGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCTCTGGCGATTGCCATCGCGGgtgt                        >  2:2183899/1‑76 (MQ=255)
                                            gCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGCTTGCCATCGCGGGtgtg                       >  2:342692/1‑76 (MQ=255)
                                            gCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGtgtg                       >  1:651206/1‑76 (MQ=255)
                                            gCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGtgtg                       >  1:2085750/1‑76 (MQ=255)
                                             ccGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                      <  1:480037/76‑1 (MQ=255)
                                              cGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTTGCTGGCGATTGCCATCGCGGGTGTGCt                     <  1:2399164/76‑1 (MQ=255)
                                              cGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt                     >  2:2987462/1‑76 (MQ=255)
                                               gCTGATGGTCGCGGCCTCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                    <  2:1326101/76‑1 (MQ=255)
                                               gCTGAGGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                    <  2:114416/76‑1 (MQ=255)
                                                 tGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                  >  2:589615/1‑76 (MQ=255)
                                                  gagggTCGCGGCCTCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                 <  2:5421724/73‑1 (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)
                                                         gcggcATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATTCGa          >  1:4319959/4‑76 (MQ=255)
                                                           ggcATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt        >  1:4783997/2‑76 (MQ=255)
                                                           ggcATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt        >  2:3235393/2‑76 (MQ=255)
                                                           ggcATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt        <  1:1621533/75‑1 (MQ=255)
                                                           ggcATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt        <  1:155021/75‑1 (MQ=255)
                                                              cATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGttt      >  1:3360025/1‑76 (MQ=255)
                                                                tCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGCGTGCTGATACATGCGCGtcttt    >  2:3460352/1‑75 (MQ=255)
                                                                 cTGCCGCCATCAAATTTAAGACCAATCTGTCGCTGGCGCTTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt   <  2:1999979/76‑1 (MQ=255)
                                                                 cTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt   >  1:5149946/1‑76 (MQ=255)
                                                                        |                                                                      
TTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATT  >  NC_000913/2837985‑2838125

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