Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,242,141	G→C	13.5%	noncoding (35/76 nt)	PP_t29 →	tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,242,141	0	G	C	13.5%	97.4 / 8.0	37	noncoding (35/76 nt)	PP_t29	tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base G (17/15);  new base C (3/2);  total (20/17)
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.27e-01

ACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCC  >  NC_002947/2241996‑2242287                                                                                                                                                  |                                                                                                                                                   acCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCATAGGGGCCGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggcatg                                                                                                                                                >  1:249419/1‑150 (MQ=255)   cGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATggcaggca                                                                                                                                              <  2:354959/150‑1 (MQ=255)    gCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCag                                                                                                                                             <  1:108633/150‑1 (MQ=255)      ccTTTCACGGCGGTAACAGGGGTTCAAGTCCCCTAGGCGACGCCACTCTTCATACGTGTTGCAGGGCTTAAACAGGATGGTCGATTATTGACGAGCCTAGTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGCCATGCaggc                                                                                                                                           <  2:1132712/150‑2 (MQ=37)        tttCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGCCATGCaggcgc                                                                                                                                         <  1:1201561/150‑4 (MQ=255)                          ggTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGAGGAGCCTTTTGTTGCCGGTCGATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCTGGTCGATGCCGc                                                                                                                       >  1:461696/1‑150 (MQ=37)                           gTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGCTTATTGATGAGCCTTTTGTTGCGGGGCTATAGATCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCt                                                                                                                      >  2:526526/1‑150 (MQ=255)                               gcgTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGc                                                                                                                  >  2:1061698/3‑150 (MQ=255)                                    cccGAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTccac                                                                                                             >  1:1174934/1‑150 (MQ=255)                                        aGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTAAGCGGTTCGATCCCGCTTAGCTCCACCAtt                                                                                                         >  1:16567/1‑150 (MQ=255)                                                      ttCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCa                                                                                           <  2:1174934/150‑1 (MQ=255)                                                          tacGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCa                                                                                       >  1:765711/3‑150 (MQ=255)                                                          tGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCa                                                                                       >  2:844447/1‑150 (MQ=255)                                                          tGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCa                                                                                       >  1:521906/1‑150 (MQ=255)                                                                 tGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGAt                                                                                <  2:249419/150‑1 (MQ=255)                                                                    aGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTc                                                                                                       >  1:800645/1‑124 (MQ=255)                                                                    aGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTc                                                                                                       <  2:800645/124‑1 (MQ=255)                                                                    aGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCtt                                                                             >  1:99171/1‑150 (MQ=255)                                                                                      tcgATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTcccc                                                           >  1:658328/1‑150 (MQ=255)                                                                                      tcgATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTcccc                                                           <  1:912283/150‑1 (MQ=255)                                                                                          ttattGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCg                                                       >  1:138202/1‑150 (MQ=255)                                                                                                atgaGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGtt                                                                   <  1:682555/132‑1 (MQ=255)                                                                                                atgaGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGtt                                                                   >  2:682555/1‑132 (MQ=255)                                                                                                         ttGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCCGCGGTTCGATCCCGCTGAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGa                                        <  1:526526/150‑1 (MQ=255)                                                                                                         ttGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTCGTGGCCTAgca                                        <  2:765711/150‑3 (MQ=255)                                                                                                                  ggCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGAAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCt                               >  1:275041/1‑150 (MQ=255)                                                                                                                      ataGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCa                           <  2:99171/150‑1 (MQ=255)                                                                                                                       taGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTg                                                                                   >  1:692222/1‑93 (MQ=255)                                                                                                                       taGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTg                                                                                   <  2:692222/93‑1 (MQ=255)                                                                                                                        aGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCAcg                         >  2:670333/1‑150 (MQ=255)                                                                                                                         gCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCGTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCAcgg                        <  1:52553/150‑1 (MQ=255)                                                                                                                              gCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTa                   >  1:314587/1‑150 (MQ=255)                                                                                                                                 gggAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACa                <  1:103932/150‑1 (MQ=255)                                                                                                                                        ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTGCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTc         >  1:86879/3‑150 (MQ=255)                                                                                                                                         tccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCg        >  2:1302725/2‑150 (MQ=255)                                                                                                                                          ccTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa       <  2:314587/150‑1 (MQ=255)                                                                                                                                          ccTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa       <  2:521906/150‑1 (MQ=255)                                                                                                                                           cTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTAAGTCCCCTTCGTCTAGCGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTCCGAg      >  1:1339093/1‑150 (MQ=255)                                                                                                                                            tGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGt     <  2:461696/150‑1 (MQ=255)                                                                                                                                               atggcatgCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTccc  <  1:844447/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   ACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCC  >  NC_002947/2241996‑2242287

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