Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,384,954	A→G	7.8%	M26T (ATG→ACG)	proS ←	proline‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,384,954	0	A	G	7.8%	143.1 / 2.7	51	M26T (ATG→ACG)	proS	proline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base G (2/2);  total (26/25)
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.20e-01

GATGCTCGGCATCAGCACTTCCAGGGCGCCGGCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCG  >  NC_002947/1384809‑1385097                                                                                                                                                  |                                                                                                                                                 gATGCTCGGCATCAGCACTTCCAGGGCGCCGGCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑t                                                                                                                                                <  1:19949/147‑1 (MQ=255)   tGCTCGGCATCAGCACTTCCAGGGCGCCGGCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAgc                                                                                                                                            <  7:57506/149‑1 (MQ=255)     cTCGGCAACAGCACTTCCAGGGCGCCGGCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAgct                                                                                                                                           >  2:63361/1‑148 (MQ=255)      tCGGCATCAGCACTTCCAGGGCGCCGGCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAgctga                                                                                                                                         >  4:291189/1‑149 (MQ=255)                   ttCCAGGGCGCCGGCGGCGTTCATTTCCTAACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCa                                                                                                                            >  1:21663/1‑149 (MQ=255)                     ccAGGGCGCCGGCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACt                                                                                                                          >  4:221718/1‑149 (MQ=255)                      cAGGGCGCCGTCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTg                                                                                                                         >  8:48722/1‑149 (MQ=255)                       aGGGCGCCGGCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGc                                                                                                                        <  3:291189/149‑1 (MQ=255)                        gggCGCCGGCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCg                                                                                                                       >  6:56211/1‑149 (MQ=255)                              cggcggCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCgg                                                                                                                   >  2:45328/1‑147 (MQ=255)                               ggcggcGTTCATTTCCTCACGCACGACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGGTGATGGCTGATGACCACTGCGTCGGCAg                                                                                                                <  4:269839/149‑1 (MQ=255)                                  ggcgTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAgggg                                                                                                             <  1:237763/149‑1 (MQ=255)                                        cATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑Tg                                                                                                                                               <  7:162111/109‑1 (MQ=255)                                        cATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑Tg                                                                                                                                               >  8:162111/1‑109 (MQ=255)                                              cTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCATTGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTg                                                                                                   <  1:63361/148‑1 (MQ=255)                                              cTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTgg                                                                                                  <  7:48722/148‑1 (MQ=255)                                              cTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTgg                                                                                                  <  1:130409/148‑1 (MQ=255)                                                 aCGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCAGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGt                                                                                                <  1:52748/147‑1 (MQ=255)                                                  cGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGc                                                                                             <  6:276393/149‑1 (MQ=255)                                                   gCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCG‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGc                                                                                             >  7:7315/1‑148 (MQ=255)                                                   gCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGc                                                                                             >  4:159857/1‑148 (MQ=255)                                                        acGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCaa                                                                                       <  5:56211/149‑1 (MQ=255)                                                         cGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCG‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCaaa                                                                                      <  8:7315/149‑1 (MQ=255)                                                          ggCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAat                                                                                     <  3:128395/149‑1 (MQ=255)                                                           gCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCg                                                                                                                    <  1:156100/117‑1 (MQ=255)                                                           gCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCg                                                                                                                    >  2:156100/1‑117 (MQ=255)                                                                 aCCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTg                                                                                                   >  3:148658/1‑128 (MQ=255)                                                                 aCCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTg                                                                                                   <  4:148658/128‑1 (MQ=255)                                                                   cTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTgg                                                                            >  5:91094/1‑149 (MQ=255)                                                                        gcATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTgcg                                                                        >  6:47455/1‑148 (MQ=255)                                                                               ccGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCG‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGc                                                                                             <  8:114620/120‑1 (MQ=255)                                                                               ccGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCG‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGc                                                                                             >  7:114620/1‑120 (MQ=255)                                                                                 gCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGc                                                              >  1:123180/1‑149 (MQ=255)                                                                                     gCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCgtg                                                          <  3:159857/149‑1 (MQ=255)                                                                                       ccATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTc                                                        >  2:68125/1‑149 (MQ=255)                                                                                       ccATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTc                                                        <  1:68125/149‑1 (MQ=255)                                                                                        cATCGGCAGCCAGGTGTACAGGCAGGAGGCCAGTTTGCGGAGCATGCCGGAAAGCAGCA‑TGATCTTATTGCTGATGAACACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGTTTAGCCTTGTc                                                        >  4:73316/1‑148 (MQ=255)                                                                                             gCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGc                                                  >  3:265364/1‑149 (MQ=255)                                                                                              cAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCgtg                                                          <  4:118368/140‑1 (MQ=255)                                                                                              cAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCgtg                                                          >  3:118368/1‑140 (MQ=255)                                                                                                        tACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGtt                                        <  5:47455/148‑1 (MQ=255)                                                                                                            ggCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGGTGAAAt                                   <  6:91094/149‑1 (MQ=255)                                                                                                               cGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGc                                >  5:299981/1‑149 (MQ=255)                                                                                                                      cAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCgtgt                                                         >  3:17746/1‑117 (MQ=255)                                                                                                                      cAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCgtgt                                                         >  7:57616/1‑117 (MQ=255)                                                                                                                      cAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATCGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGt                                         <  2:96734/133‑1 (MQ=255)                                                                                                                      cAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGAGGCGAGCAAATATCGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGt                                         >  1:96734/1‑133 (MQ=255)                                                                                                                              ggATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACggg                 >  8:45325/1‑149 (MQ=255)                                                                                                                               gATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACgggg                <  3:221718/149‑1 (MQ=255)                                                                                                                                 tCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCTGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGGGCGCATGGTTAGCCGTGGCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACgggggg              <  1:45328/149‑1 (MQ=255)                                                                                                                                     gCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGaaa          >  6:242540/1‑149 (MQ=255)                                                                                                                                             gcagca‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGcc                               >  8:4083/1‑120 (MQ=255)                                                                                                                                             gcagca‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCg  >  5:171513/1‑149 (MQ=255)                                                                                                                                             gcagca‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACCTTGAAATAGcc                               <  7:4083/120‑1 (MQ=255)                                                                                                                                               agca‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAAt                                   >  5:306174/1‑114 (MQ=255)                                                                                                                                               agca‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAAt                                   <  6:306174/114‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 GATGCTCGGCATCAGCACTTCCAGGGCGCCGGCGGCGTTCATTTCCTCACGCACCACGGCTTCGACCTTGCGCATTACCCGCAAGCCCATCGGCAGCCAGGTGTACAGGCCGGAGGCCAGTTTGCGGATCATGCCGGCACGCAGCA‑TGAGCTGATGGCTGATGACCACTGCGTCGGCAGGGGTTTCTTTCTGGGTGGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCG  >  NC_002947/1384809‑1385097

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