Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,859,573	G→T	17.3%	intergenic (‑123/+33)	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,573	0	G	T	17.3%	69.4 / 4.2	28	intergenic (‑123/+33)	PP_1662/PP_1663	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (9/14);  new base T (2/3);  total (11/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 = 1.00e+00

GGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCCCAGTTCTT  >  NC_002947/1859430‑1859719                                                                                                                                                |                                                                                                                                                   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                                                                                   <  2:637230/150‑1 (MQ=255)                                                            aaaCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGt                                                                                   <  1:664558/150‑1 (MQ=255)                                                                     tGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttc                                                                                                                                        <  1:214172/88‑4 (MQ=38)                                                                     tGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttc                                                                                                                                        >  2:214172/1‑85 (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)                                                                                                                                             cTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCCCAGTTCtt  <  2:832650/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCCCAGTTCTT  >  NC_002947/1859430‑1859719

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