Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,763,547	G→C	16.2%	intergenic (+24/+1651)	PP_3326 → / ← PP_3328	hypothetical protein/ring‑cleaving dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,763,547	0	G	C	16.2%	72.6 / 8.6	31	intergenic (+24/+1651)	PP_3326/PP_3328	hypothetical protein/ring‑cleaving dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (15/11);  new base C (3/2);  total (18/13)
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.20e-01

TACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGC  >  NC_002947/3763399‑3763686                                                                                                                                                     |                                                                                                                                            tACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGa                                                                                                                                            >  1:688546/1‑150 (MQ=255)                tCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGa                                                                                                                                     <  1:301982/142‑1 (MQ=255)                tCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGa                                                                                                                                     >  2:301982/1‑142 (MQ=255)                        ctcGCGGCACGGTCTACCTTTTCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                                     <  2:688546/149‑1 (MQ=255)                        atcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                                     <  2:695902/150‑1 (MQ=255)                              gCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTc                                                                                                               <  2:589658/150‑1 (MQ=255)                                           tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                                  <  1:200664/150‑1 (MQ=255)                                               aaTGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTACATCTCGGTCCTGGGGAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGGCATGATGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCaa                                                                                              >  1:662659/1‑150 (MQ=255)                                                               cccGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCAt                                                                              <  2:734513/150‑1 (MQ=255)                                                                                         ttcaCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCgg                                                    >  1:685921/1‑150 (MQ=255)                                                                                            aCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCca                                                 >  2:40317/1‑150 (MQ=255)                                                                                                tCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCg                                                                                                              <  2:540000/85‑1 (MQ=255)                                                                                                tCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCg                                                                                                              >  1:540000/1‑85 (MQ=255)                                                                                                  gTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGc                                           <  2:685921/150‑1 (MQ=255)                                                                                                   ttCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCa                                          >  1:749867/1‑150 (MQ=255)                                                                                                      ctagctAAGCAGCACCAGCAGTAAACTAAGCTACATGCCTCGTCCATGACGGCCGAGGTATGGAAGCGCGGGTCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATc                                       <  1:53530/145‑1 (MQ=255)                                                                                                       tGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCg                                      >  1:499819/1‑150 (MQ=255)                                                                                                         ggTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTg                                                                                                 >  2:629872/1‑89 (MQ=255)                                                                                                         ggTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTg                                                                                                 <  1:629872/89‑1 (MQ=255)                                                                                                         ggTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGaa                                    <  1:335975/150‑1 (MQ=255)                                                                                                               ccGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGAt                              >  1:768654/1‑150 (MQ=255)                                                                                                                     cacGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCCCACGGCAATCGAACACGATGAGCAt                        >  2:404569/1‑150 (MQ=255)                                                                                                                     cacGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCAt                        >  2:277255/1‑150 (MQ=255)                                                                                                                          gTAAACTAAGCGACATACCTCGGCCATGACGGCCGCGGTATGGAAGTTGGAGGCGTTCGGAATGATGGCCTGCAACCAGCATACGCACCATTACCATTGGCACAAGCCCGGCCGCGGCCACCCGGCACTCGAACACGATCGGCATGTCGa                   >  2:114299/1‑150 (MQ=255)                                                                                                                              aCTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTAc               >  2:159547/1‑150 (MQ=255)                                                                                                                                tAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCt             >  1:550764/1‑150 (MQ=255)                                                                                                                                    tgccATACCTCGGCCGTCATGGCCGAGGTATGTCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGccc                                                            >  1:147589/4‑99 (MQ=37)                                                                                                                                    tgccATACCTCGGCCGTCATGGCCGAGGTATGTCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGccc                                                            <  2:147589/96‑1 (MQ=37)                                                                                                                                       cATACCTCGGCCGTCATGGCCGAGGTATGTCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCatgat      >  1:569294/1‑150 (MQ=255)                                                                                                                                        aTACCTCGGCCGTCATGGCCGAGGTATGTCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCatgatg     <  2:569294/150‑1 (MQ=255)                                                                                                                                           ccTCGGCCGTCATGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACGAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGAATGATgggg  >  2:271597/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                            TACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGC  >  NC_002947/3763399‑3763686

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