Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,679,880	(C)10→11	coding (99/1860 nt)	PP_4983 ←	amine oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,679,870	1	.	C	75.0%	70.4 / 21.4	28	coding (109/1860 nt)	PP_4983	amine oxidase
Reads supporting (aligned to +/- strand):  ref base . (3/4);  new base C (5/16);  total (8/20)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CGGCTGTGCCTTGGCGTTGCCGATCCACTCGCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑CCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGA  >  NC_002947/5679736‑5679999                                                                                                                                         |                                                                                                                                  cGGCTGTGCCTTGGCGTTGCCGATCCACTCGCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑c                                                                                                                                  >  2:708512/1‑136 (MQ=255)  ggCTGTGCCTTGGCGTTGCCGATCCACTCGCTGGTGGCCCTGCCCGAACGCCCGCCCATGCCCGCCGTGGCCTCCAGCAGGGTCACCTGCCAGCCGTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑cc                                                                                                                                 <  2:982670/136‑1 (MQ=255)     tgtgCCTTGGCGTTGCCGATCCACTCGCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCC‑G‑ccccc                                                                                                                              >  1:1027309/1‑136 (MQ=255)     tgtgCCTTGGCGTTGCCGATCCACTCGCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑ccccc                                                                                                                              >  1:286236/1‑136 (MQ=255)      gtgCCTTGGCGTTGCCGATCCACTCGCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑cccccc                                                                                                                             >  1:247783/1‑136 (MQ=255)            tGGCGTTGCCGATCCACTCGCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑cccccccccccc                                                                                                                       >  1:44569/1‑134 (MQ=255)                 ttGCCGATCCACTCGCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑CCCCCCCCCCACGACaa                                                                                                                  >  1:1012382/1‑136 (MQ=255)                       aaccACTCGCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCATGGCCTCCATCAGGGGCAACAGCCAGCCTTTTTTCTGCAGGTCGTAGGGCGCCGTCAGGGCTTCCA‑GCCCCCCCCCCCACGACAATCGcc                                                                                                             <  2:998402/134‑1 (MQ=255)                              tgCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCATCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCg                                                                                                      <  1:467855/135‑1 (MQ=255)                                        cAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCttt                                                                                            <  2:286236/136‑1 (MQ=255)                                               gAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑CCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCcagcg                                                                                    <  1:193231/136‑1 (MQ=255)                                                aaCGCCGGCCCATGCCCGCCTTGGCCTCCTGCAGTGGCTCCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCT‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCcagcg                                                                                    <  2:822470/136‑1 (MQ=255)                                                aaCGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCcagcg                                                                                    <  1:985185/136‑1 (MQ=255)                                                aaCGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTAAGGCCTTCCATC‑CCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCcagcg                                                                                    <  2:468406/136‑1 (MQ=255)                                                   gcccgcccATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACCCCAAACGCCGGGGGGGGGGTGTGTTTTGCcagcgccc                                                                                 >  2:787950/1‑134 (MQ=255)                                                            tGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTgg                                                                        <  1:665395/136‑1 (MQ=255)                                                              cccGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCGGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGac                                                                      <  1:708512/136‑1 (MQ=255)                                                               cccccTTGGCCTCCTGCAGGGGCACCTGCCTGCCTTTGTGCTGCCCCTCGTAGCCTGCGCTCTGGCCTGCCACC‑CCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGaca                                                                     <  2:626665/133‑1 (MQ=255)                                                                     tGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACCCCAAACGCCGGGGGGGGGTTGTCTTTTGCCAGCGCCCCGCCGCGGGACCCCCcaa                                                               >  1:768688/1‑136 (MQ=255)                                                                        cagccAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACcaatca                                                            <  1:1086347/133‑1 (MQ=255)                                                                          tCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACCACAAACGCCGTGGGGGGGTTGTCTTTTGCCCGCGCCGCGCCGCGGGGCCCACcaaaaaaa                                                          >  1:1077655/1‑131 (MQ=255)                                                                             agcagGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATAcc                                                       <  2:499956/136‑1 (MQ=255)                                                                                agGGTCACCTGCCAGCCTTTGTTCTGCAGCGCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAAc                                                    <  1:29785/136‑1 (MQ=255)                                                                                  gggcACCTGCCAGCCTTTGTTCTGCAGCGCGTAGGCCGCCGTCAGGCCCGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAAcgc                                                  <  1:613168/133‑1 (MQ=255)                                                                                       cctccaaGCCTTTGTCCGGCTGCTCTTCGGCTGCCGGCTGGCCTGCCC‑C‑CCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAg                                            <  2:1275871/130‑1 (MQ=255)                                                                                        cTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAg                                            <  1:512684/136‑1 (MQ=255)                                                                                         tGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACACTCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAgg                                           <  1:204913/136‑1 (MQ=255)                                                                                         tGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGGGGGGGGGTTGTCTTTTGCCCGGGGAGCGGCGCTGGGCACCCCCATCAAAAACAACCCGCACAgg                                           >  1:1177302/1‑136 (MQ=255)                                                                                            cAGCCTATGTTCTGGAGCTCGTAGGCTGCCGGCAGGCCTGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTTGACACACCATTCAATACCAACGCGCACAGgcgc                                        <  2:1043189/136‑1 (MQ=255)                                                                                                            gCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑CCCCCCCCCCACGACAATCGCCGTGGGCGGCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACCCACCAATCAATACCAACGCGCACAGGCGCACCCAAGCACCAGCCaa                       >  1:595610/1‑135 (MQ=255)                                                                                                                 tAGGCTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCACGACATTCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCg                   <  2:1355060/136‑1 (MQ=255)                                                                                                                     cTGCCGTCAGGCCTGCCA‑GCCCCCCCCCCCCCGCCAATCGGCGGGGGGGGGTTGGGTTTTGCCGGCGCCGGGCCGGGGGGCCCCCCATGCAATACCACCGCGCCCAGGGGCCCCCAAGCGGCCGCCATGgcggcgg               >  2:246695/1‑135 (MQ=255)                                                                                                                          gTCAGGCCTGCCA‑G‑CCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTccc         >  2:311824/1‑135 (MQ=255)                                                                                                                                  tGCCA‑GCCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCaga  <  2:872653/136‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  CGGCTGTGCCTTGGCGTTGCCGATCCACTCGCTGGTGGCCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCA‑G‑CCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGA  >  NC_002947/5679736‑5679999

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