Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,983,761 C→T 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.3% 193.4 / ‑4.1 58E437K (GAA→AAA) araGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (27/30);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

ATCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983664‑1983856
                                                                                                 |                                                                                               
aTCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGt                                                                                               <  2:566901/100‑1 (MQ=255)
   cGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCt                                                                                            >  1:129720/1‑100 (MQ=255)
        ggCAACGCCGATGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCTCCGCCAGCGCATAACTTACGTTATATATTTTGTGCt                                                                                            <  2:641673/95‑1 (MQ=255)
         gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      <  1:355561/100‑1 (MQ=255)
          cAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGcc                                                                                     <  2:1296834/100‑1 (MQ=255)
             cGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                  <  1:356051/100‑1 (MQ=255)
                cGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGAGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATc                                                                               <  2:714732/100‑1 (MQ=255)
                  aGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCaa                                                                             <  2:660829/100‑1 (MQ=255)
                      cTTCAGGTAAGTCGCTGGAGGCAAACATCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGcc                                                                         <  2:1297266/100‑1 (MQ=255)
                      cTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGcc                                                                         >  2:322303/1‑100 (MQ=255)
                      cTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGcc                                                                         >  2:848206/1‑100 (MQ=255)
                         ctggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                      <  2:209977/98‑1 (MQ=255)
                         cAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                      >  1:985570/1‑100 (MQ=255)
                         cAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                      >  1:278212/1‑100 (MQ=255)
                         cAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                      >  2:486908/1‑100 (MQ=255)
                           ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                    <  1:226230/100‑1 (MQ=255)
                           ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgag                                                                    >  1:209979/1‑98 (MQ=255)
                            gTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGt                                                                   >  1:408756/1‑100 (MQ=255)
                             tAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTa                                                                  <  1:162060/100‑1 (MQ=255)
                              aaGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAg                                                                 <  2:222650/100‑1 (MQ=255)
                               aGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAgg                                                                <  1:932303/100‑1 (MQ=255)
                                gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                               >  1:973616/1‑100 (MQ=255)
                                 ncGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGtt                                                              >  2:35961/2‑100 (MQ=255)
                                   gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTCCGCTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTc                                                             >  1:280618/1‑99 (MQ=255)
                                   gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTc                                                             <  2:233086/99‑1 (MQ=255)
                                   gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                            >  2:572566/1‑100 (MQ=255)
                                    cTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCAt                                                           <  2:1212102/100‑1 (MQ=255)
                                      ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCAtag                                                         >  1:1212113/1‑98 (MQ=255)
                                      ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         <  1:53835/100‑1 (MQ=255)
                                      ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         >  1:782335/1‑100 (MQ=255)
                                      ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         >  2:430262/1‑100 (MQ=255)
                                       gAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCa                                                        >  2:112901/1‑100 (MQ=255)
                                        aGGCCAACAGCTCCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCTAACATCAATGCCGCGCGTAGGTTCACCCAg                                                       <  2:951840/100‑1 (MQ=255)
                                            aaaCAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaaa                                                   >  1:672595/1‑100 (MQ=255)
                                             aaCAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaaaa                                                  <  1:350032/100‑1 (MQ=255)
                                                aGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGa                                               <  2:291396/100‑1 (MQ=255)
                                                  cACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAcc                                             <  2:985562/100‑1 (MQ=255)
                                                  cACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATAGTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAcc                                             <  1:690226/100‑1 (MQ=255)
                                                   aCCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCt                                            >  2:332789/1‑100 (MQ=255)
                                                    ccgccccgccCTGCGCCGCCAGCGCATAAATTACGTTATATACTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGCCCt                                            >  2:506031/1‑99 (MQ=255)
                                                    ccGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCtt                                           >  1:743582/1‑100 (MQ=255)
                                                    ccGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCt                                            >  2:450891/1‑99 (MQ=255)
                                                     cGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATAGATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTc                                          <  2:285071/100‑1 (MQ=255)
                                                        cACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtc                                       >  1:32594/1‑100 (MQ=255)
                                                        cACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtc                                       <  1:35962/100‑1 (MQ=255)
                                                              cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 >  2:1301653/1‑100 (MQ=255)
                                                                  gccgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATa                             >  1:912662/1‑100 (MQ=255)
                                                                            cATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGcc                   <  2:1100969/100‑1 (MQ=255)
                                                                             aTAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGccc                  <  1:486912/100‑1 (MQ=255)
                                                                               aaaTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAg                <  2:672591/100‑1 (MQ=255)
                                                                               aaaTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAg                >  2:949009/1‑100 (MQ=255)
                                                                                     cGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGc          <  1:1181905/100‑1 (MQ=255)
                                                                                      gTTATATATTTTGTGCTTAGCGCCCACATCAATGCCGCTCGTAAGTTCATCCAGCAAAATGCCCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGc          >  2:637838/1‑99 (MQ=255)
                                                                                         atatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCtttt      <  1:430266/100‑1 (MQ=255)
                                                                                            tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc   <  1:710233/100‑1 (MQ=255)
                                                                                            tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc   <  2:938776/100‑1 (MQ=255)
                                                                                            tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc   <  1:1014945/100‑1 (MQ=255)
                                                                                             atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCGCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  >  1:1243226/1‑100 (MQ=255)
                                                                                                 |                                                                                               
ATCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983664‑1983856

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