Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,055,250	A→G	17.0%	S64G (AGC→GGC)	amn →	AMP nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,055,250	0	A	G	17.0%	106.9 / 17.4	53	S64G (AGC→GGC)	amn	AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base G (2/7);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGCCCTACTCTCTTTCGTTC  >  NC_000913/2055153‑2055344                                                                                                  |                                                                                               cATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGc                                                                                              <  2:2183488/100‑1 (MQ=255)  aTTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTGGca                                                                                             <  1:1810775/100‑1 (MQ=255)  aTTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGca                                                                                             <  1:988579/100‑1 (MQ=255)      gCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAAc                                                                                         <  1:398026/100‑1 (MQ=255)       cAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAAcc                                                                                        <  2:278800/100‑1 (MQ=255)             tatTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTCGGACGGTAGCACAACCAATccc                                                                                  <  2:1386827/100‑1 (MQ=255)                tACAAGTAGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATcccccc                                                                               <  2:361187/100‑1 (MQ=255)                               aCCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCTCCCAAAACGCGCGCAttt                                                                >  1:2247202/1‑100 (MQ=255)                                ccTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTg                                                               >  1:351493/1‑100 (MQ=255)                                ccTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTg                                                               >  2:1497657/1‑100 (MQ=255)                                   gatgaAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTc                                                            <  2:1564418/100‑1 (MQ=255)                                   gatgaAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTc                                                            >  2:347901/1‑100 (MQ=255)                                    atgaAAACGCCCGCAAACAAGGTCTTTTTGTCTATCTATCACTTTCTTTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCg                                                           >  2:72753/1‑100 (MQ=255)                                     tgaAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGc                                                          <  1:107197/100‑1 (MQ=255)                                       aaaaCGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTGGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCtt                                                        >  2:1128778/1‑100 (MQ=255)                                        aaaCGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTGGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCttt                                                       <  1:2310259/100‑1 (MQ=255)                                         aaCGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTa                                                      <  1:1667809/100‑1 (MQ=255)                                           cGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTAcc                                                    <  2:1090817/100‑1 (MQ=255)                                           cGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTAcc                                                    <  2:23635/100‑1 (MQ=255)                                             cccGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTGGCACAACCAATCCCCCCAAACCGCGCGCATTTGGTCGCTTTACCCa                                                  >  1:1781959/1‑100 (MQ=255)                                             cccGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCa                                                  >  1:1519236/1‑100 (MQ=255)                                              ccGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCAc                                                 >  2:1565211/1‑100 (MQ=255)                                               cGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACg                                                >  2:904033/1‑100 (MQ=255)                                               cGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACg                                                <  1:9888/100‑1 (MQ=255)                                                 cAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACgca                                              <  2:222240/100‑1 (MQ=255)                                                 cAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACgca                                              <  2:2268074/100‑1 (MQ=255)                                                   aaCAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACgcagg                                            <  2:1111891/100‑1 (MQ=255)                                                    aCAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTGGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACgcaggc                                           <  1:2216588/100‑1 (MQ=255)                                                    aCAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTGGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACgcaggc                                           <  1:1366609/100‑1 (MQ=255)                                                    aCAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCACTCCCCCCAACACGCTCGCATTTGGTCGCTTTACCCACgcaggc                                           <  2:1642111/100‑1 (MQ=255)                                                    aCAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACgcaggc                                           <  1:1974867/100‑1 (MQ=255)                                                        ggTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCt                                       >  2:1667527/1‑100 (MQ=255)                                                          tCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTac                                     <  2:1240402/100‑1 (MQ=255)                                                                tGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTAcaccacc                               >  2:600192/1‑100 (MQ=255)                                                                tGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTAcaccacc                               >  1:279937/1‑100 (MQ=255)                                                                  tCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCACTCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTAcaccaccac                             <  2:1810469/100‑1 (MQ=255)                                                                   ctatctatCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGAATTTGGTCGCTTTACCCACGCAGGCAGCTAcaccaccac                             >  2:251713/1‑99 (MQ=255)                                                                     aTCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGAt                          >  1:1359818/1‑100 (MQ=255)                                                                      tCCATCACTGACCGTAACCTGGGACGGTGGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGAtt                         <  1:1399014/100‑1 (MQ=255)                                                                       ccATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTa                        >  2:1909036/1‑100 (MQ=255)                                                                       ccATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTa                        >  2:1709466/1‑100 (MQ=255)                                                                        cATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTAc                       >  1:1233817/1‑100 (MQ=255)                                                                          tCACTGACCGTAACCTGGGACGGTGGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTc                     <  1:1525017/100‑1 (MQ=255)                                                                              tGACCGTAACCTGGGACGGTGGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGccc                 <  2:2139159/100‑1 (MQ=255)                                                                              tGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGccc                 >  2:623518/1‑100 (MQ=255)                                                                               gACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGCCCt                <  2:530013/100‑1 (MQ=255)                                                                                 ccGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGCCCTAc              >  2:2281662/1‑100 (MQ=255)                                                                                  cGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGCCCTAct             <  1:261640/100‑1 (MQ=255)                                                                                  cGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGCCCTAct             >  1:1823477/1‑100 (MQ=255)                                                                                        cTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGCCCTACTCTCttt       >  1:1445061/1‑100 (MQ=255)                                                                                          gggACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGCCCTACTCTCTttcg     <  1:1127549/100‑1 (MQ=255)                                                                                            gACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGCCCTACTCTCTttcgtt   <  2:1896705/100‑1 (MQ=255)                                                                                              cGGTAGCACAAACTATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCTGGCAGCTACACCACCACGATTACTCGCTCTACCCTCTttcgttc  <  2:826502/99‑1 (MQ=255)                                                                                                  |                                                                                               CATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCCCCAAAACGCGCGCATTTGGTCGCTTTACCCACGCAGGCAGCTACACCACCACGATTACTCGCCCTACTCTCTTTCGTTC  >  NC_000913/2055153‑2055344

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