Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,396,689	G→A	S369S (AGC→AGT)	rng ←	ribonuclease G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,396,689	0	G	A	95.8%	223.2 / 0.8	71	S369S (AGC→AGT)	rng	ribonuclease G
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (30/38);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 9.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

CAGTACGTGCTCAATGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTGCTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCCCGCCCAGATTACGC  >  NC_000913/3396593‑3396785                                                                                                 |                                                                                                 caGTACGTGCTCAATGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTc                                                                                               <  2:894506/100‑1 (MQ=255)    tACGTGCTCAATGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAAc                                                                                            <  1:905830/100‑1 (MQ=255)      cGTGCTCAATGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGc                                                                                          >  1:1042359/1‑100 (MQ=255)        tGCTCAATGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCt                                                                                        >  1:1562154/1‑100 (MQ=255)        tGCTCAATGCTTTCGCGGGTGAGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTGAACg                                                                                           >  1:662689/1‑97 (MQ=255)         gCTCAATGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTg                                                                                       >  1:668258/1‑100 (MQ=255)            cAATGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTc                                                                                    >  1:105256/1‑100 (MQ=255)            cAATGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTc                                                                                    <  2:854268/100‑1 (MQ=255)             aaTGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCAGTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTcc                                                                                   <  1:994764/100‑1 (MQ=255)              aTGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCa                                                                                  >  2:274246/1‑100 (MQ=255)              aTGCTGTCGCGGGTGCGTTTACGCGTCATCTCCACCAAAAACAGAAAAAAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTcc                                                                                   <  1:490652/99‑1 (MQ=255)               tGCTTTCGCGGGTGCGTTTACTCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAg                                                                                 <  2:410383/100‑1 (MQ=255)               tGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAg                                                                                 <  2:1362024/100‑1 (MQ=255)               tGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCAACCGGTCTTTACTCAACGCCTGCTCCAg                                                                                 >  2:69658/1‑100 (MQ=255)                gCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAgn                                                                                <  2:586328/100‑2 (MQ=255)                 cTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCg                                                                               <  2:67885/100‑1 (MQ=255)                 cTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCg                                                                               >  1:1249357/1‑100 (MQ=255)                     cgcgGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTg                                                                           >  2:1364834/1‑100 (MQ=255)                      gcgGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTgc                                                                          <  1:1814762/100‑1 (MQ=255)                         ggTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTgcagc                                                                       <  1:1034859/100‑1 (MQ=255)                         ggTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTgcagc                                                                       >  1:337572/1‑100 (MQ=255)                         ggTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTgcagc                                                                       <  1:1308910/100‑1 (MQ=255)                            gCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACt                                                                    <  2:1010435/100‑1 (MQ=255)                              gTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTcg                                                                  >  2:4547/1‑100 (MQ=255)                              gTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTcg                                                                  >  1:1470358/1‑100 (MQ=255)                                   cgcgTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTGCTCAACGCCTGCTCCAGCGAGTGCAGCACTcggcggc                                                             >  2:1672131/1‑100 (MQ=255)                                   cgcgTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGATTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTcggcggc                                                             >  1:1005869/1‑100 (MQ=255)                                    gcgTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCTGCACTcgctggag                                                            <  2:1372971/100‑7 (MQ=255)                                     cgTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTcggcggcgg                                                           <  1:854441/100‑1 (MQ=255)                                         aTCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGa                                                        >  2:347890/1‑99 (MQ=255)                                         aTCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGCCGGTGa                                                        >  1:219612/1‑99 (MQ=255)                                          tctcCTCCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATc                                                      <  1:1359414/100‑1 (MQ=255)                                          tctcCACCAGCCCCAGCGCCGAAAAACCATTATCGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATc                                                      <  2:115613/100‑1 (MQ=255)                                           ctcCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCt                                                     <  2:865486/100‑1 (MQ=255)                                            tcCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTTCAGCACTCGGCAGCGGTAATCtt                                                    >  2:414876/1‑100 (MQ=255)                                             ccaccaGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCt                                                     >  1:659503/1‑98 (MQ=255)                                              caccaGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCa                                                  <  2:734037/100‑1 (MQ=255)                                              caccaGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCa                                                  >  2:658311/1‑100 (MQ=255)                                               accaGCCCCAGCGCCGAAAAACCATTATCGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCat                                                 <  2:599984/100‑1 (MQ=255)                                               accaGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCat                                                 <  2:104399/100‑1 (MQ=255)                                               accaGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCat                                                 >  1:410180/1‑100 (MQ=255)                                                 caGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCat                                                 <  2:1469180/98‑1 (MQ=255)                                                 caGCCCCAGCGCCGAAAAACCATTAACGCTCGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCatta                                               <  1:359496/100‑1 (MQ=255)                                                   gCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCattatt                                             <  1:237859/100‑1 (MQ=255)                                                   gCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCattatt                                             >  2:386787/1‑100 (MQ=255)                                                    ccccAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGTTCTTTACTCATCGCCTGCTCCATCGAGTGCAGCACTCGGCGGCAGTGATCTTCATTATTc                                            <  1:1640902/100‑1 (MQ=255)                                                     cccAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCa                                           <  1:200593/100‑1 (MQ=255)                                                          cgcCTAAAAACCATTAACGCTTGTTTTCACGCGTTCTTCACTCAACACCTGCTCCAGCGAGTGCATCACTCGGCTCCGGTGATCTACATTATTCAGATCg                                      <  1:733444/100‑1 (MQ=255)                                                          cgcCTAAAAACCATTAACGCTTGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCTACTGCATCACTCGGCGGCGGTTATCTTCATTATTCata                                         <  1:367746/97‑1 (MQ=255)                                                            ccGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTACAGAACACGACGGCGGTGATCTTCATTATTCATATCGAt                                    >  2:1510572/1‑100 (MQ=255)                                                             cGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTGCTCAACGCCTGCTCCAGCGAGTGCAGCACTCGACGGCGGTGATCTTCATTATTCATATCGATg                                   >  2:988070/1‑100 (MQ=255)                                                             cGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATg                                   <  2:337574/100‑1 (MQ=255)                                                              gAAAAACCATTAACGCTGGATTTCAACCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGAACTTCATTATTCATATCGATGa                                  >  1:276993/1‑100 (MQ=255)                                                                  aaCCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATc                              <  2:1470369/100‑1 (MQ=255)                                                                  aaCCATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATc                              <  2:1735371/100‑1 (MQ=255)                                                                     cATTAACGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCaat                           >  2:741186/1‑100 (MQ=255)                                                                         aaCGCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATg                       >  1:710486/1‑100 (MQ=255)                                                                            gCTGGTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATa                    <  2:982408/100‑1 (MQ=255)                                                                               ggTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATc                 >  2:673661/1‑100 (MQ=255)                                                                               ggTTTTCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATATTCATTATTCATATCGATGAAATCAATAATGATAAt                  <  1:338998/99‑1 (MQ=255)                                                                                  tttCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATccc               >  2:899610/1‑99 (MQ=255)                                                                                  tttCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCCCg              >  1:331157/1‑100 (MQ=255)                                                                                   ttCACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCCCGc             <  1:1759712/100‑1 (MQ=255)                                                                                     cACCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCCCGccc           <  2:947102/100‑1 (MQ=255)                                                                                      aCCCGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCCAGCCCa          <  2:276994/100‑1 (MQ=255)                                                                                       cccGGTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCCCGCCCAg         >  2:461193/1‑100 (MQ=255)                                                                                        ccGGTCTTTGCTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCACGcaaaga        >  2:786338/1‑95 (MQ=255)                                                                                          ggTCTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCTTATCGATGAAATCAATAATGATAATCCCGCCCAGAtt      >  1:188912/1‑100 (MQ=255)                                                                                             cTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCCCGCCCAGATTACg   <  2:314687/100‑1 (MQ=255)                                                                                             cTTTACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCCCGCCCAGATTACg   >  2:1009396/1‑100 (MQ=255)                                                                                              tttACTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAGTAATGATAATCCCGCCCAGATTACGc  <  1:1642891/100‑1 (MQ=255)                                                                                                 |                                                                                                 CAGTACGTGCTCAATGCTTTCGCGGGTGCGTTTACGCGTCATCTCCACCAGCCCCAGCGCCGAAAAACCATTAACGCTGGTTTTCACCCGGTCTTTGCTCAACGCCTGCTCCAGCGAGTGCAGCACTCGGCGGCGGTGATCTTCATTATTCATATCGATGAAATCAATAATGATAATCCCGCCCAGATTACGC  >  NC_000913/3396593‑3396785

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑