Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,673,184	G→T	9.4%	M166I (ATG→ATT) ‡	PP_2340 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,673,184	0	G	T	9.4%	158.4 / 7.3	53	M166I (ATG→ATT) ‡	PP_2340	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base T (3/2);  total (28/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 = 4.63e-01

TGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCG  >  NC_002947/2673039‑2673327                                                                                                                                                    |                                                                                                                                                tGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATggg                                                                                                                                               >  6:12060/1‑148 (MQ=255)   gAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAg                                                                                                                                            >  7:119020/1‑149 (MQ=255)           tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCa                                                                                                                                      <  2:97563/147‑1 (MQ=255)           tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCa                                                                                                                                      <  1:4811/147‑1 (MQ=255)                    gtTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGt                                                                                                                              >  2:315350/1‑146 (MQ=255)                    gtTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGa                                                                                                                            <  1:315350/148‑1 (MQ=255)                          gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                     >  3:193593/1‑149 (MQ=255)                           ccGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCggg                                                                                                                      <  3:252721/147‑1 (MQ=255)                           ccGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCggg                                                                                                                      >  4:252721/1‑147 (MQ=255)                                gTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGaa                                                                                                                  <  4:259766/146‑1 (MQ=255)                                gTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGaa                                                                                                                  >  3:259766/1‑146 (MQ=255)                                      tGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGTGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGcc                                                                                                         >  6:161867/1‑149 (MQ=255)                                          tGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCa                                                                                                                                      <  4:213371/116‑1 (MQ=255)                                          tGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCa                                                                                                                                      >  3:213371/1‑116 (MQ=255)                                          tGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGgtg                                                                                                     >  5:237562/1‑149 (MQ=255)                                               aCTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCTGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTc                                                                                                <  8:119020/149‑1 (MQ=255)                                                   ccccAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCaa                                                                                            >  8:223334/1‑149 (MQ=255)                                                      cAATTGCGCGATGCGATCG‑ATCGCGAGGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAAcc                                                                                          >  1:172139/1‑149 (MQ=255)                                                       gaactcgcgAT‑CGATCGCATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGc                                                                                        >  2:24024/6‑149 (MQ=255)                                                              cgatgcgatCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCa                                                                                 <  5:12060/149‑1 (MQ=255)                                                               gatgcgatCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAg                                                                                <  3:40214/149‑1 (MQ=255)                                                                      tcg‑atcgCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                         <  4:193593/149‑1 (MQ=255)                                                                        g‑atcgCGA‑GTTCTTCCTGCTTTCAGAGGCACACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTg                                                                        <  7:160423/148‑1 (MQ=255)                                                                          atcgCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGa                                                                      <  2:172139/149‑1 (MQ=255)                                                                           tcgCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                     >  2:269212/1‑149 (MQ=255)                                                                                   ttcttcCTGCTTTCAGCGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAAcc                                                              <  7:31357/149‑1 (MQ=255)                                                                                       tcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCg                                                                  <  3:243169/141‑1 (MQ=255)                                                                                       tcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCg                                                                  >  4:243169/1‑141 (MQ=255)                                                                                        ccTGCTTTCAGTGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTgg                                                          <  5:161867/148‑1 (MQ=255)                                                                                         cTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGa                                                         <  7:223334/148‑1 (MQ=255)                                                                                           gCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGACACTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGTATGCGAACCTTGGACaa                                                      >  6:85634/1‑148 (MQ=255)                                                                                             tttCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAAc                                                               <  3:85613/138‑1 (MQ=255)                                                                                             tttCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAAc                                                               >  4:85613/1‑138 (MQ=255)                                                                                               tCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTcc                                                  >  3:101893/1‑149 (MQ=255)                                                                                                           tACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGc                                      <  4:101893/149‑1 (MQ=255)                                                                                                               aGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGa                                  >  5:124642/1‑149 (MQ=255)                                                                                                               aGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGa                                  >  6:44965/1‑149 (MQ=255)                                                                                                                    tataTCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCa                                           <  1:120102/135‑1 (MQ=255)                                                                                                                    tataTCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCa                                           >  2:120102/1‑135 (MQ=255)                                                                                                                    tataTCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCt                               <  2:48109/147‑1 (MQ=255)                                                                                                                    tataTCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCt                               >  1:48109/1‑147 (MQ=255)                                                                                                                    tataTCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTgg                             >  6:166850/1‑149 (MQ=255)                                                                                                                     atatCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTcc                                                  <  8:5418/127‑1 (MQ=255)                                                                                                                     atatCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTcc                                                  >  7:5418/1‑127 (MQ=255)                                                                                                                     atatCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGt                            >  6:282380/1‑149 (MQ=255)                                                                                                                        tCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                          <  6:124642/148‑1 (MQ=255)                                                                                                                           cccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                          >  3:85443/1‑145 (MQ=255)                                                                                                                                  cTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACtt               >  7:130804/1‑149 (MQ=255)                                                                                                                                   gcccTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGt                            >  2:180312/3‑135 (MQ=255)                                                                                                                                   gcccTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGt                            <  1:180312/133‑1 (MQ=255)                                                                                                                                   gcccTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTg              >  4:173835/3‑149 (MQ=255)                                                                                                                                         cgcgATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt         <  6:237562/148‑1 (MQ=255)                                                                                                                                          cctctGCCAGGGGCAGCCCTCAGGGCCCGAGAGCCTGGCGAATTTGTCGTCGGGGTATTAGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt         <  4:18337/143‑1 (MQ=255)                                                                                                                                               gCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCg  >  4:154162/1‑149 (MQ=255)                                                                                                                                                    gggCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                                      <  3:240773/60‑1 (MQ=255)                                                                                                                                                    gggCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                                      >  4:240773/1‑60 (MQ=255)                                                                                                                                                    |                                                                                                                                                TGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGA‑GTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCG  >  NC_002947/2673039‑2673327

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