Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,151,938	A→T	8.9%	intergenic (+128/+23)	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,938	0	A	T	8.9%	125.7 / 3.2	45	intergenic (+128/+23)	PP_3655/PP_3656	cytosine/purine/uracil/thiamine/allantoin permease family protein/aromatic compound‑specific porin
Reads supporting (aligned to +/- strand):  ref base A (17/24);  new base T (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 = 9.33e-01

CTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTC  >  NC_002947/4151793‑4152075                                                                                                                                                    |                                                                                                                                          cTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAg                                                                                                                                          >  4:290757/1‑147 (MQ=255)        cGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAACGGTCCGGCCTTCTGAa‑gccggc                                                                                                                                 <  3:12531/149‑1 (MQ=255)              cTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGttt                                                                                                                           <  3:290757/149‑1 (MQ=255)              cTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGttt                                                                                                                           >  7:144165/1‑149 (MQ=255)                 aaTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTc                                                                                                                          <  7:34169/147‑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                     >  7:130852/1‑137 (MQ=255)                                                                                                                                 gAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaac                     <  8:130852/137‑1 (MQ=255)                                                                                                                                   aaaaTGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGAc       <  2:244276/149‑1 (MQ=255)                                                                                                                                      acgggCCGGCCTTCTGAA‑GCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAacac                   <  4:304123/132‑1 (MQ=255)                                                                                                                                      acgggCCGGCCTTCTGAA‑GCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAacac                   >  3:304123/3‑134 (MQ=255)                                                                                                                                       tGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCAct   <  5:288279/149‑1 (MQ=255)                                                                                                                                        gggCCGGCCTTCTGAA‑GCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCActc  >  1:296071/1‑149 (MQ=255)                                                                                                                                                 ttCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaa                      <  2:214123/121‑1 (MQ=255)                                                                                                                                                 ttCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaa                      >  1:214123/1‑121 (MQ=255)                                                                                                                                                    |                                                                                                                                          CTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGG‑CGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTC  >  NC_002947/4151793‑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: ‑