Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,750,733	C→G	9.7%	A31G (GCC→GGC)	PP_1551 →	replication protein O

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,750,733	0	C	G	9.7%	147.5 / 3.6	49	A31G (GCC→GGC)	PP_1551	replication protein O
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base G (2/3);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGC  >  NC_002947/1750585‑1750876                                                                                                                                                     |                                                                                                                                                cGAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggc                                                                                                                                                 >  5:269950/1‑149 (MQ=255)   aaaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggccc                                                                                                                                               >  8:13438/1‑149 (MQ=255)      aaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCggg                                                                                                                                            <  7:13438/149‑1 (MQ=255)          aCCGCCTGGCAGTGCTGCTTTATCTACAGCATTAGAACTTGTTAAAGCATGACAAACATCGTCCAACTTGACAAGTCCAGTGGGTTCACCCGGATGGACAACCAGCTCATGTATTGCCTGCTGGCTATCGATCTCCCGTGCCGGGAGAt                                                                                                                                        <  2:285613/149‑1 (MQ=255)               cTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATg                                                                                                                                       >  1:11734/1‑145 (MQ=255)                  gCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAtt                                                                                                                                 <  5:289540/148‑1 (MQ=255)                  gCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATtg                                                                                                                                >  3:114655/1‑149 (MQ=255)                     gggCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGc                                                                                                                             <  1:42590/149‑1 (MQ=255)                        cggcTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGGCCGGGAGAtcgat                                                                                                                                    <  6:231980/139‑5 (MQ=255)                        cggcTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGGCCGGGAGAtcgat                                                                                                                                    >  5:231980/1‑135 (MQ=255)                         ggcTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTa                                                                                                                         <  5:291130/149‑1 (MQ=255)                                tctACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGc                                                                                                                   >  5:51140/1‑148 (MQ=255)                                tctACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGc                                                                                                                   >  5:232212/1‑148 (MQ=255)                                 ctACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGcc                                                                                                                  >  3:303711/1‑148 (MQ=255)                                  tACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATg                                                                                                                                       >  2:153083/1‑126 (MQ=255)                                  tACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATg                                                                                                                                       <  1:153083/126‑1 (MQ=255)                                    cAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAg                                                                                                               <  6:269950/148‑1 (MQ=255)                                      gCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCA‑GGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGcca                                                                                                           <  5:318927/149‑1 (MQ=255)                                       cATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGcca                                                                                                           >  8:109495/1‑149 (MQ=255)                                       cATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGcca                                                                                                           <  4:155987/149‑1 (MQ=255)                                        aTTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGccac                                                                                                          <  4:188359/149‑1 (MQ=255)                                            gAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATc                                                                                                      <  1:115342/149‑1 (MQ=255)                                                  gTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTc                                                                                                <  2:11734/149‑1 (MQ=255)                                                       aGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGtag                                                                                           >  8:282332/1‑147 (MQ=255)                                                                 aCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCaa                                                                                                    >  4:225636/1‑130 (MQ=255)                                                                 aCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCaa                                                                                                    <  3:225636/130‑1 (MQ=255)                                                                 aCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCaa                                                                                 <  7:277303/149‑1 (MQ=255)                                                                  cATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                <  6:51140/149‑1 (MQ=255)                                                                       tCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAAATTCGGTGCGGGCGCTCAACGCATc                                                                           >  7:287225/1‑149 (MQ=255)                                                                         ccACTTGACAAGTGCAGTGGGTTCACCAGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCgg                                                                                              <  8:292840/128‑1 (MQ=255)                                                                         ccACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCGAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATccc                                                                         >  5:280584/1‑149 (MQ=255)                                                                         ccACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCgg                                                                                              >  7:292840/1‑128 (MQ=255)                                                                                    gTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGAc                                                               >  3:112086/1‑148 (MQ=255)                                                                                            ggTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGaaa                                                      <  7:251413/149‑1 (MQ=255)                                                                                                cACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGcc                                                   <  3:164230/148‑1 (MQ=255)                                                                                                 aCCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCAt                                                 <  6:280584/149‑1 (MQ=255)                                                                                                  cccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATc                                                >  6:48091/1‑149 (MQ=255)                                                                                                   ccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATc                                                >  4:103761/1‑148 (MQ=255)                                                                                                         ggACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTc                                         >  7:109938/1‑149 (MQ=255)                                                                                                             aaCCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAccc                                     >  4:232683/1‑149 (MQ=255)                                                                                                              aCCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCt                                    >  8:17957/1‑149 (MQ=255)                                                                                                                       tggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtg                           >  7:195802/1‑149 (MQ=255)                                                                                                                        ggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtgt                          >  4:142845/1‑149 (MQ=255)                                                                                                                         gatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTc                         <  7:109495/149‑1 (MQ=255)                                                                                                                                  ctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCgg                 <  8:39475/148‑1 (MQ=255)                                                                                                                                   tgGCTATCGATCTCCCGGCCAGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTAAACGCATCCCGGCTACCGACATAGCGAAAAACATCCACGCTCACCCTGACACCGTGTCCAAGGCGGtt               >  2:138779/1‑149 (MQ=255)                                                                                                                                     gCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGttt              <  4:112086/148‑1 (MQ=255)                                                                                                                                        aTCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAg          >  6:111632/1‑149 (MQ=255)                                                                                                                                            aTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGcc        >  2:56805/1‑147 (MQ=255)                                                                                                                                            aTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGcc        <  1:56805/147‑1 (MQ=255)                                                                                                                                             tctcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgc     >  6:74936/1‑149 (MQ=255)                                                                                                                                                cccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACAcc                              >  6:253567/1‑121 (MQ=255)                                                                                                                                                cccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCtgctgc  >  7:69573/1‑149 (MQ=255)                                                                                                                                                 ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgctg   <  3:142845/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CGAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGC  >  NC_002947/1750585‑1750876

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