Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,243,370 C→T 55.5% L49L (CTG→CTA malG ← maltose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,243,3700CT55.5% 19.0 / 91.2 75L49L (CTG→CTAmalGmaltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (14/19);  new base T (19/22);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

AGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGTAGCAGTAGCAGGTG  >  NC_000913/4243272‑4243468
                                                                                                  |                                                                                                  
agTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCa                                                                                                   <  2:2142343/100‑1 (MQ=255)
 gTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAg                                                                                                  <  1:856126/100‑1 (MQ=255)
   aCCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACTCCAGTTTCCATTGATCCCAGGAGATTTGTTCCGGGATCAgcc                                                                                                <  1:1512976/100‑3 (MQ=255)
   aCCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGc                                                                                                <  1:437663/100‑1 (MQ=255)
      gggAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGc                                                                                             <  1:1242928/100‑1 (MQ=255)
      gggAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGc                                                                                             >  1:1978393/1‑100 (MQ=255)
       ggAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGcc                                                                                            <  1:1347004/100‑1 (MQ=255)
       ggAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGcc                                                                                            <  1:1397795/100‑1 (MQ=255)
       ggAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGcc                                                                                            <  1:760907/100‑1 (MQ=255)
            ggCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCg                                                                                       <  2:1611234/100‑1 (MQ=255)
            ggCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCg                                                                                       <  1:1907119/100‑1 (MQ=255)
             gCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGc                                                                                      >  1:664410/1‑100 (MQ=255)
              cGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCa                                                                                     <  2:677438/100‑1 (MQ=255)
                   gCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGtt                                                                                <  1:293907/100‑1 (MQ=255)
                    cGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGttt                                                                               >  2:1955510/1‑100 (MQ=255)
                     gTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTc                                                                              <  2:145038/100‑1 (MQ=255)
                      tAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCTGGATCAGGCTGCCGGTCGCAAAGTTTcc                                                                             <  1:847552/100‑1 (MQ=255)
                       aaTGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTccc                                                                            >  2:895803/1‑100 (MQ=255)
                       aaTGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTccc                                                                            >  2:685291/1‑100 (MQ=255)
                        aTGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCt                                                                           <  1:1095687/100‑1 (MQ=255)
                          gCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGa                                                                         <  1:871926/100‑1 (MQ=255)
                          gCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGa                                                                         >  2:1280476/1‑100 (MQ=255)
                           cGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGAc                                                                        >  1:2213633/1‑100 (MQ=255)
                               cATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAg                                                                    >  2:1344748/1‑100 (MQ=255)
                                    gCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCgaga                                                               >  1:162773/1‑100 (MQ=255)
                                    gCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCgag                                                                <  2:310449/99‑1 (MQ=255)
                                    gCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCgaga                                                               >  2:493480/1‑100 (MQ=255)
                                      cTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATa                                                             >  2:1885298/1‑100 (MQ=255)
                                       tGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGAt                                                              <  2:765230/98‑1 (MQ=255)
                                        gtacAACGCTAAAGCCTACCGCCAGTTTCCAGTGATCCCAGGAGATGTGCTCGGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAg                                                            <  2:208592/96‑1 (MQ=255)
                                        gTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGACTTGCTCCGGGATGAGGCTGCCGGCCGCAACGTTTCCCTGACGCAGCGAGATa                                                             >  2:13979/1‑98 (MQ=255)
                                          tCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGcga                                                         >  1:572263/1‑100 (MQ=255)
                                           cAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGTAGATTTGCTCCGGGATCAGGCTGCCGGTCACAAAGTTTCCCTGACGCAGCGAGATAGcgac                                                        <  2:243789/100‑1 (MQ=255)
                                           cAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGcgac                                                        >  2:58565/1‑100 (MQ=255)
                                           cAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGcgac                                                        >  2:535148/1‑100 (MQ=255)
                                               gCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGAcc                                                    <  2:716829/100‑1 (MQ=255)
                                                  aaaaCCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATc                                                 >  1:286492/1‑100 (MQ=255)
                                                  aaaaCCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATc                                                 >  1:195964/1‑100 (MQ=255)
                                                  aaaaCCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATc                                                 <  2:1350931/100‑1 (MQ=255)
                                                     agcTAACGCCGGTTTCCAGTGATTCCAGGAGATTTGCTCCAGGATCATGCTGCCTGTCGCAAAGTTTCCCTGACGCAGTGATATAGCGACGACCATcag                                               <  2:852756/97‑1 (MQ=255)
                                                      ccTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATcagca                                             >  1:1258148/1‑100 (MQ=255)
                                                      ccTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATcagca                                             >  1:2241844/1‑100 (MQ=255)
                                                       cTAACGCCAGTTTCCAGTGATCGAAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCGGCGAGATAGCGACGACCATcagcag                                            <  1:1205253/100‑1 (MQ=255)
                                                        tAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCagcagc                                           <  2:396538/100‑1 (MQ=255)
                                                         aaCGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCATCg                                          <  2:1016488/100‑1 (MQ=255)
                                                            gCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGa                                       >  1:1299057/1‑100 (MQ=255)
                                                            gCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGa                                       <  2:650275/100‑1 (MQ=255)
                                                               aGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCATCGAGATAGCGACGACCATCAGCAGCGGGAAc                                     <  1:1388241/99‑1 (MQ=255)
                                                               aGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACa                                    >  1:2247018/1‑100 (MQ=255)
                                                                    ccAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATc                               <  2:96992/100‑1 (MQ=255)
                                                                     naGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCg                              >  1:723528/2‑100 (MQ=255)
                                                                     cAGTGATCCCAGGAGATTTGCTCCGGGATTCGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCg                              >  1:1385977/1‑100 (MQ=255)
                                                                     cAGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCg                              >  2:1148407/1‑100 (MQ=255)
                                                                      aGTGATCCCAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGc                             >  1:687137/1‑100 (MQ=255)
                                                                      aGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGc                             >  2:1994001/1‑100 (MQ=255)
                                                                          aTCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCg                         >  1:430761/1‑100 (MQ=255)
                                                                          aTCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCg                         <  1:1980510/100‑1 (MQ=255)
                                                                           ttccAGTAGATTTGCTCCGGGATTAGTCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCg                         <  1:92814/97‑1 (MQ=255)
                                                                            cccAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGAt                       >  1:49272/1‑100 (MQ=255)
                                                                            cccAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGAt                       <  1:183822/100‑1 (MQ=255)
                                                                             ccAGGAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATa                      >  1:224772/1‑100 (MQ=255)
                                                                                ggAGATTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaaa                   <  2:49272/100‑1 (MQ=255)
                                                                                ggAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaaa                   >  1:1606329/1‑100 (MQ=255)
                                                                                 gagaTTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaaaa                  <  2:1268232/100‑1 (MQ=255)
                                                                                   gaTTTGCTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagt                <  1:67486/100‑1 (MQ=255)
                                                                                   gaTTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagt                <  2:28311/100‑1 (MQ=255)
                                                                                   gaTTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagt                >  2:1422819/1‑100 (MQ=255)
                                                                                       tGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagca            >  2:1323322/1‑100 (MQ=255)
                                                                                         cTCCGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagcagt          >  2:1152939/1‑100 (MQ=255)
                                                                                         cTCCGGGATCAGGCTGCCGGTCGCAAAGTTTGCCTGACGCAGCGAGATAGCGACGACAATCAGCAGCGGGAACATAATCGCTGCGAGAAAAag                 >  2:1894634/1‑93 (MQ=255)
                                                                                         cTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagcagt          <  1:1721400/100‑1 (MQ=255)
                                                                                            cGGGATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagcagtagc       <  2:1960623/100‑1 (MQ=255)
                                                                                            cGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagcagtagc       <  1:422531/100‑1 (MQ=255)
                                                                                              ggATTAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGtagcagtagcag     <  1:895804/100‑1 (MQ=255)
                                                                                                 tCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGTAGCAGTAGCAGGTg  <  1:501662/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
AGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGTAGCAGTAGCAGGTG  >  NC_000913/4243272‑4243468

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: