Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,731,969	A→T	7.4%	L196* (TTG→TAG) ‡	hutI ←	imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,731,969	0	A	T	7.4%	160.9 / 3.2	54	L196* (TTG→TAG) ‡	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base T (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTT  >  NC_002947/5731831‑5732117                                                                                                                                           |                                                                                                                                                     gATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCaa                                                                                                                                            <  5:241008/149‑1 (MQ=255) gATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCaa                                                                                                                                            <  5:224997/149‑1 (MQ=255)  aTGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCTTGGCCGGcgaa                                                                                                                                           <  4:121636/149‑4 (MQ=255)  aTGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAAc                                                                                                                                           <  8:101525/149‑1 (MQ=255)    gACAGTTGTTCGTCATGCAGCTTGACCGGCAGGCCCAGTTCGCTCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCATCGCCGGCCTTGGCCGGcgaag                                                                                                                                          <  7:318152/148‑5 (MQ=255)    gACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCGGTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGGATCCACCGCCACCACCAGGCCTTCGCCGGCCACGGCCGGCAACAt                                                                                                                                         >  4:270243/1‑149 (MQ=255)       aGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAtctc                                                                                                                                      <  5:145014/149‑1 (MQ=255)        gttgttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCg                                                                                                                                     >  5:308195/1‑149 (MQ=255)        gttgttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCg                                                                                                                                     <  3:270243/149‑1 (MQ=255)         ttgttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGt                                                                                                                                    <  4:50327/149‑1 (MQ=255)            ttCGGGATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAAGAGGCGCTCGACCTGGGCCGGCGAGAAGGCGAGGTGTTCGCAGAAGTCATGCACCCGATCCACCACGCCTTCGCCGGCCATGGCCGTCAACATCTCGTca                                                                                                                                  <  2:139822/148‑1 (MQ=255)               ggCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTca                                                                                                                                  >  7:55002/1‑145 (MQ=255)                gCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAt                                                                                                                             <  4:228079/149‑1 (MQ=255)                  aTGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAtgt                                                                                                                           <  6:308195/149‑1 (MQ=255)                  aTGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAtgt                                                                                                                           <  8:314606/149‑1 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAt                                                                                                                             >  7:73161/1‑138 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAt                                                                                                                             <  8:73161/138‑1 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATg                                                                                                                   <  2:108459/148‑1 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATg                                                                                                                   >  5:311983/1‑148 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGCCAACATCTCGTCACAGATGTGCGCGATg                                                                                                                   <  1:209420/148‑1 (MQ=255)                               ggCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtc                                                                                                              >  3:26457/1‑149 (MQ=255)                                     cccAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGa                                                                                                                     >  1:145261/1‑136 (MQ=255)                                       cAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGccc                                                                                                      >  8:2236/1‑149 (MQ=255)                                                cgcCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCACGTCACACATGTGCGCGATGTAGTCGTCGGCCAGACCcgcg                                                                                              >  4:153288/1‑148 (MQ=255)                                                          aaCAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCgggg                                                                                    <  8:55002/148‑1 (MQ=255)                                                                cGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTTGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGc                                                                             >  8:239530/1‑149 (MQ=255)                                                                cGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGc                                                                             >  8:239529/1‑149 (MQ=255)                                                                 gCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCg                                                                            >  3:172870/1‑149 (MQ=255)                                                                           ggCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCa                                                                  <  7:239529/149‑1 (MQ=255)                                                                           ggCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCa                                                                  <  7:2236/149‑1 (MQ=255)                                                                           ggCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCa                                                                  <  7:239530/149‑1 (MQ=255)                                                                            gCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAg                                                                 <  4:172870/149‑1 (MQ=255)                                                                                  gagaACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGt                                                           >  1:240415/1‑149 (MQ=255)                                                                                      aCGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgc                                                        <  4:26457/148‑1 (MQ=255)                                                                                        gcgAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTgc                                                                         >  8:201878/1‑129 (MQ=255)                                                                                        gcgAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTgc                                                                         <  7:201878/129‑1 (MQ=255)                                                                                          gAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgcacgca                                                   >  5:365640/1‑149 (MQ=255)                                                                                              tgtTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGt                                               >  3:297256/1‑149 (MQ=255)                                                                                                     aGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTg                                                          >  1:309620/1‑131 (MQ=255)                                                                                                     aGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTg                                                          <  2:309620/131‑1 (MQ=255)                                                                                                         ggCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCt                                    >  5:314586/1‑149 (MQ=255)                                                                                                              ccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgtcg                               >  6:39399/1‑149 (MQ=255)                                                                                                                accgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGc                             >  6:22033/1‑149 (MQ=255)                                                                                                                    catccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgt                                 <  3:201498/141‑1 (MQ=255)                                                                                                                    catccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgt                                 >  4:201498/1‑141 (MQ=255)                                                                                                                         accaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGgcggcg                    >  4:76048/1‑149 (MQ=255)                                                                                                                          ccaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGcggcgg                   <  6:71538/149‑1 (MQ=255)                                                                                                                              gatgttgCCGGCCTTGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGGCGGCCAGGCGGCGGGCCa               >  7:293451/7‑149 (MQ=255)                                                                                                                                cTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACt             >  1:285731/1‑149 (MQ=255)                                                                                                                                  ttgCCGGCCTTGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCg           >  6:193576/3‑149 (MQ=255)                                                                                                                                    gCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCgca         <  3:76048/149‑1 (MQ=255)                                                                                                                                    gCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCgc          <  5:22033/148‑1 (MQ=255)                                                                                                                                        gCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCAt     >  7:178659/1‑149 (MQ=255)                                                                                                                                         ccAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATc    <  5:250596/149‑1 (MQ=255)                                                                                                                                           aaGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTCCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCtt  >  3:122813/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                     GATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTT  >  NC_002947/5731831‑5732117

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