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.8%	160.3 / ‑4.2	46	A498A (GCC→GCT)	polA	DNA polymerase I
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (23/22);  total (23/23)
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.85e-01

TTGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCCGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGTTGA  >  NC_000913/4048383‑4048535                                                                             |                                                                             ttGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACgctgcc                                                                           <  2:815025/80‑1 (MQ=255) ttGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACgctgcc                                                                           <  2:815024/80‑1 (MQ=255)  tGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGctgccg                                                                          <  1:815026/80‑1 (MQ=255)  tGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGctgccg                                                                          >  1:815079/1‑80 (MQ=255)   gAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGctgccg                                                                          >  2:815083/1‑79 (MQ=255)   gAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGctgccg                                                                          >  1:815103/1‑79 (MQ=255)   gAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGa                                                                         <  2:815027/80‑1 (MQ=255)    aaGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGaa                                                                        <  1:815028/80‑1 (MQ=255)        gaTTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATg                                                                    >  1:815099/1‑80 (MQ=255)        gaTTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGc                                                                   <  2:815029/81‑1 (MQ=255)             cTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATg                                                              <  2:815030/81‑1 (MQ=255)                gTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGt                                                             <  1:815032/79‑1 (MQ=255)                gTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTc                                                            <  2:815031/80‑1 (MQ=255)                 tAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCAc                                                          <  1:815033/81‑1 (MQ=255)                  aaaGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCAcc                                                         >  2:815086/1‑81 (MQ=255)                   aaGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCt                                                        >  1:815089/1‑81 (MQ=255)                   aaGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCt                                                        >  2:815093/1‑81 (MQ=255)                   aaGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCt                                                        >  1:815105/1‑81 (MQ=255)                        aaaaaTCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAg                                                   <  2:815034/81‑1 (MQ=255)                         aaaaTCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGt                                                  >  2:815082/1‑81 (MQ=255)                          aaaTCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGtt                                                 >  2:815084/1‑81 (MQ=255)                          aaaTCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGtt                                                 >  2:815087/1‑81 (MQ=255)                           aaTCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTg                                                >  1:815118/1‑81 (MQ=255)                             tCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCa                                              <  2:815035/81‑1 (MQ=255)                               aaCTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATc                                            <  1:815037/81‑1 (MQ=255)                               aaCTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATc                                            <  1:815036/81‑1 (MQ=255)                                 cTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCCGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTg                                          <  2:815038/81‑1 (MQ=255)                                   gACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGa                                         >  2:815072/1‑80 (MQ=255)                                    aCCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGAAGTCACCTTGCAGTTGCATCTGaa                                        <  2:815039/80‑1 (MQ=255)                                      cTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAAt                                     <  1:815040/81‑1 (MQ=255)                                         tAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAAtgtg                                  >  2:815098/1‑81 (MQ=255)                                          aaCCAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAAtgtg                                  <  1:815041/80‑1 (MQ=255)                                            ccAGATTGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGcc                               <  1:815042/81‑1 (MQ=255)                                                 ttGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATc                          >  2:815114/1‑81 (MQ=255)                                                  tGCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCt                         >  2:815110/1‑81 (MQ=255)                                                   gCCCTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTg                        >  1:815106/1‑81 (MQ=255)                                                    cccTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGc                       <  2:815043/81‑1 (MQ=255)                                                      cTCGAAGAAGCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCa                      >  2:815117/1‑80 (MQ=255)                                                          aagaagCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCaaaaa                  >  2:815094/1‑80 (MQ=255)                                                          aagaagCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAac                 <  2:815044/81‑1 (MQ=255)                                                            gaagCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAacac               <  1:815045/81‑1 (MQ=255)                                                              agCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAAcaca              >  1:815123/1‑80 (MQ=255)                                                              agCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAAcaca              >  1:815091/1‑80 (MQ=255)                                                               gCCGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAAcaca              >  1:815100/1‑79 (MQ=255)                                                                ccGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAg           <  1:815046/81‑1 (MQ=255)                                                                 cGGACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAgg          >  2:815107/1‑81 (MQ=255)                                                                  ggACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAggg         >  1:815085/1‑81 (MQ=255)                                                                   gACGTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGG‑CAGATCTGCAAAAACACAAAGGGcc       >  2:815146/1‑81 (MQ=255)                                                                      gTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGt     <  2:815048/81‑1 (MQ=255)                                                                      gTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGt     >  1:815095/1‑81 (MQ=255)                                                                      gTTACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGt     <  1:815047/81‑1 (MQ=255)                                                                       ttACGCTGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGtt    <  2:815049/81‑1 (MQ=255)                                                                          cgctgcCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGTTGa  <  2:815050/80‑1 (MQ=255)                                                                             |                                                                             TTGAAGAGATTGCTGGTAAAGGCAAAAATCAACTGACCTTTAACCAGATTGCCCTCGAAGAAGCCGGACGTTACGCCGCCGAAGATGCAGATGTCACCTTGCAGTTGCATCTGAAAATGTGGCCGGATCTGCAAAAACACAAAGGGCCGTTGA  >  NC_000913/4048383‑4048535

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