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_000913	1,983,761	0	C	T	92.8%	211.1 / 4.9	69	E437K (GAA→AAA)	araG	L‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base T (23/41);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGC  >  NC_000913/1983666‑1983855                                                                                                |                                                                                               ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                            <  1:1077479/100‑1 (MQ=255) ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                            <  1:1498523/100‑1 (MQ=255)  cGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAACTACGTTATATATTTTGTGCt                                                                                           <  2:958634/100‑1 (MQ=255)   ggTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCtt                                                                                          >  2:206985/1‑100 (MQ=255)   ggTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCtt                                                                                          <  1:206986/100‑1 (MQ=255)   ggTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCtt                                                                                          >  1:1158869/1‑100 (MQ=255)   ggTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCtt                                                                                          >  1:1256422/1‑100 (MQ=255)   ggTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCtt                                                                                          <  2:806027/100‑1 (MQ=255)        gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                     >  2:1165622/1‑100 (MQ=255)        gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                     <  1:194514/100‑1 (MQ=255)            cGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                 <  1:1462517/100‑1 (MQ=255)            cGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                 <  2:1012716/100‑1 (MQ=255)             gCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACa                                                                                <  1:790662/100‑1 (MQ=255)             gCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACCTTATATAATTTGTGCTTAGCGCCAACa                                                                                >  2:948033/1‑100 (MQ=255)              ccGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACAt                                                                               >  2:1064903/1‑100 (MQ=255)               cGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATc                                                                              >  1:786334/1‑100 (MQ=255)               cGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATc                                                                              >  1:527305/1‑100 (MQ=255)                    aCTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGc                                                                         <  2:1026403/100‑1 (MQ=255)                       tCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgc                                                                      >  1:1135525/1‑100 (MQ=255)                       tCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgc                                                                      <  1:535084/100‑1 (MQ=255)                        cAGGTAAGTCGCTGGAGGCATACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                     >  2:349446/1‑100 (MQ=255)                        cAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                     <  1:1455134/100‑1 (MQ=255)                          ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                   <  2:1465864/100‑1 (MQ=255)                           gTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGt                                                                  <  2:471341/100‑1 (MQ=255)                            tAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTa                                                                 <  2:34149/100‑1 (MQ=255)                             aaGTCGCTGGAGGCAAACAGCCCCGCCACGCCCTGCGCCGCCAGCGCATACCCTACGTTCTCTATTTTGTGCTTAGCGCCAACAGCAATGCCGCGCGTAg                                                                <  2:521876/100‑1 (MQ=255)                              aGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAgg                                                               >  1:244696/1‑100 (MQ=255)                              aGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAgg                                                               >  2:114024/1‑100 (MQ=255)                               gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCCACATCAATGCCGCGCGTa                                                                 >  2:1573438/1‑97 (MQ=255)                                 cGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGTGTAGGTTc                                                            >  2:62564/1‑100 (MQ=255)                                 cGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTc                                                            <  2:1063029/100‑1 (MQ=255)                                  gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                           <  2:1331849/100‑1 (MQ=255)                                  gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                           <  2:1455674/100‑1 (MQ=255)                                    tGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATc                                                         <  1:1336704/100‑1 (MQ=255)                                     ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                        <  1:46139/100‑1 (MQ=255)                                     ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                        <  1:232061/100‑1 (MQ=255)                                      gAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCa                                                       >  2:982572/1‑100 (MQ=255)                                       aGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAg                                                      >  1:1301976/1‑100 (MQ=255)                                       aGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAg                                                      <  1:458624/100‑1 (MQ=255)                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                   >  2:854527/1‑100 (MQ=255)                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                   >  2:1293113/1‑100 (MQ=255)                                              cAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATg                                               >  1:777366/1‑100 (MQ=255)                                              cAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATg                                               <  2:438670/100‑1 (MQ=255)                                                gCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAc                                             >  2:157227/1‑100 (MQ=255)                                                  aCCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCt                                           <  2:33210/100‑1 (MQ=255)                                                   ccGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCtt                                          <  1:86221/100‑1 (MQ=255)                                                     gCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                        >  1:696643/1‑100 (MQ=255)                                                       cACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtc                                      <  1:82381/100‑1 (MQ=255)                                                       cACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtc                                      <  2:202056/100‑1 (MQ=255)                                                        aCGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtc                                      <  2:836101/99‑1 (MQ=255)                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                <  2:786329/100‑1 (MQ=255)                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                <  1:1013702/100‑1 (MQ=255)                                                               gcgcCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                <  2:79298/98‑1 (MQ=255)                                                               gcgcCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGa                              <  1:1069814/100‑1 (MQ=255)                                                                   cgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAc                          <  1:1145573/100‑1 (MQ=255)                                                                   cgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAc                          <  1:1146277/100‑1 (MQ=255)                                                                    gccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAcc                         >  2:122487/1‑100 (MQ=255)                                                                      cAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAg                       <  1:30976/100‑1 (MQ=255)                                                                         cgcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCACCTCTTCCGATAACCAGCgg                    <  1:1662/100‑1 (MQ=255)                                                                              aaaTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAg               <  1:603241/100‑1 (MQ=255)                                                                                aTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGaa             <  1:851078/100‑1 (MQ=255)                                                                                 ttACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAAt            >  2:801699/1‑100 (MQ=255)                                                                                     gTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCt        <  2:781034/100‑1 (MQ=255)                                                                                        atatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCtttt     <  1:854533/100‑1 (MQ=255)                                                                                         tatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCttttt    >  2:742134/1‑100 (MQ=255)                                                                                         tatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCttttt    <  2:1144104/100‑1 (MQ=255)                                                                                          atatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTg   <  1:122487/100‑1 (MQ=255)                                                                                           tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc  <  2:527302/100‑1 (MQ=255)                                                                                           tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc  <  1:1424058/100‑1 (MQ=255)                                                                                                |                                                                                               CCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGC  >  NC_000913/1983666‑1983855

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