Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,097,231	G→C	11.8%	intergenic (+222/+70)	PP_mr66 → / ← acsA‑I	P6/acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,097,231	0	G	C	11.8%	96.2 / 3.2	33	intergenic (+222/+70)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (15/14);  new base C (2/2);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

GCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGT  >  NC_002947/5097083‑5097362                                                                                                                                                     |                                                                                                                                     gCACAAACTGCTGACTCCGATGCTCCGAAACGTGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAAGAGTAACCGATGATTGTCGGGTAACAATGGAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAggg                                                                                                                                      <  6:268333/149‑1 (MQ=255)   acaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCGAATGCGTGAGGTCAGGGGa                                                                                                                                    >  5:27372/1‑149 (MQ=255)   acaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGa                                                                                                                                    <  2:322917/149‑1 (MQ=255)   acaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGa                                                                                                                                    >  3:243160/1‑149 (MQ=255)    caAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGAc                                                                                                                                   <  3:247130/149‑1 (MQ=255)        aTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCCAATGCGTGAGGTGAGGGGACACCt                                                                                                                               >  5:221158/1‑149 (MQ=255)         tGCTGACTCCTATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCACCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCt‑t                                                                                                                             <  6:292621/149‑1 (MQ=255)           cTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑Ta                                                                                                                            <  3:91309/148‑1 (MQ=255)              aCTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAGCCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCt                                                                                                                        >  1:85658/1‑149 (MQ=255)              aCTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAggg                                                                                                                                      >  8:197332/1‑136 (MQ=255)              aCTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAggg                                                                                                                                      <  7:197332/136‑1 (MQ=255)              aCTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGc                                                                                                                          <  5:186920/147‑1 (MQ=255)                 ccGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGa                                                                                                                      >  8:45339/1‑148 (MQ=255)                       cACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACg                                                                                                               <  6:172714/149‑1 (MQ=255)                        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGAGACCTGACGGCACTCCT‑AAGCCACAACTGACg                                                                                                               >  2:156880/1‑148 (MQ=255)                         ccGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCa                                                                                                             >  1:26587/1‑149 (MQ=255)                            aaaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTc                                                                                                          <  2:26587/149‑1 (MQ=255)                             aaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcg                                                                                                         <  1:156880/149‑1 (MQ=255)                              aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgc                                                                                                        <  2:63227/149‑1 (MQ=21)                               cGGGCCGGAATTGTTTCCTGGCGCCCCGAAGCTGTGCGGGGCCGGTCACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGCGCCCGCGAATGCGGCACTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgcg                                                                                                       <  1:35360/149‑1 (MQ=18)                                gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                      >  5:130013/1‑149 (MQ=255)                                gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                      <  7:45339/149‑1 (MQ=255)                                gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                      <  7:174104/149‑1 (MQ=255)                                gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                      <  1:272572/149‑1 (MQ=255)                                gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGGAGTGGTGTGGTCAGGGGGCAGCT‑TAGCCGGAACTGACGCCTTCGCgg                                                                                                      >  6:50671/1‑149 (MQ=255)                                    cGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGc                                                                                                    >  3:145069/1‑147 (MQ=255)                                    cGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCaa                                                                                                  <  4:243160/149‑1 (MQ=255)                                      gAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGc                                                                                                <  2:85658/149‑1 (MQ=255)                                                       cccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCGAATGCGTGAGGTCAGGGGACAGGT‑TAGGCTGGACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGt                                                                               >  7:150534/1‑149 (MQ=255)                                                                                             atgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGCAATGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCg                                         >  5:185338/1‑149 (MQ=255)                                                                                                                  aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGGATCGATCAGGt                    >  7:95089/1‑149 (MQ=255)                                                                                                                                ccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc      >  1:325905/1‑149 (MQ=255)                                                                                                                                ccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc      >  4:110225/1‑149 (MQ=255)                                                                                                                                  gcgAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTTTTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg    >  4:373404/1‑149 (MQ=18)                                                                                                                                  gcgAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTTCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg    >  5:293599/1‑149 (MQ=21)                                                                                                                                    gAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGt  <  3:110225/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                     GCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGT  >  NC_002947/5097083‑5097362

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