Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,056,572	(G)6→7	coding (1180/1320 nt)	PP_3573 →	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,056,566	1	.	G	97.4%	137.3 / ‑3.0	39	coding (1174/1320 nt)	PP_3573	monooxygenase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (14/24);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

TGGACTTCAGCTGGCTGCAGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCT‑GGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCA  >  NC_002947/4056432‑4056692                                                                                                                                        |                                                                                                                               tGGACTTCAGCTGGCTGCAGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCT‑g                                                                                                                               <  2:984184/136‑1 (MQ=255) gggACTTCAGCTGGCTGCAGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGATCCTGGCGTGTACTTCCTCGGCCTGCGGTGGTTGTCGCGCCGTGGTTCATCGTTCATCT‑g                                                                                                                               <  2:640795/135‑1 (MQ=255) gggACTTCAGCTGGCTGCAGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCT‑g                                                                                                                               <  1:990077/135‑1 (MQ=255)           ctggctgCAGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCT‑gggggggtgtg                                                                                                                     <  1:1159992/136‑6 (MQ=255)             ggctgCAGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGc                                                                                                                    <  2:318138/136‑1 (MQ=255)             ggctgCAGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGc                                                                                                                    <  2:1077435/136‑1 (MQ=255)                 gCAGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTTCTTCCTCGGCCATCCGTTGGTGTCGCGCCGTGGTTGATCGTTGCTCTGGGGGGGTGTgggcagc                                                                                                                >  1:793215/1‑131 (MQ=255)                   aGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTACAAGGCCAGCCGGGGGTGTCGCGCCCTGGTTCATCGGTAATATGGGGGGGTGTTGgacgacc                                                                                                              >  1:904463/1‑135 (MQ=255)                                   tggcAACGGCAAGCCCCAGCACCAGCGTGGGGTGGGCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCgg                                                                                              <  1:794474/133‑1 (MQ=255)                                   tGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCgg                                                                                              <  2:1326236/136‑1 (MQ=255)                                   tGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCgg                                                                                              <  2:1313329/136‑1 (MQ=255)                                   tGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCgg                                                                                              <  1:13023/136‑1 (MQ=255)                                      cAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCAAAGCACGTCGCCGGCca                                                                                           <  2:731382/136‑1 (MQ=255)                                      cAACGGCAAGCCCCAGCACCAGCGTGGTGTGACCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCca                                                                                           <  1:1040011/136‑1 (MQ=255)                                          ggCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACATACTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACCCCGCCAAGCACCTCGCCCGCcccaaa                                                                                       >  1:749387/1‑134 (MQ=255)                                             aaGCCCCAGCACCAGCGTGGGGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGATCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGcc                                                                                    <  1:459332/136‑1 (MQ=255)                                                ccccAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCcaca                                                                                 <  1:121015/136‑1 (MQ=255)                                                ccccAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCcaca                                                                                 <  1:752497/136‑1 (MQ=255)                                                ccccAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCcaca                                                                                 <  2:969596/136‑1 (MQ=255)                                                ccccAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCcaca                                                                                 <  1:114687/136‑1 (MQ=255)                                                  ccAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCcacaca                                                                               <  2:268380/136‑1 (MQ=255)                                                         cAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCAcc                                                                        <  2:904463/136‑1 (MQ=255)                                                           gCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTa                                                                      <  2:152237/136‑1 (MQ=255)                                                              tggtgTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGGTTATCGTTGATCTGGGGGGGTGTGGCACGACGCCCAGCACGTCGCCGGCCACATAGCCACACCACGCCCCTaaac                                                                   >  2:931258/1‑136 (MQ=255)                                                                         gTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGGGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTaccgcg                                                        >  1:360573/1‑136 (MQ=255)                                                                             aCCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATAGTTGATCTGGGGGGGTGTGGCACGCCGCCACGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTaccgcgaccg                                                    >  2:291076/1‑136 (MQ=255)                                                                              ccTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTTGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTaccgcgaccgc                                                   >  1:289482/1‑136 (MQ=255)                                                                                    gtgtACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTTGCACGACGCCAAGCACGTCGCCCGCCACATCGCCACACAACGCACCTACACCGCCAACCGCCACCGCCAGCaa                                             >  2:1182010/1‑136 (MQ=255)                                                                                         cTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCTAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACgcgc                                        <  2:1368314/136‑1 (MQ=255)                                                                                               cGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGa                                  >  1:1204501/1‑136 (MQ=255)                                                                                                    tGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGagc                             >  2:906075/1‑136 (MQ=255)                                                                                                     gCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGagca                            <  1:642855/136‑1 (MQ=255)                                                                                                     gCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGagca                            <  1:298714/136‑1 (MQ=255)                                                                                                     gCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGagca                            <  1:960096/136‑1 (MQ=255)                                                                                                       cGTGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGagcagc                          <  1:323615/136‑1 (MQ=255)                                                                                                         tGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACCCAACGCACCTACACCGCCTACCGCGACCCCGAGCAACGCGCAGCGGATGAGCAGCaa                        >  1:133184/1‑136 (MQ=255)                                                                                                         tGGTTGTCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCaa                        >  1:496250/1‑136 (MQ=255)                                                                                                         tGGTTGTCGCGCCGTGGTTCATCGTTCATCT‑GGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCaaa                       >  1:288637/1‑136 (MQ=255)                                                                                                               tCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGcc                  >  2:1279183/1‑136 (MQ=255)                                                                                                               tCGCGCCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGcc                  <  1:437161/136‑1 (MQ=255)                                                                                                                 gcgcCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCaa                <  1:260881/136‑1 (MQ=255)                                                                                                                   gcCGTGGTTCATCGTTCATCTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGa              >  1:454919/1‑136 (MQ=255)                                                                                                                               cgttcatcTGGGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCa  <  1:313434/136‑1 (MQ=255)                                                                                                                                        |                                                                                                                               TGGACTTCAGCTGGCTGCAGGTAGACACCTTCGATGCCAACGGCAAGCCCCAGCACCAGCGTGGTGTGGCCCGTGAACCTGGCGTGTACTTCCTCGGCCTGCCGTGGTTGTCGCGCCGTGGTTCATCGTTCATCT‑GGGGGGTGTGGCACGACGCCAAGCACGTCGCCGGCCACATCGCCACACAACGCACCTACACCGCCTACCGCGACCGCGAGCAACGCGCAGCGGATGAGCAGCAAACCGCCAAGATCAGCAATGTCA  >  NC_002947/4056432‑4056692

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