Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,813,216	A→T	6.9%	N46I (AAC→ATC)	argS →	arginine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,813,216	0	A	T	6.9%	171.3 / 3.1	58	N46I (AAC→ATC)	argS	arginine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base T (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

AAAGACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTCCAGAACACCGACG  >  NC_002947/5813083‑5813362                                                                                                                                      |                                                                                                                                                    aaaGACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGGGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGc                                                                                                                                       <  7:187696/148‑1 (MQ=255)  aaGACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTg                                                                                                                                     <  2:238517/149‑1 (MQ=255)  aaGACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTg                                                                                                                                     <  8:243806/149‑1 (MQ=255)       acCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCATCATGGCGATGTTGCTGGCGaa                                                                                                                                <  1:9143/149‑1 (MQ=255)       acCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGc                                                                                                                                               <  2:5125/134‑1 (MQ=255)       acCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGc                                                                                                                                               >  1:5125/1‑134 (MQ=255)          aTTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGAGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAgcc                                                                                                                             >  8:38208/1‑149 (MQ=255)          aTTCGCCAGCTGATCCAGCAACCCCTCACCCAACTCGTCACCGACGGTCTGCTGCCTGAAGGGCTGACGCCCGCGAACCAGGTGGAAAACGCCCGGGACAACACCCACGCCGAAATAGCCAGCAACATCGCCATGAAGCTGGCCAAgcc                                                                                                                             >  8:183472/1‑149 (MQ=255)                 aGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATCATGCTGGCCAAGccggccgg                                                                                                                       >  6:32747/1‑148 (MQ=255)                      aTCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGa                                                                                                                  <  5:341830/148‑1 (MQ=255)                         cAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑Gcc                                                                                                             >  3:194806/1‑149 (MQ=255)                         cAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCAACTTCGCCAGCATCATGGCGATGTTGCTgg                                                                                                                                    >  5:85864/1‑127 (MQ=255)                         cAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCAACTTCGCCAGCATCATGGCGATGTTGCTgg                                                                                                                                    <  6:85864/127‑1 (MQ=255)                             aaGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCcgcgc                                                                                                         <  7:226816/149‑1 (MQ=255)                               gCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGAGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGa                                                                                                       <  7:38208/149‑1 (MQ=255)                                 ccTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGGTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGAcc                                                                                                     >  2:240053/1‑149 (MQ=255)                                  cTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCt                                                                                                    >  8:78424/1‑149 (MQ=255)                                     aCCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGc                                                                                                 <  3:337303/149‑1 (MQ=255)                                     aCCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGc                                                                                                 >  1:308631/1‑149 (MQ=255)                                       ccAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGcc                                                                                                >  6:241879/1‑148 (MQ=255)                                       ccAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGcc                                                                                                >  4:69932/1‑148 (MQ=255)                                                 aCCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGccggccgg                                                                                                                       <  3:157280/116‑1 (MQ=255)                                                 aCCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGccggccgg                                                                                                                       >  4:157280/1‑116 (MQ=255)                                                 aCCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGAt                                                                                     <  1:312596/149‑1 (MQ=255)                                                  ccGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAACGCCGCGCGACCTCGCCGAAAAACTGAt                                                                                     <  1:198635/149‑1 (MQ=255)                                                  ccGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATc                                                                                    <  5:32747/149‑1 (MQ=255)                                                                  tGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGc                                                                      >  8:51916/1‑147 (MQ=255)                                                                      gggCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAgcgc                                                                <  4:132857/149‑1 (MQ=255)                                                                      gggCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAgcgc                                                                <  7:280005/149‑1 (MQ=255)                                                                       ggCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGcc                                                               <  7:78424/149‑1 (MQ=255)                                                                        gCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGcc                                                               >  8:19409/1‑148 (MQ=255)                                                                        gCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCACCGcc                                                               >  7:63282/1‑148 (MQ=255)                                                                            tCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGAc                                                            <  4:194806/147‑1 (MQ=255)                                                                              gCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATc                                                         <  7:51916/148‑1 (MQ=255)                                                                               ccGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAg                                                       >  6:122058/1‑149 (MQ=255)                                                                                         aGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTc                                                >  7:3986/1‑146 (MQ=255)                                                                                             ggAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTg                                         >  2:57257/1‑149 (MQ=255)                                                                                                aaaCGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAg                                      >  1:72801/1‑149 (MQ=255)                                                                                                    gCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGccc                                  <  1:57257/149‑1 (MQ=255)                                                                                                        gggACAAGACCCACGGCGACTTCGCCAGCAACATCGcc                                                                                                                                              >  8:154663/1‑38 (MQ=255)                                                                                                        gggACAAGACCCACGGCGACTTCGCCAGCAACATCGcc                                                                                                                                              <  7:154663/38‑1 (MQ=255)                                                                                                        gggACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGc                              <  7:19409/149‑1 (MQ=255)                                                                                                             aaGACCCACGGCGACTACGCCAGCAACATCGCCATGAGGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTCATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCt                         <  2:296160/149‑1 (MQ=255)                                                                                                               gACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTc                        >  3:270405/1‑148 (MQ=255)                                                                                                                 cccACGGCGACTTCGCCAGCAACATCGCCATGATGGTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAAc                     <  1:240053/149‑1 (MQ=255)                                                                                                                  ccACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAAct                    >  8:39841/1‑149 (MQ=255)                                                                                                                  ccACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAAct                    >  5:189818/1‑149 (MQ=255)                                                                                                                    aCGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAActtc                  >  2:177700/1‑149 (MQ=255)                                                                                                                     cGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCg                                >  8:210637/1‑134 (MQ=255)                                                                                                                     cGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCg                                <  7:210637/134‑1 (MQ=255)                                                                                                                     cGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAActtct                 >  4:155569/1‑149 (MQ=255)                                                                                                                       gCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACttcttc               <  2:308631/149‑1 (MQ=255)                                                                                                                        cGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTcc              <  1:177700/149‑1 (MQ=255)                                                                                                                        ataCTTCGCCAGCTACATCGCCATGATGCTGGCCAAGCCTGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCACTTACGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTcc              <  7:183472/147‑1 (MQ=255)                                                                                                                          gcTTGGCCAGCATCATGGCGATGTTGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACttcttc               >  5:325246/2‑146 (MQ=255)                                                                                                                            ttCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGACATTCCAGGCCCCGGCTTCCTCAACTTCTTCCAGaa          >  7:112914/1‑149 (MQ=255)                                                                                                                              cGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCCACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGGGGAGATTGCTGGCCCGGGCTTCCTCAACTGCTTCGAGcaca        >  4:250268/1‑149 (MQ=255)                                                                                                                                    cAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTCCAGAACACCGACg  <  3:155569/149‑1 (MQ=255)                                                                                                                                      |                                                                                                                                                    AAAGACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAA‑GCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTCCAGAACACCGACG  >  NC_002947/5813083‑5813362

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