Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,325,532	G→A	15.3%	intergenic (‑140/‑35)	PP_3795 ← / → pvdA	hypothetical protein/L‑ornithine 5‑monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,325,532	0	G	A	15.3%	57.9 / 3.1	27	intergenic (‑140/‑35)	PP_3795/pvdA	hypothetical protein/L‑ornithine 5‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (10/12);  new base A (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 = 9.86e-01

GCTCATGTCGATTCCCTGCCTGAAAGATGGGTCCTTGCCGTAACCAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTGGGCTTCGGCCCCTCCAACCTGGCCTTGGCCATCGCCCTGGAAGAGCTCGCC  >  NC_002947/4325387‑4325665                                                                                                                                                  |                                                                                                                                      gCTCATGTCGATTCCCTGCCTGAAAGATGGGTCCTTGCCGTAACCAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGt                                                                                                                                   >  2:512098/1‑150 (MQ=255)   tCATGTCGATTCCCTGCCTGAAAGATGGGTCCTTGCCGTAACCAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTcc                                                                                                                                 >  2:73875/1‑150 (MQ=255)                    ctgAAAGATGGGTCCTTGCCGTAACCAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAgg                                                                                                                >  1:22828/1‑150 (MQ=255)                    ctgAAAGATGGGTCCTTGCCGTAACCAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTCGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAgg                                                                                                                <  1:131691/150‑1 (MQ=255)                        aaGATGGGTCCTTGCCGTAACCAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAAt                                                                                                            <  1:73875/150‑1 (MQ=255)                               gTCCTTGCCGTAACCAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGcc                                                                                                     >  2:245582/1‑150 (MQ=255)                                       cGTAACCAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCa                                                                                             <  2:22828/150‑1 (MQ=255)                                            ccAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCt                                                                                        >  2:501524/1‑150 (MQ=255)                                             cAACACATCGGCCGCATCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTc                                                                                       >  1:101612/1‑150 (MQ=255)                                                   aTCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCAcagca                                                                                 >  1:28183/1‑150 (MQ=255)                                                             tCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATc                                                                       >  2:60202/1‑150 (MQ=255)                                                                      tCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATg                                                              >  2:516873/1‑150 (MQ=255)                                                                             cccGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGa                                                                                                                         <  1:393393/84‑1 (MQ=255)                                                                             cccGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGa                                                                                                                         >  2:393393/1‑84 (MQ=255)                                                                                 tttAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTggg                                                   <  2:520550/149‑1 (MQ=255)                                                                                 cTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGAGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAGGACATGATCGGCGTgg                                                    >  1:69946/1‑149 (MQ=255)                                                                                 cTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGAGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTgg                                                    >  1:115731/1‑149 (MQ=255)                                                                                 cTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGAGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTgg                                                    <  2:69946/149‑1 (MQ=255)                                                                                 cTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGAGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTgg                                                    <  2:115731/149‑1 (MQ=255)                                                                                    aaaTTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGCGCGTCCATATTCGACCATTCAGGGAAGCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTGGGCtt                                                <  1:478440/150‑1 (MQ=255)                                                                                     aaTTCCCCCGAGCCTGCCTCGTTCCAACTCATGTACCTTCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTGGGCTTc                                               <  1:505023/150‑1 (MQ=255)                                                                                        tCCCCCTCGGCTGGCTCGTTCCAACTGATGTACTTGCGCCGGTGCCGGAGCCCCGGCGGGCGTCCATATTCGCCCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTGGGCTTCGGc                                            <  2:397645/150‑1 (MQ=255)                                                                                                    ggctCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTGGGCTTCGGCCCCTCCAACCTg                                <  1:271148/150‑1 (MQ=255)                                                                                                           tCCAACTCATGTACTAGCGCAGAAGCCGGATCCCCGGCTGCCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGGCAGTCCTCACAGCTGGAACCCATCAAAGACATGATCGGCGTGGGCTTCGGCCCCTCCTGCCTGGCCTTgg                         <  2:39191/150‑1 (MQ=255)                                                                                                                     gTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTGGGCTTCGGCCCCTCCAACCTGGCCTTGGCCATCGCCCt               <  1:245582/150‑1 (MQ=255)                                                                                                                         ttGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTGGGCTTCGGCCCCTCCAACCTGGCCTTGGCCATCGCCCTGGaa           <  1:501524/150‑1 (MQ=255)                                                                                                                                  agccggagccCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTGGGCTTCGGCCCCTCCAACCTGGCCTTGGCCATCGCCCTGGAAGAGCTCGcc  <  1:512098/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                      GCTCATGTCGATTCCCTGCCTGAAAGATGGGTCCTTGCCGTAACCAACACATCGGCCGCCTCACGCTGTTCTTTAGCCCGCTTAAATTCCCCCTCGGCTGGCTCGTTCCAACTCATGTACTTGCGCAGGAGCCGGAGCCCCGGCGGGCGTCCATATTCGACCATTCAGGGAATCGCAGCCATGAGCCAGTCCTCACAGCAGGAAACCATCAAAGACATGATCGGCGTGGGCTTCGGCCCCTCCAACCTGGCCTTGGCCATCGCCCTGGAAGAGCTCGCC  >  NC_002947/4325387‑4325665

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