Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	344,141	T→A	9.7%	noncoding (35/76 nt)	PP_t72 ←	tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	344,141	0	T	A	9.7%	136.1 / 3.7	49	noncoding (35/76 nt)	PP_t72	tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base A (3/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGATTTTTGGGGTCATTCGGCCGGTGGAACG  >  NC_002947/343994‑344284                                                                                                                                                     |                                                                                                                                                cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                                                                                                                                 <  6:14505/148‑1 (MQ=255) cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                                                                                                                                 <  7:8669/148‑1 (MQ=255)    cGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCaa                                                                                                                                              <  1:9041/148‑1 (MQ=255)    cGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCa                                                                                                                                               <  6:71907/147‑1 (MQ=255)          cgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTGAAAATCctcgt                                                                                                                                       <  8:62154/149‑3 (MQ=255)          cgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTGAAAATCctc                                                                                                                                         <  6:184467/147‑1 (MQ=255)          cgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTGAAAATCctc                                                                                                                                         >  5:184467/1‑147 (MQ=255)               gAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTa                                                                                                                                  >  5:131050/1‑149 (MQ=255)                aaGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAc                                                                                                                                 <  2:238369/149‑1 (MQ=255)                aaGTTAAAAATCACAGACAACAAAAAACACTCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTa                                                                                                                                  <  7:187877/148‑1 (MQ=255)                     aaaaaTCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACt                                                                                                                            >  1:250893/1‑149 (MQ=255)                                 caacaaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGTTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTggg                                                                                                                <  2:250893/149‑1 (MQ=255)                                 caacaaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTggg                                                                                                                <  2:180528/149‑1 (MQ=255)                                  aacaaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc                                                                                                               >  8:220207/1‑149 (MQ=255)                                    caaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTg                                                                                                                           <  1:231712/135‑1 (MQ=255)                                    caaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTg                                                                                                                           >  2:231712/1‑135 (MQ=255)                                     aagaaaCCCGC‑CTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACGAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac                                                                                                            >  2:209311/4‑149 (MQ=255)                                     aagaaaCCCGC‑CTAGTGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACAACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTACAATCCTCTGATCTACCGACTGAGCTATCTGGGcgac                                                                                                            >  6:76130/4‑149 (MQ=255)                                      aaaaaCCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgacg                                                                                                           <  2:36335/149‑1 (MQ=255)                                            ccGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCTCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTg                                                                                                      >  5:12312/1‑148 (MQ=255)                                            ccGCACTAG‑GCGCGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTg                                                                                                      >  7:208431/1‑148 (MQ=255)                                                  tAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTAt                                                                                                                     <  7:4115/127‑1 (MQ=255)                                                  tAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTAt                                                                                                                     >  8:4115/1‑127 (MQ=255)                                                           tttCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTc                                                                                       <  4:189449/149‑1 (MQ=255)                                                                tgttgtCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATACTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGAc                                                                                  >  7:189891/1‑149 (MQ=255)                                                                    gtCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                                                                                                                                 <  7:21038/82‑1 (MQ=255)                                                                    gtCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                                                                                                                                 >  8:21038/1‑82 (MQ=255)                                                                        cTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAATTAGATTTCCGGACCGCTCGTc                                                                          <  6:131050/149‑1 (MQ=255)                                                                            gggTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAAc                                                                       <  6:12312/148‑1 (MQ=255)                                                                             ggTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGa                                                                     <  7:220207/149‑1 (MQ=255)                                                                                    ttGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTATACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACttttt                                                               >  4:162409/1‑148 (MQ=255)                                                                                          ccaccaGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTTGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGaaaaaa                                                        >  8:102693/1‑149 (MQ=255)                                                                                           caccaGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAAt                                                       >  6:190269/1‑149 (MQ=255)                                                                                             ccaGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACttttt                                                               >  6:126246/1‑139 (MQ=255)                                                                                             ccaGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACttttt                                                               <  5:126246/139‑1 (MQ=255)                                                                                                gAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTa                                                   >  7:144189/1‑148 (MQ=255)                                                                                                 aaGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTTGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTaaa                                                 <  5:76130/149‑1 (MQ=255)                                                                                                         aGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTcc                                         <  8:215216/149‑1 (MQ=255)                                                                                                                 cccAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTa                                 <  8:144189/149‑1 (MQ=255)                                                                                                                         cGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGAttttt                         <  8:208431/149‑1 (MQ=255)                                                                                                                             gTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTACGTCGCTTACGATTTTTgggg                     >  7:1096/1‑149 (MQ=255)                                                                                                                                gAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCt                                   >  1:175538/1‑132 (MQ=255)                                                                                                                                gAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCt                                   <  2:175538/132‑1 (MQ=255)                                                                                                                                gAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGATTTTTGGGGTCa                  >  1:194845/1‑149 (MQ=255)                                                                                                                                           acGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGATTTTTGGGGt                    >  5:53710/1‑136 (MQ=255)                                                                                                                                           acGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGATTTTTGGGGt                    <  6:53710/136‑1 (MQ=255)                                                                                                                                            agAGGATTGAAAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCAGTCGCTTACGATTTTTTGGGTCATTCGGCCGGTg       >  1:197936/2‑148 (MQ=255)                                                                                                                                             gAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGATTTTTGGGGTCATTCGGCCGGTgg      >  3:201654/1‑148 (MQ=255)                                                                                                                                              aGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGAGCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGATTTTTGGGGTCATTCGGCCGGTGGaa    <  5:190269/149‑1 (MQ=255)                                                                                                                                               ggATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAAt                                                       >  5:168808/1‑97 (MQ=255)                                                                                                                                               ggATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAAt                                                       <  6:168808/97‑1 (MQ=255)                                                                                                                                                gATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGATTTTTGGGGTCATTCGGCCGGTGGAACg  >  3:117637/1‑149 (MQ=255)                                                                                                                                                gATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGATTTTTGGGGTCATTCGGCCGGTGGAACg  >  1:246247/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                CTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAATAGATTTCCGGACCGCTCGTCAACGACTTTTTGAAAAAATTTTAAATTAATTCCGTCGCTTACGATTTTTGGGGTCATTCGGCCGGTGGAACG  >  NC_002947/343994‑344284

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