Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,650,990	G→T	13.8%	intergenic (‑31/+23)	PP_4113 ← / ← PP_4114	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,650,990	0	G	T	13.8%	73.5 / 3.4	29	intergenic (‑31/+23)	PP_4113/PP_4114	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base T (2/2);  total (13/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 = 8.38e-01

CTGGCTATCGCTGAACAAGTCGCCGCAACGCTGCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTGGCCTTCCATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTC  >  NC_002947/4650843‑4651139                                                                                                                                                    |                                                                                                                                                      ctgGCTATCGCTGAACAAGTCGCCGCAACGCTGCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGAc                                                                                                                                                     <  1:639303/150‑1 (MQ=255)       aTCGCTGAACAAGTCGCCGCAACGCTGCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACgg                                                                                                                                               <  2:400654/150‑1 (MQ=255)              aaCAAGTCGCCGCAACGCTGCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCt                                                                                                                                        <  2:132554/150‑1 (MQ=255)                    tCGCCGCAACGCTGCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAgg                                                                                                                                  <  2:774735/150‑1 (MQ=255)                      gccgcAACGCTGCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGtt                                                                                                                                >  1:160069/1‑150 (MQ=255)                      gccgcAACGCTGCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGtt                                                                                                                                >  2:463909/1‑150 (MQ=255)                           aaCGCTGCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCt                                                                                                                           <  1:663438/150‑1 (MQ=255)                                 gCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGcc                                                                                                                     <  1:780466/150‑1 (MQ=255)                                                        tCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACgggg                                                                                                                                             >  1:272891/1‑103 (MQ=255)                                                        tCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACgggg                                                                                                                                             <  2:272891/103‑1 (MQ=255)                                                          tgctgcAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGc                                                                                            <  1:526779/150‑1 (MQ=255)                                                          tgctgcAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGc                                                                                            >  2:526779/1‑150 (MQ=255)                                                                                   cAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGc                                                                                            <  1:320957/125‑1 (MQ=255)                                                                                   cAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGc                                                                                            >  2:320957/1‑125 (MQ=255)                                                                                        ccATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAAcgcg                                                              >  2:44053/1‑150 (MQ=255)                                                                                          aTAACAGCAAGATCAGGACACCCTTCTTGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTa                                                            >  1:796214/1‑150 (MQ=255)                                                                                               aGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGc                                                       >  1:87458/1‑150 (MQ=255)                                                                                               aGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGc                                                       <  2:160069/150‑1 (MQ=255)                                                                                                       cAGGACACCCTTCTTGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCCCCACGATTg                                               <  2:796214/150‑1 (MQ=255)                                                                                                              cccTTAATGCACTCCCCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTGGCCTTcc                                        <  2:87458/150‑1 (MQ=255)                                                                                                                       cAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTGGCCTTCCATACCAGCt                               >  1:18008/1‑150 (MQ=255)                                                                                                                                       tCCCTGTGAGTGCTACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTGGCCTTCCATAcc                                   <  1:353710/130‑1 (MQ=255)                                                                                                                                       tCCCTGTGAGTGCTACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTGGCCTTCCATAcc                                   >  2:353710/1‑130 (MQ=255)                                                                                                                                       tCCCTGTGAGTGCTACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTGGCCTTCCATAcc                                   >  2:353697/1‑130 (MQ=255)                                                                                                                                       tCCCTGTGAGTGCTACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTGGCCTTCCATAcc                                   <  1:353697/130‑1 (MQ=255)                                                                                                                                       tCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTg                                               >  1:206880/1‑118 (MQ=255)                                                                                                                                       tCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTg                                               <  2:206880/118‑1 (MQ=255)                                                                                                                                                gTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCg                                                                                      <  2:135230/70‑1 (MQ=255)                                                                                                                                                gTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCg                                                                                      >  1:135230/1‑70 (MQ=255)                                                                                                                                                    gACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTGGCCTTCCATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTc  >  1:748645/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CTGGCTATCGCTGAACAAGTCGCCGCAACGCTGCAGCTGCTTGATGCGCAGCGTGTCTGCTGCAGCCTGTGCTCCAAGACTGCAAAGCCATAACAGCAAGATCAGGACACCCTTCATGCAATCACCCCACGCGCTCCCTGTGAGTGCGACTTACGGGGCACCTCAAAGGTTCACCTGTGGCCGGCGACTGCCTGGTCGCGGGCGAGCCATTCGGCCATGGCAGTATCCTCCAACGCGTACTCGCCACGATTGGCCTTCCATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTC  >  NC_002947/4650843‑4651139

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