Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,200,456	A→G	10.4%	F183L (TTC→CTC)	yafJ ←	amidotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,200,456	0	A	G	10.4%	124.1 / 5.2	48	F183L (TTC→CTC)	yafJ	amidotransferase
Reads supporting (aligned to +/- strand):  ref base A (19/24);  new base G (2/3);  total (21/27)
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.40e-01

CTCGGTCAGGGCCTCGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGT  >  NC_002947/5200310‑5200599                                                                                                                                                    |                                                                                                                                                ctcGGTCAGGGCCTCGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGagcagaag                                                                                                                                                <  5:76665/148‑1 (MQ=255)    ggTCAGGGCCTCGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGagcagaagcag                                                                                                                                             <  8:165968/148‑1 (MQ=255)     gTCAGGGCCTCGGTGGCAATCACCGTGACCACGTCGTTGGCGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAgcagaagcagaa                                                                                                                                           >  2:309990/1‑149 (MQ=255)      tCAGGGCCTCGGTGGCAATCAGCGTGACCACGTAGTTGGGGGTGGTTTCGGTATGAAAATAGGCGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGAGCGTAAGGTTCACCAGCGTTGTGCAgcagaagcagaa                                                                                                                                           >  8:234901/1‑148 (MQ=255)        aGGGCCTCGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAGGCAGAACa                                                                                                                                         <  6:115533/148‑1 (MQ=255)        aGGGCCTCGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACa                                                                                                                                         <  7:70538/148‑1 (MQ=255)                 gTGGCAATTACCGTGACCTCGTCGTTGGGGGTGGTTTCGGTATTAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGc                                                                                                                                <  8:46548/148‑1 (MQ=255)                  tGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCa                                                                                                                              <  5:298377/149‑1 (MQ=255)                      aaTCACCGTGACCACGTCGATCGGGCTCGATTCGCTATGAAAATCGACGATCAGGTCGACA‑TCCAACAACCGTGCGGCACCGAATGGCGCACGGCCGGTAATGTCCACCAGCTTTCTCGAGCAGAAGCAGAACAGCCATTCTCCAACGc                                                                                                                          >  6:240487/1‑149 (MQ=255)                          aCCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGTTCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGa                                                                                                                       <  8:285022/148‑1 (MQ=255)                          aCCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAg                                                                                                                      <  3:2204/149‑1 (MQ=255)                             gTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCAt                                                                                                                   <  1:48765/149‑1 (MQ=255)                                aCCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGc                                                                                                                 <  1:309990/148‑1 (MQ=255)                                  cACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACAGTCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACAAGCTTTGTCGAGCAGAGGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAg                                                                                                               >  1:243320/1‑149 (MQ=255)                                  cACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGt                                                                                                              >  3:97972/1‑149 (MQ=255)                                  cACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTGCAACCGTGCGGCACAGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGt                                                                                                              >  1:87641/1‑149 (MQ=255)                                        gtTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTc                                                                                                                                  <  8:235417/123‑1 (MQ=255)                                        gtTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTc                                                                                                                                  >  7:235417/1‑123 (MQ=255)                                                        gTATGAAAATCGACGATCAGGTCTACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAGGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTa                                                                                        <  2:243320/149‑1 (MQ=255)                                                        gTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGc                                                                                                                     <  6:201132/120‑1 (MQ=255)                                                        gTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGc                                                                                                                     >  5:201132/1‑120 (MQ=255)                                                         tATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAg                                                                                       >  6:37799/1‑149 (MQ=255)                                                              aaaTCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCgg                                                                                   <  2:87641/148‑1 (MQ=255)                                                                        tCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGtt                                                                                                             >  4:142988/1‑112 (MQ=255)                                                                            gTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCa                                                                    <  6:40225/149‑1 (MQ=255)                                                                               gACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCa                                                                 >  1:153815/1‑149 (MQ=255)                                                                                 cA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACt                                                               >  7:24677/1‑149 (MQ=255)                                                                                  a‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTg                                                              >  1:199292/1‑149 (MQ=255)                                                                                      cTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGgc                                                            <  4:97972/148‑1 (MQ=255)                                                                                      cTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGgc                                                            <  7:103013/148‑1 (MQ=255)                                                                                       ttCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGgcag                                                          >  4:18914/1‑149 (MQ=255)                                                                                         cAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGgcag                                                          <  1:200672/147‑1 (MQ=255)                                                                                                     aCCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCcaccg                                            >  5:139609/1‑149 (MQ=255)                                                                                                      ccGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCa                                                                                                                >  7:268912/1‑80 (MQ=255)                                                                                                      ccGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCa                                                                                                                <  8:268912/80‑1 (MQ=255)                                                                                                       cGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGgcagca                                                        >  8:114310/1‑135 (MQ=255)                                                                                                       cGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGgcagca                                                        <  7:114310/135‑1 (MQ=255)                                                                                                           tGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCcaccggcaccg                                      <  8:24677/149‑1 (MQ=255)                                                                                                               gcACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCg                                                                        >  7:201920/1‑111 (MQ=255)                                                                                                               gcACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCg                                                                        <  8:201920/111‑1 (MQ=255)                                                                                                                  cggcgtGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCggg                               <  2:153815/149‑1 (MQ=255)                                                                                                                              tgCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGc                    <  5:37799/148‑1 (MQ=255)                                                                                                                                 aCCAGCTTTGTCTAGCAGAAGCAGAACAGCCAGTCGCCATCGCGGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGAt                <  5:240487/149‑1 (MQ=255)                                                                                                                                   cAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGc              >  1:208506/1‑149 (MQ=255)                                                                                                                                     gCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCgg            <  3:293014/149‑1 (MQ=255)                                                                                                                                       tttGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGtt          >  8:204549/1‑149 (MQ=255)                                                                                                                                         tGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACCCCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTg         <  7:234901/148‑1 (MQ=255)                                                                                                                                         tGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTc                                  <  3:56080/123‑1 (MQ=255)                                                                                                                                         tGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTc                                  >  4:56080/1‑123 (MQ=255)                                                                                                                                           tCGAGCAGAGGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGa        <  1:111446/147‑1 (MQ=255)                                                                                                                                           tCGAGCAGAGGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGa        >  2:111446/1‑147 (MQ=255)                                                                                                                                             gAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGc      >  6:84097/1‑147 (MQ=255)                                                                                                                                              agcagaagcagaACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAg   >  5:73796/1‑149 (MQ=255)                                                                                                                                               gcagaagcagaaCAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGt  <  2:199292/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                CTCGGTCAGGGCCTCGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACA‑TCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGT  >  NC_002947/5200310‑5200599

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