Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,151,933:1	+C	8.9%	intergenic (+123/+28)	PP_3655 → / ← PP_3656	cytosine/purine/uracil/thiamine/allantoin permease family protein/aromatic compound‑specific porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,151,933	1	.	C	8.9%	165.9 / 10.8	45	intergenic (+123/+28)	PP_3655/PP_3656	cytosine/purine/uracil/thiamine/allantoin permease family protein/aromatic compound‑specific porin
Reads supporting (aligned to +/- strand):  ref base . (17/24);  new base C (2/2);  total (19/26)
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.55e-01

GCGCCAAGCTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTC  >  NC_002947/4151785‑4152075                                                                                                                                                       |                                                                                                                                               gcgcCAAGCTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCgg                                                                                                                                                 >  2:83283/1‑149 (MQ=255)         cTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAg                                                                                                                                          >  4:290757/1‑147 (MQ=255)                cGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAACGGTCCGGCCTTCTGAa‑gccggc                                                                                                                                 <  3:12531/149‑1 (MQ=255)                      cTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGttt                                                                                                                           >  7:144165/1‑149 (MQ=255)                      cTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGttt                                                                                                                           <  3:290757/149‑1 (MQ=255)                         aaTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTc                                                                                                                          <  7:34169/147‑1 (MQ=255)                          aTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑Ctt                                                                                                                                             >  4:17390/1‑127 (MQ=255)                          aTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑Ctt                                                                                                                                             <  3:17390/127‑1 (MQ=255)                           tCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGtt                                                                                                                            <  4:222895/143‑1 (MQ=255)                           tCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGtt                                                                                                                            >  3:222895/1‑143 (MQ=255)                                  tGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGaa                                                                                                               <  8:144165/149‑1 (MQ=255)                                    gCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGt                                                                                                             <  8:262632/149‑1 (MQ=255)                                          cGTTGCAAAAGATGAAACGGCTGGGGGCAGGCGTCCGAACGCCTGGGGCAATTCACGCGGAAAACACCGAACCCTTCGATTGGGAATAAATGGCCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGaa                                                                                                         <  1:83283/148‑1 (MQ=255)                                          agTTTCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAg                                                                                                        <  8:124991/147‑1 (MQ=255)                                            ttGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGgtg                                                                                                     <  4:77149/149‑1 (MQ=255)                                               cAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACTTTTGTTTAAATTCATGCGGAAAACAGCTAACCTTTAGATTGTTAATAAATTT‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTa                                                                                                   <  6:145871/148‑1 (MQ=255)                                                  aaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAgcgc                                                                                               <  8:98184/149‑1 (MQ=255)                                                                 ggTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGTGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAg                                                                                                  <  7:146423/131‑1 (MQ=255)                                                                 ggTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGTGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAg                                                                                                  >  8:146423/1‑131 (MQ=255)                                                                 gcaccAGGCGTTCGGACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAttt                                                                                >  2:172956/5‑149 (MQ=255)                                                                  caccAGGCGTTCGGACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAtttt                                                                               >  6:288279/4‑149 (MQ=255)                                                                            tcggaCGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCt                                                                     >  8:41469/5‑149 (MQ=255)                                                                                 cGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAg                                                                >  4:65323/1‑149 (MQ=255)                                                                                     tGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCa                                                                             >  1:259055/1‑132 (MQ=255)                                                                                     tGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCa                                                                             <  2:259055/132‑1 (MQ=255)                                                                                             ttCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGc                                                    <  7:41469/149‑1 (MQ=255)                                                                                                   cGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCg                                              >  7:212479/1‑149 (MQ=255)                                                                                                    ggAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGt                                             <  4:158919/149‑1 (MQ=255)                                                                                                       aaaCAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGa                                          <  3:65323/149‑1 (MQ=255)                                                                                                       aaaCAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGATATTTAGAAGTTGAAGGTGTAGAGCATGATCAGGCGATTTTAATTTGCGATGTCAGCGTAGGTAGAGCGGTTCGTg                                            >  7:240979/1‑147 (MQ=255)                                                                                                        aaCAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTAAGAAGGCAGGCCCGTTTCGCGATTTAGAAGTTTAACGTGGAGCACATGATCAGCCTATTTTCATTAGCGCAGTCAGCGTAGTTCGAGCGGTTCGTGGAg                                         >  6:46846/1‑149 (MQ=255)                                                                                                            gCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTAc                                     <  8:212479/149‑1 (MQ=255)                                                                                                                 cccTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCAc                                <  1:172956/149‑1 (MQ=255)                                                                                                                  ccTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACt                               >  7:112930/1‑149 (MQ=255)                                                                                                                    ttCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCa                                                                                              >  1:267152/1‑84 (MQ=255)                                                                                                                    ttCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCa                                                                                              <  2:267152/84‑1 (MQ=255)                                                                                                                    ttCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTg                              >  3:162969/1‑148 (MQ=255)                                                                                                                         ttGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacga                        <  4:162969/149‑1 (MQ=255)                                                                                                                           tgtaaTAAATGG‑CGTAAAATGGGCCGG‑CTTCATAAGGTCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCACATGATGAGGCGATTTTCATTATCGCTGTCAGCGTAGGTCGATCTGTTCGTGGAGTTACGCCACTgaaacgaa                       <  5:46846/145‑1 (MQ=255)                                                                                                                                         gAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaac                     <  8:130852/137‑1 (MQ=255)                                                                                                                                         gAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaac                     >  7:130852/1‑137 (MQ=255)                                                                                                                                           aaaaTGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGAc       <  2:244276/149‑1 (MQ=255)                                                                                                                                              acgggCCGGCCTTCTGAA‑GCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAacac                   >  3:304123/3‑134 (MQ=255)                                                                                                                                              acgggCCGGCCTTCTGAA‑GCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAacac                   <  4:304123/132‑1 (MQ=255)                                                                                                                                               tGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCAct   <  5:288279/149‑1 (MQ=255)                                                                                                                                                gggCCGGCCTTCTGAA‑GCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCActc  >  1:296071/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                               GCGCCAAGCTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTC  >  NC_002947/4151785‑4152075

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