Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,241,611	G→C	8.8%	noncoding (35/76 nt)	PP_t25 →	tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,241,611	0	G	C	8.8%	125.2 / 2.8	44	noncoding (35/76 nt)	PP_t25	tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base C (2/2);  total (23/21)
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.76e-01

TTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACAC  >  NC_002947/2241472‑2241752                                                                                                                                            |                                                                                                                                              tttGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCaggag                                                                                                                                     >  1:835817/1‑150 (MQ=255)    gACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTc                                                                                                                                  >  1:992948/1‑150 (MQ=255)       aGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGc                                                                                                                               >  1:1061120/1‑150 (MQ=255)         cccTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCgg                                                                                                                             >  2:913764/1‑150 (MQ=255)                    gAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGc                                                                                                                  >  1:524778/1‑150 (MQ=255)                       cTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTa                                                                                                               <  2:992948/150‑1 (MQ=255)                        ttAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAg                                                                                                              >  2:70399/1‑150 (MQ=255)                             atGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTcca                                                                                                         <  2:1061120/150‑1 (MQ=255)                                        ccACAGCGATGCACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGc                                                                                              <  1:848914/150‑1 (MQ=255)                                           caGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTa                                                                                           <  2:890419/150‑1 (MQ=255)                                               gATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagcag                                                                                       >  2:129172/1‑150 (MQ=255)                                               gATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGAATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagcag                                                                                       >  2:396232/1‑150 (MQ=255)                                                  gAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTCGCAGCTa                                                                                    >  1:50582/1‑150 (MQ=255)                                                     cgaaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCag                                                                                                                                        <  1:1019612/95‑1 (MQ=255)                                                     cgaaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCag                                                                                                                                        >  2:1019612/1‑95 (MQ=255)                                                      gaaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGc                                                                                <  1:129172/150‑1 (MQ=255)                                                      gaaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGc                                                                                <  2:835817/150‑1 (MQ=255)                                                       aaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCt                                                                               <  2:1121619/150‑1 (MQ=255)                                                            gaaacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATaa                                                                          >  2:929049/1‑150 (MQ=255)                                                            gaaacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATaa                                                                          <  2:524778/150‑1 (MQ=255)                                                                cgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCaa                                                                      <  1:913764/150‑1 (MQ=255)                                                                  ccAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGa                                                                    >  2:607734/1‑150 (MQ=255)                                                                              ccagcGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCt                                                                               >  1:544654/1‑127 (MQ=255)                                                                              ccagcGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCt                                                                               <  2:544654/127‑1 (MQ=255)                                                                              ccagcGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGAATGCAGGAGGTCAGCGGTTCGATCCCGATTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGCCCGGGTCCAgac                                                        >  2:208411/1‑148 (MQ=255)                                                                                 gcGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAAcc                                                     >  1:214738/1‑150 (MQ=255)                                                                                   aaTTCAGTCCACAGCCACATTGTGGGGCTATAGCTCACGTGGGAGAACGCCTGCATGGCATGCAGGAGGTCAGCGTTTCGATCCCGCTTAGCACCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCgg                                                   <  2:214738/149‑1 (MQ=25)                                                                                      tCAGTTCAAAGAGACATAGGCCGTCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTc                                                <  1:208411/150‑1 (MQ=25)                                                                                              aaGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCCTGATCTGCAGGAGGTGTGCGCGTCGATCCAGCTTCGGTCCACCAAATTCCGCATGGCTGCTAATGCGGATGAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGt                                        <  1:929049/150‑1 (MQ=21)                                                                                               aGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTa                                       >  2:862449/1‑150 (MQ=255)                                                                                                        gtgGGGCTATAGCTCAGCTGGGAGAGCGTCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTc                              <  1:743957/150‑1 (MQ=255)                                                                                                           gggCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGtt                               <  1:958993/146‑1 (MQ=255)                                                                                                           gggCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGtt                               >  2:958993/1‑146 (MQ=255)                                                                                                           gggCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTc                           <  2:1151594/150‑1 (MQ=255)                                                                                                            ggCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTcc                          >  1:961901/1‑150 (MQ=255)                                                                                                                    cTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTaa                                      <  2:120698/130‑1 (MQ=255)                                                                                                                    cTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTaa                                      >  1:120698/1‑130 (MQ=255)                                                                                                                     tCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCt                 >  1:770734/1‑150 (MQ=255)                                                                                                                        gCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGt              <  2:961901/150‑1 (MQ=255)                                                                                                                                 cctccTTCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGt              >  2:924071/4‑141 (MQ=255)                                                                                                                                 cctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGt              <  1:924071/138‑1 (MQ=255)                                                                                                                                  ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCt                                                                                             >  2:654858/3‑61 (MQ=255)                                                                                                                                  ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCt                                                                                             <  1:654858/59‑1 (MQ=255)                                                                                                                                    ccTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGacac  >  2:397231/1‑150 (MQ=255)                                                                                                                                            |                                                                                                                                              TTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACAC  >  NC_002947/2241472‑2241752

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