Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,983,761 C→T 100% E437K (GAA→AAA)  araG ← L‑arabinose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,983,7610CT98.8% 255.8 / ‑5.1 77E437K (GAA→AAA) araGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (41/35);  total (42/35)
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.95e-01

TCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983665‑1983856
                                                                                                |                                                                                               
tctggTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGGGGCAAACAGCACCGCCACGCCCTTCGTCGCCAGCGCATAAGTTACGTTATATATTTTGt                                                                                               <  2:216673/96‑1 (MQ=255)
 ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             <  1:1981631/100‑1 (MQ=255)
 ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             <  1:1535613/100‑1 (MQ=255)
  cGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCt                                                                                            >  1:2273950/1‑100 (MQ=255)
  cGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCt                                                                                            >  1:554258/1‑100 (MQ=255)
    ttCGGCAACGCCGAGTACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             <  1:2003827/96‑1 (MQ=255)
    gTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACTGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCtt                                                                                           <  1:847320/99‑1 (MQ=255)
     tCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAg                                                                                         <  1:730137/100‑1 (MQ=255)
       ggCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcg                                                                                       >  1:1646170/1‑100 (MQ=255)
        gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      >  1:1401868/1‑100 (MQ=255)
         cAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGcc                                                                                     <  1:1043476/100‑1 (MQ=255)
         cAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGcc                                                                                     >  2:1410978/1‑100 (MQ=255)
          aaCGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCa                                                                                    <  2:1115857/100‑1 (MQ=255)
                gAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCa                                                                              <  1:2131745/100‑1 (MQ=255)
                  ggACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAAt                                                                            <  1:231496/100‑1 (MQ=255)
                  ggACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAAt                                                                            <  2:1608122/100‑1 (MQ=255)
                  ggACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAAt                                                                            <  1:970197/100‑1 (MQ=255)
                   gACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATg                                                                           <  2:2138067/100‑1 (MQ=255)
                       tCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCcgc                                                                       >  2:578996/1‑100 (MQ=255)
                        cAGTTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAAGGCcgc                                                                       <  2:486689/99‑1 (MQ=255)
                          ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                    >  1:1646567/1‑100 (MQ=255)
                          ggTAAGTAGCTGTAGGCACAAAGAGCTTCCTCGCCATGCGGCGCCAGCGCATATCTTCCGTTAGATCTTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                    <  2:1697076/100‑1 (MQ=255)
                             aaGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAg                                                                 >  1:940523/1‑100 (MQ=255)
                             aaGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAg                                                                 <  2:378651/100‑1 (MQ=255)
                              aGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAgg                                                                >  1:2180602/1‑100 (MQ=255)
                              aGTCGCTGGAGGCAAACAGCACCGCCAAGCCCTGCGCCGCCAGCGCATCAATTACGTTATATATATTGTGCTTAACGCAACCCTCAATGCTGCTCGTa                                                                  >  2:1098005/1‑98 (MQ=255)
                               gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                               >  2:1520659/1‑100 (MQ=255)
                               gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                               <  2:587652/100‑1 (MQ=255)
                                tCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGtt                                                              <  2:2047392/100‑1 (MQ=255)
                                 cGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTc                                                             >  2:631917/1‑100 (MQ=255)
                                 cGCTGGAGGCAAACAGCACAGCAACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGtt                                                              >  2:1495683/1‑99 (MQ=255)
                                  gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                            >  2:1797183/1‑100 (MQ=255)
                                     ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         <  1:1797204/100‑1 (MQ=255)
                                     ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         <  1:1581044/100‑1 (MQ=255)
                                     ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         <  2:885897/100‑1 (MQ=255)
                                        ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAg                                                       >  1:238807/1‑99 (MQ=255)
                                        ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGc                                                      >  2:492327/1‑100 (MQ=255)
                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                    >  1:1985392/1‑100 (MQ=255)
                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                    >  1:838360/1‑100 (MQ=255)
                                           aaaCAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaaa                                                   >  2:1028376/1‑100 (MQ=255)
                                              cAGCTCCGACACGACCTGCACCGCCAGCGCAAAAATAACGTGATATGTTGTTGGCTTAGCGCCAACATCGAGGCGGCGCGTGTGTTGATCCGGCAATATg                                                >  2:196608/1‑100 (MQ=255)
                                              cAGCACGGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAAt                                                 >  2:1997759/1‑99 (MQ=255)
                                                gCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAc                                              >  2:1426454/1‑100 (MQ=255)
                                                 cACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTCGCGGCAAAATGAATGCAGGGCGTAGGTTAATACGACAAAATGAct                                             >  2:2170844/1‑99 (MQ=255)
                                                 cACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAcc                                             >  1:1917020/1‑100 (MQ=255)
                                                    cgccccgcccTGCGCCGCCAGCGCATACATTACGTTATCTATTTTGTGCCTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTc                                          <  2:2396972/100‑1 (MQ=255)
                                                     gCCGCGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATAAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                         >  1:1556436/1‑100 (MQ=255)
                                                     gCCACGTCCTGCTCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                         <  2:119369/100‑1 (MQ=255)
                                                     gCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                         >  1:379298/1‑100 (MQ=255)
                                                        aCGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtct                                      >  2:775410/1‑100 (MQ=255)
                                                         cacccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtctc                                     >  1:1618652/3‑100 (MQ=255)
                                                           cccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCtt                                   >  2:1592966/1‑100 (MQ=255)
                                                            ccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTc                                  <  2:271161/100‑1 (MQ=255)
                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 <  1:1028387/100‑1 (MQ=255)
                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 <  2:1662207/100‑1 (MQ=255)
                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 <  1:659125/100‑1 (MQ=255)
                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 >  1:592738/1‑100 (MQ=255)
                                                                  ccgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATaa                            >  1:1314305/1‑100 (MQ=255)
                                                                  ccgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATaa                            <  1:1603657/100‑1 (MQ=255)
                                                                   cgccAGCGCATACATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAc                           >  1:1575303/1‑100 (MQ=255)
                                                                    gccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAcc                          <  1:1133025/100‑1 (MQ=255)
                                                                     ccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAcc                          <  1:2128941/99‑1 (MQ=255)
                                                                         cgcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCgg                     >  2:747933/1‑100 (MQ=255)
                                                                           cATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGcc                   >  2:2094227/1‑100 (MQ=255)
                                                                           cATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATACAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGcc                   >  2:1485987/1‑100 (MQ=255)
                                                                             tAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAGCCAGCGGCCCa                 >  2:875317/1‑100 (MQ=255)
                                                                             tAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCa                 <  2:1625593/100‑1 (MQ=255)
                                                                             tAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCa                 <  2:46418/100‑1 (MQ=255)
                                                                             tAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCa                 >  1:2145266/1‑100 (MQ=255)
                                                                              aaaTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAg                >  2:98212/1‑100 (MQ=255)
                                                                                 ttACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAAt             <  1:473882/100‑1 (MQ=255)
                                                                                  tACGTTATATCTTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATg            <  2:1575283/100‑1 (MQ=255)
                                                                                    cGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGc          >  2:331781/1‑100 (MQ=255)
                                                                                    cGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGc          >  1:2225988/1‑100 (MQ=255)
                                                                                         tatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCtttt      <  1:1349393/99‑1 (MQ=255)
                                                                                          atatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTg    >  1:1510928/1‑100 (MQ=255)
                                                                                            atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  >  2:885221/1‑100 (MQ=255)
                                                                                                |                                                                                               
TCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983665‑1983856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑