Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,976,460	C→A	12.3%	L69I (CTC→ATC)	yifB →	ATP‑dependent chelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,976,460	0	C	A	12.3%	99.8 / 2.5	38	L69I (CTC→ATC)	yifB	ATP‑dependent chelatase
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base A (3/2);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAAACCCACCTGGCCAACGGCCTGCCTCATCTCACCCTGGTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTGGCCCTGTCGGGGAAATTGCGCCCG  >  NC_002947/5976317‑5976609                                                                                                                                                |                                                                                                                                                      tCGAAACCCACCTGGCCAACGGCCTGCCTCATCTCACCCTGGTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGccc                                                                                                                                                 >  1:517300/1‑150 (MQ=255)             tGGCCAACGGCCTGCCTCATCTCACCCTGGTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGc                                                                                                                                     <  2:517300/150‑1 (MQ=255)                           cctCATCTCACCCTGGTCGGCCTGCCGGAAACCACGATCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATggcggc                                                                                                                       >  1:520193/1‑150 (MQ=255)                              cATCTCACCCTGGTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGa                                                                                                                    <  1:50853/150‑1 (MQ=255)                              cATCTCACCCTGGTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACATCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGa                                                                                                                    >  2:121212/1‑150 (MQ=255)                              cATCTCACCCTGCTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACATCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGa                                                                                                                    >  2:425574/1‑150 (MQ=255)                                 cgcACCCTGGTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCCGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGatat                                                                                                                 <  2:148323/148‑1 (MQ=255)                                    aCCCTGGTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGAc                                                                                                              >  2:142288/1‑150 (MQ=255)                                       cTGGTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTg                                                                                                           <  1:68808/150‑1 (MQ=255)                                          gTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGcc                                                                                                        <  2:137646/150‑1 (MQ=255)                                             ggCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATc                                                                                                     <  1:588250/150‑1 (MQ=255)                                                cTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGcc                                                                                                  >  1:600630/1‑150 (MQ=255)                                                 tGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGccc                                                                                                 >  2:510034/1‑150 (MQ=255)                                                     ggAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTggg                                                                                             >  1:542589/1‑150 (MQ=255)                                                              ggTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTcgc                                                                                    >  2:543515/1‑150 (MQ=255)                                                              ggTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACATCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGATCTGGCCATCGCCCTGGGCATCCTcgc                                                                                    <  1:425574/150‑1 (MQ=255)                                                              ggTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACATCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTcgc                                                                                    <  1:121212/150‑1 (MQ=255)                                                                          cAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCa                                                                        <  1:391235/150‑1 (MQ=255)                                                                           aaGGACCGGGTGCGCAGTGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAg                                                                       >  2:412460/1‑150 (MQ=255)                                                                                ccGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAg                                                                       >  1:30831/1‑145 (MQ=255)                                                                                ccGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAg                                                                       <  2:30831/145‑1 (MQ=255)                                                                                ccGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGcc                                                                  <  2:600630/150‑1 (MQ=255)                                                                                  gggTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGa                                                                >  2:124503/1‑150 (MQ=255)                                                                                                      aaCTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATg                                            <  2:438699/150‑1 (MQ=255)                                                                                                      aaCTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATg                                            >  1:283139/1‑150 (MQ=255)                                                                                                         tCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGaa                                         <  2:283139/150‑1 (MQ=255)                                                                                                             ggCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGcc                                     <  1:326937/150‑1 (MQ=255)                                                                                                                  aaGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGc                                >  2:49308/1‑150 (MQ=255)                                                                                                                     tACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGaa                             >  2:300753/1‑150 (MQ=255)                                                                                                                        ccACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTCCACATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTg                          >  1:497654/1‑150 (MQ=255)                                                                                                                        ccACAGCGGCGCATTACCCAGAACATGGCCCCCGCCGACCTGACCAAGGATGGCGGCCGATATGACCTGGCCCTCGCCCTGGGCATCCTCGCCGCAGACGACCAGGTCCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGaagtc                          >  2:276356/1‑147 (MQ=255)                                                                                                                                  gcATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTGGCCCTGTCgg                <  2:542589/150‑1 (MQ=255)                                                                                                                                   cATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTGGCCCTGTCggg               <  1:543515/150‑1 (MQ=255)                                                                                                                                   cATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCATAGGTGCACGTTCTCGCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTGGCCCTGTCggg               <  2:177418/150‑1 (MQ=255)                                                                                                                                       aCCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTGGCCCTGTCGGGGaaa           >  2:296845/1‑150 (MQ=255)                                                                                                                                         ccAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTGGCCCTGTCGGGGAAAtt         >  2:105839/1‑150 (MQ=255)                                                                                                                                          cAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATCGCTAGCCGGACCGATAGGGAATGCCTGGGCGAACTGGCCCTGTCGGGGAAATTg        >  1:291569/1‑150 (MQ=255)                                                                                                                                              aCCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTGGCCCTGTCGGGGAAATTGCGcc    <  1:510034/150‑1 (MQ=255)                                                                                                                                                cTCGCCCCCGCAGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTGGCCCTGTCGGGGAAATTGCGCCCg  >  2:89227/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                      TCGAAACCCACCTGGCCAACGGCCTGCCTCATCTCACCCTGGTCGGCCTGCCGGAAACCACGGTCAAGGAAAGCAAGGACCGGGTGCGCAGCGCCATCGTCAACTCCGGGCTGAAGTACCCACAGCGGCGCATTACCCAGAACCTCGCCCCCGCCGACCTGCCCAAGGATGGCGGCCGATATGACCTGGCCATCGCCCTGGGCATCCTCGCCGCCGACGGCCAGGTGCCGATTGCCACCCTGACCGATATGGAATGCCTGGGCGAACTGGCCCTGTCGGGGAAATTGCGCCCG  >  NC_002947/5976317‑5976609

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