Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,751,425	C→G	11.2%	H262D (CAT→GAT)	PP_1551 →	replication protein O

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,751,425	0	C	G	11.2%	139.6 / 7.9	52	H262D (CAT→GAT)	PP_1551	replication protein O
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base G (3/3);  total (26/26)
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

CAGTGCCTGACCTTCGGTATCCAGCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGTTTCCGCGTCCAGCCTGCGCACTGCGAGC  >  NC_002947/1751287‑1751564                                                                                                                                            |                                                                                                                                            cAGTGCCTGACCTTCGGTATCCAGCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTGATGACGGCtg                                                                                                                                   >  4:76396/1‑148 (MQ=255)    tGCCTGACCTTCGGTATCCAGCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTGATGACGGCtgttc                                                                                                                                <  3:76396/149‑5 (MQ=255)          aCCTTCGGTATCCAGCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACcg                                                                                                                          >  3:8339/1‑149 (MQ=255)            cTTCGGTATCCAGCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTGATGACGGCtgtt                                                                                                                                 <  8:43443/140‑4 (MQ=255)            cTTCGGTATCCAGCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTGATGACGGCtgtt                                                                                                                                 >  7:43443/1‑137 (MQ=255)                  tATCCAGCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTaca                                                                                                                  >  3:81649/1‑149 (MQ=255)                      cAGCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACAccgc                                                                                                              >  6:294198/1‑149 (MQ=255)                        gCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAAc                                                                                                                               >  4:70827/1‑130 (MQ=255)                        gCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAAc                                                                                                                               <  3:70827/130‑1 (MQ=255)                          cTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGTCTTAATCACGGCTTCAACGACCTCGACTACACCGCCGGc                                                                                                          <  4:186469/149‑1 (MQ=255)                             cccAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCTAATATATTTAATCACTGCTTGAACTACGTCTACTACACCGCCGGCCt                                                                                                        <  4:53110/148‑1 (MQ=255)                                   tCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAg                                                                                                 >  8:278086/1‑149 (MQ=255)                                   tCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAg                                                                                                 >  8:230519/1‑149 (MQ=255)                                    cATGGAAATCGTCATCGAGAGCGGGTGGGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAg                                                                                                 >  1:166446/1‑149 (MQ=255)                                           aTCGTGATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCTAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGa                                                                                         <  3:19092/149‑1 (MQ=255)                                           aTCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGa                                                                                         >  1:185711/1‑149 (MQ=255)                                                cATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCa                                                                                    >  5:174282/1‑149 (MQ=255)                                                      gagCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATg                                                                              >  8:57977/1‑149 (MQ=255)                                                        gCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCg                                                                            <  2:185711/149‑1 (MQ=255)                                                                gcgcAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGTCTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATc                                                                         <  1:70710/145‑1 (MQ=255)                                                                gcgcAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATc                                                                         >  2:70710/1‑145 (MQ=255)                                                                 cgcAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAAt                                                                    >  6:92282/1‑149 (MQ=255)                                                                 cgcAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAAt                                                                    <  5:294198/149‑1 (MQ=255)                                                                   ctgCTTCGTGGTGGAATGGTTCACCTAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGTTCAGCCGTCATCTCGGCTTCAACGACCGCGACTACAGCGCCGGCCTGTCCGAGCGTGAGGACGGCACGTATGCGATCTGAATc                                                                   <  5:200824/146‑1 (MQ=255)                                                                      cTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCgg                                                                                                           <  4:11463/105‑1 (MQ=255)                                                                      cTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCgg                                                                                                           >  3:11463/1‑105 (MQ=255)                                                                      cTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAACCCGAACAGCCGTCATCACGGCTTCAACGACCGCCACTACACCGCCGGCCTCTCCGAGCGAGAGGACGCCACCTATGCGATCTGAATCgag                                                                >  6:200824/1‑146 (MQ=255)                                                                            ggTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACt                                                         <  7:57977/149‑1 (MQ=255)                                                                                gAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGt                                                     <  7:278086/149‑1 (MQ=255)                                                                                    ggATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGa                                                              <  3:166294/136‑1 (MQ=255)                                                                                    ggATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGa                                                              >  4:166294/1‑136 (MQ=255)                                                                                         aCCAAGCGATTCGCCGCCAAATCGCCCGCACAGTGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTg                                             <  5:136326/148‑1 (MQ=255)                                                                                         aCCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAg                                                                                           <  3:210904/102‑1 (MQ=255)                                                                                         aCCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAg                                                                                           >  4:210904/1‑102 (MQ=255)                                                                                            aaGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGa                                          >  5:263682/1‑148 (MQ=255)                                                                                            aaGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGa                                          >  4:187047/1‑148 (MQ=255)                                                                                               cGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAAcg                                      <  6:174282/149‑1 (MQ=255)                                                                                                 aTTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAAcgcc                                    <  4:8339/149‑1 (MQ=255)                                                                                                     gccgccAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCCCCTATGCGATCTGAAGCGGTGATTACTATGTCCGAc                                                >  8:202681/1‑133 (MQ=255)                                                                                                     gccgccAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGt                                                     >  4:146277/1‑128 (MQ=255)                                                                                                     gccgccAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGt                                                     <  3:146277/128‑1 (MQ=255)                                                                                                     gccgccAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACgccgcc                                 <  5:92282/148‑1 (MQ=255)                                                                                                        gccAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGtt                             <  7:230519/149‑1 (MQ=255)                                                                                                                  ccGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGTTTCCGCGTCCa                   <  2:166446/149‑1 (MQ=255)                                                                                                                      acaGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGTTTCCGCGTCCAGCCt               >  1:10291/1‑149 (MQ=255)                                                                                                                            cAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGccgccg                                <  7:124743/126‑1 (MQ=255)                                                                                                                            cAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGccgccg                                >  8:124743/1‑126 (MQ=255)                                                                                                                               gCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGTTTCCGCGTCCAGCCTGCGCACTGc      <  3:187047/149‑1 (MQ=255)                                                                                                                                  ttgaaGCCGTGATGACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGTTTCcgc                        >  7:157865/5‑128 (MQ=255)                                                                                                                                  ttgaaGCCGTGATGACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGTTTCcgc                        <  8:157865/124‑1 (MQ=255)                                                                                                                                  gAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGTTTCCGCGTCCAGCCTGCGCACTGCGAg   >  3:81217/1‑149 (MQ=255)                                                                                                                                    aCAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGTTTCCGCGTCCAGCCTGCGCACTGCGAGc  <  2:10291/148‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            CAGTGCCTGACCTTCGGTATCCAGCCTGCCCAGGTCATGGAAATCGTCATCGAGAGCGGGT‑GGCGCAGCTTCGAGGTGGAATGGATCACCAAGCGATTCGCCGCCAAATCGCCCGCACAGGGCAAGCCGAACAGCCGTCATCACGGCTTCAACGACCGCGACTACACCGCCGGCCTGTCCGAGCGTGAGGACGGCACCTATGCGATCTGAATCGGTGATTACTATGTCCGACGTGCGAAACGCCGCCGGTTTCCGCGTCCAGCCTGCGCACTGCGAGC  >  NC_002947/1751287‑1751564

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