Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,048,459	C→T	A498A (GCC→GCT)	polA →	DNA polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,048,459	0	C	T	97.7%	148.4 / ‑4.2	43	A498A (GCC→GCT)	polA	DNA polymerase I
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/16);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TTGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCCGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGTTG  >  NC_000913/4048383‑4048534                                                                             |                                                                            ttGAAGAGATTGCTGGTAAAGGCAAAAATCGACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACgctgcc                                                                          <  2:747534/80‑1 (MQ=255)  tGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGctgccg                                                                         <  1:747537/80‑1 (MQ=255)  tGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGctgccg                                                                         <  1:747536/80‑1 (MQ=255)      gagaTTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGAt                                                                    >  2:747602/1‑81 (MQ=255)      gagaTTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGAt                                                                    >  1:747584/1‑81 (MQ=255)       agaTTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATg                                                                   <  1:747538/81‑1 (MQ=255)         aTTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATg                                                                   >  1:747607/1‑79 (MQ=255)              tGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGt                                                            >  2:747584/1‑81 (MQ=255)               ggTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGt                                                            >  2:747605/1‑80 (MQ=255)               ggTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGt                                                            >  2:747598/1‑80 (MQ=255)                gTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTc                                                           <  1:747539/80‑1 (MQ=255)                gTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTc                                                           <  2:747540/80‑1 (MQ=255)                gTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCa                                                          <  2:747539/81‑1 (MQ=255)                      gCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGc                                                    <  1:747542/81‑1 (MQ=255)                            aTCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTg                                               >  1:747602/1‑80 (MQ=255)                            aTCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTg                                               >  2:747586/1‑80 (MQ=255)                            aTCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGc                                              <  2:747542/81‑1 (MQ=255)                             tCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCa                                             >  2:747591/1‑81 (MQ=255)                               aaCTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATc                                           >  2:747593/1‑81 (MQ=255)                               aaCTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATc                                           <  2:747544/81‑1 (MQ=255)                               aaCTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATc                                           <  2:747543/81‑1 (MQ=255)                                aCTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCt                                          >  2:747581/1‑81 (MQ=255)                                aCTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCt                                          <  2:747545/81‑1 (MQ=255)                                   gACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGa                                        >  1:747593/1‑80 (MQ=255)                                      cTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAAt                                    >  2:747600/1‑81 (MQ=255)                                          aaCCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTgg                                >  1:747609/1‑81 (MQ=255)                                             cAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCg                             >  1:747596/1‑81 (MQ=255)                                               gATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGa                           <  2:747546/81‑1 (MQ=255)                                                 ttGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATc                         <  1:747548/81‑1 (MQ=255)                                                   gCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCt                        <  1:747549/80‑1 (MQ=255)                                                   gCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTg                       >  2:747611/1‑81 (MQ=255)                                                   gCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTg                       >  2:747599/1‑81 (MQ=255)                                                     ccTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCa                     <  1:747550/81‑1 (MQ=255)                                                      cTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCaa                    <  2:747550/81‑1 (MQ=255)                                                         gaagaaGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCaaaaa                 >  1:747586/1‑81 (MQ=255)                                                          aagaagCCGGACGTTACGCCGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCaaaaa                 >  1:747591/1‑80 (MQ=255)                                                            gaagCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAaca               >  1:747604/1‑80 (MQ=255)                                                             aagCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAacac              >  2:747615/1‑80 (MQ=255)                                                              agCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAAcaca             >  1:747589/1‑80 (MQ=255)                                                                ccGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAg          <  2:747551/81‑1 (MQ=255)                                                                 cGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAgg         >  2:747621/1‑81 (MQ=255)                                                                  ggACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAggg        >  1:747621/1‑81 (MQ=255)                                                                   gACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGc       >  1:747617/1‑81 (MQ=255)                                                                        tACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGtt   >  1:747620/1‑80 (MQ=255)                                                                         aCGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGTTg  >  1:747619/1‑80 (MQ=255)                                                                         aCGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGTTg  <  1:747553/80‑1 (MQ=255)                                                                             |                                                                            TTGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCCGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGTTG  >  NC_000913/4048383‑4048534

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