Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,380,723	T→A	7.7%	intergenic (‑111/+24)	PP_t67 ← / ← PP_t66	tRNA‑Met/tRNA‑Met

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,380,723	0	T	A	7.7%	149.4 / 3.0	52	intergenic (‑111/+24)	PP_t67/PP_t66	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base A (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTC  >  NC_002947/1380580‑1380861                                                                                                                                                |                                                                                                                                            gCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACa                                                                                                                                        <  3:22091/149‑1 (MQ=255)    cgacgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTg                                                                                                                                     <  5:353340/149‑1 (MQ=255)      acgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTACGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGc                                                                                                                                   <  7:198417/149‑1 (MQ=18)         aGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑tt                                                                                                                               <  3:215863/149‑1 (MQ=255)         aGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑tt                                                                                                                               <  7:195621/149‑1 (MQ=255)             aCCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCt                                                                                                                           >  5:219422/1‑149 (MQ=255)              ccAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCtt                                                                                                                          <  3:228375/149‑1 (MQ=255)                  gctgctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCC‑tttttctttgc                                                                                                                      <  4:260007/149‑8 (MQ=18)                     gctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTg                                                                                                                   >  3:262986/1‑149 (MQ=255)                          aTCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGc                                                                                                              <  2:254736/149‑1 (MQ=255)                           tCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCg                                                                                                             <  7:173361/149‑1 (MQ=255)                            cccGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCgg                                                                                                            >  8:195142/1‑149 (MQ=255)                             ccGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCgg                                                                                                            >  6:337706/1‑148 (MQ=255)                               gcgTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAg                                                                                                         >  2:155455/1‑149 (MQ=255)                                        gggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTg                                                                                                                                     <  8:239869/113‑1 (MQ=255)                                               ttCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACcgac                                                                                         >  5:175651/1‑149 (MQ=255)                                                tCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACcgac                                                                                         <  6:219422/148‑1 (MQ=255)                                                  tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGAtt                                                                                                  >  1:123031/1‑137 (MQ=255)                                                  tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGAtt                                                                                                  <  2:123031/137‑1 (MQ=255)                                                         tgtgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCggg                                                                               >  8:182401/1‑149 (MQ=255)                                                          gtgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGt                                                                              >  2:246466/1‑149 (MQ=255)                                                           tgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGtt                                                                             <  6:175651/149‑1 (MQ=255)                                                             gTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt                                                                           <  4:262986/149‑1 (MQ=255)                                                             gTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt                                                                           <  5:337706/149‑1 (MQ=255)                                                                    tGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGccc                                                                     <  1:246466/148‑1 (MQ=255)                                                                       cAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgac                                                                  <  7:195142/148‑1 (MQ=255)                                                                       cAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgac                                                                  >  1:255147/1‑148 (MQ=255)                                                                         aGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAg                                                               >  4:85452/1‑149 (MQ=255)                                                                              tAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACcgac                                                                                         >  7:98777/1‑118 (MQ=255)                                                                              tAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACcgac                                                                                         <  8:98777/118‑1 (MQ=255)                                                                                  cTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGCCAAACAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTCAATCTGGTTGCGGGAGCAGCATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACAAGgc                                                      >  2:372023/1‑149 (MQ=17)                                                                                   tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGgct                                                     >  6:154055/1‑149 (MQ=255)                                                                                   tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGgct                                                     >  6:302519/1‑149 (MQ=255)                                                                                      gATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGgctgc                                                   <  2:255147/148‑1 (MQ=255)                                                                                       aTTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGgctgct                                                  <  7:182401/148‑1 (MQ=255)                                                                                              gACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCcg                                          >  6:136570/1‑149 (MQ=255)                                                                                              gACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCcg                                          >  7:15314/1‑149 (MQ=255)                                                                                                 gTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGt                                       >  3:187154/1‑149 (MQ=255)                                                                                                   tttttCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGt                                       >  3:24490/1‑147 (MQ=255)                                                                                                      ttCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGt                                   <  1:155455/148‑1 (MQ=255)                                                                                                               gCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCggg                                                                               >  6:285160/1‑95 (MQ=255)                                                                                                               gCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCggg                                                                               <  5:285160/95‑1 (MQ=255)                                                                                                               gCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTCGAATCTGGTTGCGGGAGCCGGAt                                                                                                   <  2:321827/75‑1 (MQ=255)                                                                                                                 ttATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTc                        >  2:232842/1‑148 (MQ=255)                                                                                                                               aagaaaaGGCCACTGTTAACACAGTGGCCTTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGgct                                                     >  6:369349/4‑106 (MQ=11)                                                                                                                                 aaaaaGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTtat       <  3:85452/149‑1 (MQ=255)                                                                                                                                    aaGGCCACTGTTAACACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATg    >  1:228271/1‑149 (MQ=255)                                                                                                                                      ggCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTc  <  8:15314/149‑1 (MQ=255)                                                                                                                                       gCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGACAGACGAGCTACCAGGCTGCTCCATCCAGCGTCAGTGGTGCGGATTCTACAGCGTGTAGTTTATATGTc  >  1:310824/1‑148 (MQ=255)                                                                                                                                        ccACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTc  >  2:212088/1‑147 (MQ=255)                                                                                                                                         cACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATc                                             >  3:330590/1‑103 (MQ=35)                                                                                                                                         cACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATc                                             <  4:330590/103‑1 (MQ=35)                                                                                                                                                |                                                                                                                                            GCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTC  >  NC_002947/1380580‑1380861

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