Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	49,855	T→G	28.7%	D11E (GAT→GAG)	folA →	dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,855	0	T	G	28.7%	84.7 / 59.4	70	D11E (GAT→GAG)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (21/29);  new base G (11/9);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

CACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGG  >  NC_000913/49758‑49949                                                                                                  |                                                                                               cACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATcg                                                                                              >  1:393167/1‑100 (MQ=255)   ccAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATcgcg                                                                                            <  2:1934810/100‑1 (MQ=255)   ccAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGcgcg                                                                                            >  1:223522/1‑100 (MQ=255)        cgacgaCGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATc                                                                                       >  2:682428/1‑100 (MQ=255)          acgacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCgg                                                                                     <  1:32473/100‑1 (MQ=255)          acgacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCgg                                                                                     <  1:918082/100‑1 (MQ=255)           cgacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATc                                                                                       <  1:2271404/97‑1 (MQ=255)            gacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCa                                                                                   <  2:395647/100‑1 (MQ=255)               ggTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATgg                                                                                <  2:1868419/100‑1 (MQ=255)                   tACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaaa                                                                            >  1:1962864/1‑100 (MQ=255)                     cGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAAc                                                                           >  2:2062440/1‑99 (MQ=255)                        tttACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                                                                       <  2:878700/100‑1 (MQ=255)                         ttACGTATAGTGGCGACAATTTTTTTTATTTTTTTATCTCAATGATCAGTCTGATTGAGGCGTTAGCGATAGATCTCGTTATCGGCATGCaacaaaccat                                                                      >  1:2367212/1‑92 (MQ=255)                         ttACGTATAGTGGCGACAATTTTTTTTATCGGTAAATCTCTATGATCTGACTCATTGCGGCTTTAGCAGTAGATCGCCTTATCTCCATTAAAAATGCCAt                                                                      >  1:1065668/1‑100 (MQ=255)                            cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGcc                                                                   >  1:709836/1‑100 (MQ=255)                               ataGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTg                                                                <  2:1419333/100‑1 (MQ=255)                               ataGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTg                                                                <  1:1009103/100‑1 (MQ=255)                                taGTGGCTACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTgg                                                               <  1:306934/100‑1 (MQ=255)                                taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTgg                                                               <  1:2480169/100‑1 (MQ=255)                                 aGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGa                                                              <  2:1474601/100‑1 (MQ=255)                                   tGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCTGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAc                                                            >  1:472130/1‑100 (MQ=255)                                    ggCGACGATTTTTTTTATCGGGAAATCTCAATGAGCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAc                                                            <  2:1121291/99‑1 (MQ=255)                                    ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcc                                                           <  1:113118/100‑1 (MQ=255)                                    ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcc                                                           >  1:110373/1‑100 (MQ=255)                                    ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcc                                                           >  2:730412/1‑100 (MQ=255)                                    ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcc                                                           >  1:498445/1‑100 (MQ=255)                                    ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcc                                                           <  1:82128/100‑1 (MQ=255)                                       gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgc                                                        >  1:1550539/1‑100 (MQ=255)                                       gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgc                                                        <  2:2008439/100‑1 (MQ=255)                                       gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgc                                                        <  2:978103/100‑1 (MQ=255)                                        aCAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcc                                                       <  2:959736/100‑1 (MQ=255)                                         cAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcct                                                      >  2:2336467/1‑100 (MQ=255)                                          aaTTTTTTTTATCGGTAAATCTCAATGATCAATCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcctg                                                     >  1:231659/1‑100 (MQ=255)                                          aaTTTTTGTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcctg                                                     <  1:2486023/100‑1 (MQ=255)                                          aaTTTTTATTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcctg                                                     >  2:1358983/1‑100 (MQ=255)                                           aTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACctgcctgc                                                    <  1:168954/100‑1 (MQ=255)                                            ttttttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCtgcctgcc                                                   <  2:1101860/100‑1 (MQ=255)                                             tttttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCtgcctgcc                                                   <  1:348862/99‑1 (MQ=255)                                              ttttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGa                                                 <  2:525909/100‑1 (MQ=255)                                                ttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtc                                               <  1:183486/100‑1 (MQ=255)                                                 tttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAATCTGCCTGCCGAtct                                              >  1:1760712/1‑100 (MQ=255)                                                 tttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtct                                              >  1:1726157/1‑100 (MQ=255)                                                  ttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGACTGCCGAtct                                              >  1:949284/1‑99 (MQ=255)                                                    aTCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGc                                           >  2:865262/1‑100 (MQ=255)                                                     tCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGcc                                          <  1:324245/100‑1 (MQ=255)                                                           aaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGCAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGtt                                     <  1:1483094/99‑1 (MQ=255)                                                           aaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGttt                                    <  2:1628144/100‑1 (MQ=255)                                                              ctcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGT‑GAACCTGCCTGCCGATCTCGCCTGGTTTAAAc                                >  1:1084835/1‑100 (MQ=255)                                                                cAATGATCAGTCTGTTTGCGGCTCTGGCGGTAGATCGCGTTATCGGCATGGAAACCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCTTGGTTTAAAc                                <  1:1605666/99‑1 (MQ=255)                                                                cAATGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACg                               >  2:1613/1‑100 (MQ=255)                                                                  aTGATCAGTCTGATTGCGGCGTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                             <  2:2459176/100‑1 (MQ=255)                                                                    gATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                           <  2:889011/100‑1 (MQ=255)                                                                       cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACAcc                        >  2:2327977/1‑100 (MQ=255)                                                                              aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaata                 <  1:2556735/100‑1 (MQ=255)                                                                                tGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataaa               >  1:2006869/1‑100 (MQ=255)                                                                                 gcggcgTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAc              >  2:2188396/1‑100 (MQ=255)                                                                                  cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAcc             >  2:2116089/1‑100 (MQ=255)                                                                                  cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAcc             <  1:763836/100‑1 (MQ=255)                                                                                   ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc            <  2:954260/100‑1 (MQ=255)                                                                                   ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAcc             <  1:842977/99‑1 (MQ=255)                                                                                    gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCg           <  1:2327963/100‑1 (MQ=255)                                                                                      gTTAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg         <  1:1068689/100‑1 (MQ=255)                                                                                       ttAGCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAGCCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGa        <  2:835146/100‑1 (MQ=255)                                                                                        tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt       <  2:709839/100‑1 (MQ=255)                                                                                        tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt       >  2:913992/1‑100 (MQ=255)                                                                                          gCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa     >  1:1877099/1‑100 (MQ=255)                                                                                          gCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa     >  1:817508/1‑100 (MQ=255)                                                                                          gCGGTAGAGCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCAGTGATTa     >  1:1387759/1‑100 (MQ=255)                                                                                            ggTAGAGCGCGTTATCGGCATGGCAAACGCCATGCCGTGGAACCTGCCTGCCGATCCCGCCTGGTTTCAACGCATCACCTTACATAAACCCCTg         >  1:2413634/1‑94 (MQ=255)                                                                                             gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg  >  2:1500327/1‑100 (MQ=255)                                                                                                  |                                                                                               CACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGG  >  NC_000913/49758‑49949

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: ‑