Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,380,725	A→C	7.7%	intergenic (‑113/+22)	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,725	0	A	C	7.7%	157.2 / 4.7	52	intergenic (‑113/+22)	PP_t67/PP_t66	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base C (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 = 5.82e-01

GCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCTAGGCTT  >  NC_002947/1380580‑1380874                                                                                                                                                  |                                                                                                                                                       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)                                                                                                                                                  aaCAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCT‑CATCCCGCGTCACTGGGGCGGATTCTACAGCGTTTAGTTTATATGTAAAGCGCAAGGCtt  >  6:16517/1‑149 (MQ=17)                                                                                                                                                  |                                                                                                                                                       GCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCC‑TTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCTAGGCTT  >  NC_002947/1380580‑1380874

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