Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,898	G→A	A26T (GCC→ACC)	folA →	dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,898	0	G	A	92.6%	79.4 / ‑1.6	27	A26T (GCC→ACC)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (14/11);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

TGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCA  >  NC_000913/49837‑49956                                                              |                                                           tGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCg                                                         >  1:743430/1‑65 (MQ=255)  gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCCTGCCGTGGAACCTGCCTACCGa                                                        <  1:1232463/65‑1 (MQ=255)  gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGa                                                        <  1:2409240/65‑1 (MQ=255)     gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAACCGCCATGCCGTGGAACCTGCCTACCGAtct                                                     <  1:2041261/65‑1 (MQ=255)       gTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCg                                                   <  1:348309/65‑1 (MQ=255)              gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGtt                                            >  1:269506/1‑65 (MQ=255)                 gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTaa                                         <  1:12390/65‑1 (MQ=255)                  aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGGTTcaa                                        >  1:2450573/1‑65 (MQ=255)                  aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCCACCGCTCTCGCCAGGGTTaaa                                        >  1:2807118/1‑65 (MQ=255)                  aTCGCGTTATCGGCATGGAAAAAGCCATGCCGGGGAACCCGCCTACCGATCTCGCCAGGTTTaaa                                        >  1:1231428/1‑65 (MQ=255)                   tCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAAc                                       >  1:2312358/1‑65 (MQ=255)                     gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGc                                     >  1:449771/1‑65 (MQ=255)                         tATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAaca                                 >  1:3005932/1‑65 (MQ=255)                          aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAacac                                <  1:369262/65‑1 (MQ=255)                          aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAacac                                <  1:300934/65‑1 (MQ=255)                          aTCGGCATGGAAAACGCCATGCCGGGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                <  1:554746/65‑1 (MQ=255)                           tccgCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACAcc                               <  1:217777/62‑1 (MQ=255)                           tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACAcc                               <  1:1684096/65‑1 (MQ=255)                           tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACAcc                               <  1:2766016/65‑1 (MQ=255)                              gCATGGAAAACGCCATGCCGTGGACCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTta                            <  1:1617691/65‑1 (MQ=255)                               cATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTtaa                           <  1:698083/65‑1 (MQ=255)                                    aaaaCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTtaaataaa                      >  1:297565/1‑65 (MQ=255)                                        cGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCg                  >  1:1050943/1‑65 (MQ=255)                                         gCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGt                 >  1:851175/1‑65 (MQ=255)                                         gCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGt                 >  1:1522890/1‑65 (MQ=255)                                               ccGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGCTTAt           >  1:1426974/1‑65 (MQ=255)                                                     aaCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg     >  1:2275479/1‑65 (MQ=255)                                                      aCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc    >  1:2452010/1‑65 (MQ=255)                                                      aCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc    >  1:695492/1‑65 (MQ=255)                                                      aCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCCTAAATAAACCCGTGATTATGGgccgc    >  1:1920070/1‑65 (MQ=255)                                                        ctgcctacCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcta  >  1:1124728/1‑63 (MQ=255)                                                              |                                                           TGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCA  >  NC_000913/49837‑49956

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑