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.5% 239.3 / ‑4.3 68E437K (GAA→AAA) araGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/38);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983666‑1983856
                                                                                               |                                                                                               
ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             <  2:948033/100‑1 (MQ=255)
ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             <  1:1254783/100‑1 (MQ=255)
 cGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCt                                                                                            <  1:105692/100‑1 (MQ=255)
 cGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCt                                                                                            <  1:778758/100‑1 (MQ=255)
   gTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTa                                                                                          >  2:821055/1‑100 (MQ=255)
   gTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTa                                                                                          >  2:1271744/1‑100 (MQ=255)
       gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      >  1:1071239/1‑100 (MQ=255)
       gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      >  1:1086036/1‑100 (MQ=255)
       gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      >  2:886078/1‑100 (MQ=255)
       gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      >  1:696254/1‑100 (MQ=255)
         aaCGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCa                                                                                    <  1:821058/100‑1 (MQ=255)
           cGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTATGTTATATATTTTGTGCTTAGCGCCAAc                                                                                  <  1:29445/100‑1 (MQ=255)
            gCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACa                                                                                 >  1:72884/1‑100 (MQ=255)
            gCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACa                                                                                 <  1:1117052/100‑1 (MQ=255)
              cGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATc                                                                               <  1:1725538/100‑1 (MQ=255)
                 ggACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAAt                                                                            <  1:1418401/100‑1 (MQ=255)
                 ggACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAAt                                                                            <  1:998812/100‑1 (MQ=255)
                   aCTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGc                                                                          <  1:1742223/100‑1 (MQ=255)
                   aCTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGc                                                                          <  2:899966/100‑1 (MQ=255)
                    cTTCATGTAAGGCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGcc                                                                         <  1:641875/100‑1 (MQ=255)
                      tCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgc                                                                       >  1:1931825/1‑100 (MQ=255)
                         ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                    <  1:1271750/100‑1 (MQ=255)
                         ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                    <  1:185772/100‑1 (MQ=255)
                           tAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGTCGCCAGCGCATAAATTACGATATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTa                                                                  <  2:1115676/100‑1 (MQ=255)
                             aGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAgg                                                                <  1:1788007/100‑1 (MQ=255)
                              gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                               >  2:1365032/1‑100 (MQ=255)
                              gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCCTAGCGCCAACATCACTGCCGCGCGTAgg                                                                >  2:1203656/1‑99 (MQ=255)
                               tCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGtt                                                              <  2:1109605/100‑1 (MQ=255)
                                    ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         <  1:1444703/100‑1 (MQ=255)
                                    ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         >  1:1333849/1‑100 (MQ=255)
                                    ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         <  1:1821477/100‑1 (MQ=255)
                                    ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         >  2:175201/1‑100 (MQ=255)
                                      aGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAg                                                       <  2:1892509/100‑1 (MQ=255)
                                           aaCAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaaaa                                                  <  1:175202/100‑1 (MQ=255)
                                             cAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATg                                                >  2:779408/1‑100 (MQ=255)
                                             cAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATg                                                >  2:1326969/1‑100 (MQ=255)
                                                cACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAcc                                             >  1:409834/1‑100 (MQ=255)
                                                    gCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                         <  2:361438/100‑1 (MQ=255)
                                                    gCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                         <  1:1647001/100‑1 (MQ=255)
                                                     ccACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAt                                        >  2:451798/1‑100 (MQ=255)
                                                      cACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtc                                       >  1:860891/1‑100 (MQ=255)
                                                      cACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtc                                       >  1:1311243/1‑100 (MQ=255)
                                                          cccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCtt                                   >  1:1041189/1‑100 (MQ=255)
                                                          cccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCtt                                   <  2:1086031/100‑1 (MQ=255)
                                                           ccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTc                                  <  2:1452377/100‑1 (MQ=255)
                                                            cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 >  1:1457317/1‑100 (MQ=255)
                                                            cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 >  2:1920095/1‑100 (MQ=255)
                                                              gcgcCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGa                               <  1:208410/100‑1 (MQ=255)
                                                                gccgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATa                             >  2:808490/1‑100 (MQ=255)
                                                                   gccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAcc                          <  2:1333840/100‑1 (MQ=255)
                                                                   gccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAcc                          <  2:467159/100‑1 (MQ=255)
                                                                        cgcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCgg                     >  2:591614/1‑100 (MQ=255)
                                                                        cgcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCgg                     >  1:1018061/1‑100 (MQ=255)
                                                                          cATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGcc                   <  1:259636/100‑1 (MQ=255)
                                                                          cATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGcc                   <  2:860888/100‑1 (MQ=255)
                                                                           aTAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCGTCTAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGccc                  >  2:1432530/1‑100 (MQ=255)
                                                                            tAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGcccn                 <  2:1228156/100‑2 (MQ=255)
                                                                            tAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCa                 >  1:1518180/1‑100 (MQ=255)
                                                                            tAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCa                 >  1:342527/1‑100 (MQ=255)
                                                                               aTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACTAGCGGCCCAGaa              <  1:1341974/100‑1 (MQ=255)
                                                                               aTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGaa              >  2:1870854/1‑100 (MQ=255)
                                                                                    gTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCt         >  2:1128299/1‑100 (MQ=255)
                                                                                         atatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTg    <  2:228852/100‑1 (MQ=255)
                                                                                         atatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTg    <  2:1041184/100‑1 (MQ=255)
                                                                                         atatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTg    <  1:180993/100‑1 (MQ=255)
                                                                                          tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc   <  2:26092/100‑1 (MQ=255)
                                                                                           atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  >  2:695003/1‑100 (MQ=255)
                                                                                           atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  <  1:1365041/100‑1 (MQ=255)
                                                                                               |                                                                                               
CCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983666‑1983856

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