Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,859,569	A→C	16.7%	intergenic (‑119/+37)	PP_1662 ← / ← PP_1663	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,859,569	0	A	C	16.7%	61.4 / 9.1	30	intergenic (‑119/+37)	PP_1662/PP_1663	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/14);  new base C (2/3);  total (13/17)
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.08e-01

TTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGT  >  NC_002947/1859422‑1859706                                                                                                                                                    |                                                                                                                                          ttcAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACt                                                                                                                                         >  2:342657/1‑150 (MQ=255)         ggCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgggg                                                                                                                                 <  1:342657/150‑1 (MQ=255)                aCAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttct                                                                                                                          <  1:271474/150‑5 (MQ=255)                aCAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTgg                                                                                                                                       <  2:774068/137‑1 (MQ=255)                aCAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTgg                                                                                                                                       >  1:774068/1‑137 (MQ=255)                 cAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGt                                                                                                                         <  1:674857/150‑1 (MQ=255)                          gCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcaggc                                                                                                                >  1:637230/1‑150 (MQ=255)                            aCGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAg                                                                                                                           <  1:456698/137‑1 (MQ=255)                            aCGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAg                                                                                                                           >  2:456698/1‑137 (MQ=255)                                        taGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGt                                                                                                  >  2:664558/1‑150 (MQ=255)                                           agTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                                                                                               <  1:250579/150‑1 (MQ=255)                                            gTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGcc                                                                                              <  2:848353/150‑1 (MQ=255)                                            gTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGcc                                                                                              <  2:495089/150‑1 (MQ=255)                                                        aCTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAg                                                                                                   <  1:665190/133‑1 (MQ=255)                                                        aCTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAg                                                                                                   >  2:665190/1‑133 (MQ=255)                                                                    aaaCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGt                                                                      <  1:664558/150‑1 (MQ=255)                                                                    aaaCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGt                                                                      <  2:637230/150‑1 (MQ=255)                                                                             tGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttc                                                                                                                           >  2:214172/1‑85 (MQ=38)                                                                             tGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttc                                                                                                                           <  1:214172/88‑4 (MQ=38)                                                                                          aTGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATcaca                                                <  2:489635/150‑1 (MQ=255)                                                                                                cTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTcc                                          >  2:263318/1‑150 (MQ=255)                                                                                                   gAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTg                                       >  1:750129/1‑150 (MQ=255)                                                                                                        gATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATc                                  <  2:750129/150‑1 (MQ=255)                                                                                                             tgcCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGtct                             <  2:72383/150‑1 (MQ=255)                                                                                                              gcCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGtctc                            <  1:263318/150‑1 (MQ=255)                                                                                                                         cTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGGGGGGCTTCAGGGATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCTGTCTCTACCTGGCGCa                 >  1:324791/1‑150 (MQ=255)                                                                                                                          tttttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAg                >  2:417133/1‑150 (MQ=255)                                                                                                                                        gAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGt  >  1:832650/1‑150 (MQ=255)                                                                                                                                         tgaaGCCCCGCCCAGTGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACcagcttgt  <  2:578560/147‑5 (MQ=255)                                                                                                                                         tgaaGCCCCGCCCAGTACGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACcagcttgt  >  1:578560/3‑145 (MQ=255)                                                                                                                                                    |                                                                                                                                          TTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGT  >  NC_002947/1859422‑1859706

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