Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,303,967	T→C	8.9%	G201G (GGT→GGC)	cysS →	cysteine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,303,967	0	T	C	8.9%	149.2 / 5.4	56	G201G (GGT→GGC)	cysS	cysteine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (23/28);  new base C (3/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

CTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCA  >  NC_002947/3303827‑3304110                                                                                                                                             |                                                                                                                                                cTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCggg                                                                                                                                         <  1:124056/149‑1 (MQ=255)    tCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGc                                                                                                                                        <  2:87177/147‑1 (MQ=255)    tCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGCGACCCGGGcc                                                                                                                                       <  6:227716/148‑2 (MQ=255)         cAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGcaca                                                                                                                                 >  4:5587/1‑149 (MQ=255)              ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAg                                                                                                                            <  7:40255/149‑1 (MQ=255)              ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAg                                                                                                                            <  7:100503/149‑1 (MQ=255)              ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGCGACCCGGGc                                                                                                                                        <  5:117997/137‑1 (MQ=255)              ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGCGACCCGGGc                                                                                                                                        >  6:117997/1‑137 (MQ=255)               gATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGt                                                                                                                           <  3:69328/149‑1 (MQ=255)                aTCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGAC‑TAGCCAAGGATGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGTGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCTCGGGCTGGCACATCGAGTGc                                                                                                                         <  2:251375/149‑1 (MQ=255)                       aCCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTc                                                                                                                   <  3:106645/149‑1 (MQ=255)                        ccTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTc                                                                                                                   <  5:143394/148‑1 (MQ=255)                           gCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                                                                                               >  3:292478/1‑149 (MQ=255)                              tATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTg                                                                                                                                      <  3:10121/123‑1 (MQ=255)                              tATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTg                                                                                                                                      >  4:10121/1‑123 (MQ=255)                               aTCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGa                                                                                                                             <  6:61467/131‑1 (MQ=255)                               aTCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGa                                                                                                                             >  5:61467/1‑131 (MQ=255)                                tCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCAcc                                                                                                          >  6:44843/1‑149 (MQ=255)                                tCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGGCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCAc                                                                                                           <  4:9294/148‑1 (MQ=255)                                   gTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACctgc                                                                                                       <  4:292478/149‑1 (MQ=255)                                     gCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACctgctg                                                                                                     >  3:72721/1‑149 (MQ=255)                                          cATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTgg                                                                                                <  3:285293/149‑1 (MQ=255)                                                 gTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                                                                                               <  3:37686/127‑1 (MQ=255)                                                 gTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                                                                                               >  4:37686/1‑127 (MQ=255)                                                      cgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCAt                                                                                                                 >  3:150417/1‑120 (MQ=255)                                                      cgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCAt                                                                                                                 <  4:150417/120‑1 (MQ=255)                                                       gaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGa                                                                                   >  5:104405/1‑149 (MQ=255)                                                                     cccGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggcg                                                                     >  8:61428/1‑149 (MQ=255)                                                                     cccGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggcg                                                                     >  3:220542/1‑149 (MQ=255)                                                                        gCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGGGCTCGGTCATGGCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCa                                                                  <  3:308258/149‑1 (MQ=255)                                                                              cTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGAcc                                                            >  6:14358/1‑149 (MQ=255)                                                                                   tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCgaga                                                                                          <  8:308720/114‑1 (MQ=255)                                                                                   tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCgaga                                                                                          >  7:308720/1‑114 (MQ=255)                                                                                   tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCga                                                                                            <  4:22531/112‑1 (MQ=255)                                                                                   tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCga                                                                                            >  3:22531/1‑112 (MQ=255)                                                                                   tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAg                                                       >  7:210756/1‑149 (MQ=255)                                                                                   tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCACGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggc                                                                      >  5:240935/1‑134 (MQ=255)                                                                                   tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCACGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggc                                                                      <  6:240935/134‑1 (MQ=255)                                                                                             gggCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGAcc                                                            >  1:242134/1‑134 (MQ=255)                                                                                             gggCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGAcc                                                            <  2:242134/134‑1 (MQ=255)                                                                                                   cAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCacgag                                        <  6:104405/148‑1 (MQ=255)                                                                                                                ccAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCa                          <  5:56244/149‑1 (MQ=255)                                                                                                                      tGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGc                      >  6:169554/1‑147 (MQ=255)                                                                                                                          aaTCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAggcg                <  4:220542/149‑1 (MQ=255)                                                                                                                           aTCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTgg                                                         <  2:313729/107‑1 (MQ=255)                                                                                                                           aTCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTgg                                                         >  1:313729/1‑107 (MQ=255)                                                                                                                             cGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAggcggc              >  6:78167/1‑148 (MQ=255)                                                                                                                               ccATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGTTTGCCCAGAGCGAggcg                <  5:65493/144‑1 (MQ=255)                                                                                                                               ccATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGTTTGCCCAGAGCGAggcg                >  6:65493/1‑144 (MQ=255)                                                                                                                                 tgccagCCCGGGCGACCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCGTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGCTGGCCCAGAGCGAGGCGGCCCCCg         >  6:200191/6‑149 (MQ=255)                                                                                                                                 aTGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTg                                                          >  7:15170/1‑100 (MQ=255)                                                                                                                                 aTGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTg                                                          <  8:15170/100‑1 (MQ=255)                                                                                                                                    cagCCCGGGCGACCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGc       >  5:264721/3‑148 (MQ=255)                                                                                                                                     ggCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCaa     <  5:44843/149‑1 (MQ=255)                                                                                                                                       cccGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAg    >  7:153036/1‑148 (MQ=255)                                                                                                                                        ccGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCa  <  5:14358/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                CTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCA  >  NC_002947/3303827‑3304110

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