Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,159,147	T→C	14.8%	intergenic (+51/‑158)	PP_3661 → / → PP_3662	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,159,147	0	T	C	14.8%	64.3 / 3.3	27	intergenic (+51/‑158)	PP_3661/PP_3662	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/13);  new base C (2/2);  total (12/15)
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.91e-01

GCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAG  >  NC_002947/4158998‑4159293                                                                                                                                                      |                                                                                                                                                   gCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGATCCTACAggtgggt                                                                                                                                                    <  2:85685/150‑1 (MQ=255)                cTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACg                                                                                                                                     <  1:464937/150‑1 (MQ=255)                              gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                                                            >  2:41378/1‑145 (MQ=255)                              gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                                                            <  1:41378/145‑1 (MQ=255)                                       ccTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                                                                    >  1:111298/1‑144 (MQ=255)                                       ccTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                                                                    <  2:111298/144‑1 (MQ=255)                                             gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc                                                                                                        >  2:330951/1‑150 (MQ=255)                                               tGGTTGGCGGGGGGCTGGAGATTAACATCGGGGTGCAGAAGCTGTTTTGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                                                      <  2:387123/150‑1 (MQ=38)                                                                gCGATTCCCCTGGGGGTGCAGAAGCTGTTGTGCTGATAGGAGAGCGAGTCGGCTGTTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCTCAAGACGGCAAGGTCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                                                                     <  2:163103/150‑1 (MQ=255)                                                                   aTTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                                               >  2:161134/1‑121 (MQ=255)                                                                   aTTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                                               <  1:161134/121‑1 (MQ=255)                                                                    ttACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCa                                                                                 <  2:508081/150‑1 (MQ=255)                                                                                    gAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccc                                                                 >  1:158675/1‑150 (MQ=255)                                                                                            tGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGtt                                                         <  2:391819/150‑1 (MQ=255)                                                                                              gtcTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCAGCTCATACCGGCGGGTGCTGCCCTGAAGACGTCAAGGCCACTGTCGAAACGGATTCATTCTTTAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTga                                                       >  1:391819/3‑149 (MQ=255)                                                                                                     ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                                                                    >  2:191310/3‑82 (MQ=17)                                                                                                     ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                                                                    <  1:191310/80‑1 (MQ=17)                                                                                                     ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                                                                    <  1:191311/80‑1 (MQ=18)                                                                                                     ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                                                                    >  2:191311/3‑82 (MQ=17)                                                                                                           ccGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtct                                          <  2:405229/150‑1 (MQ=255)                                                                                                              aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTg                                                        >  2:76112/1‑133 (MQ=255)                                                                                                              aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTg                                                        <  1:76112/133‑1 (MQ=255)                                                                                                                        ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                                                                     >  1:502473/1‑94 (MQ=255)                                                                                                                        ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                                                                     <  2:502473/94‑1 (MQ=255)                                                                                                                                     ccGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaa                     >  2:54721/1‑145 (MQ=255)                                                                                                                                     ccGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaa                     <  1:54721/145‑1 (MQ=255)                                                                                                                                           cTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg                                  >  1:491035/1‑126 (MQ=255)                                                                                                                                           cTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg                                  <  2:491035/126‑1 (MQ=255)                                                                                                                                                 gtgggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCTcac    >  1:21518/1‑150 (MQ=255)                                                                                                                                                 gtgggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGcac    >  1:21566/1‑150 (MQ=255)                                                                                                                                                   gggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAg  >  2:376482/1‑150 (MQ=255)                                                                                                                                                   gggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTCCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAg  <  1:239984/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAG  >  NC_002947/4158998‑4159293

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