Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,966,349	T→A	8.9%	I2N (ATT→AAT)	PP_1762 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,966,349	0	T	A	8.9%	126.2 / 3.1	45	I2N (ATT→AAT)	PP_1762	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base A (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ACGCCCAGGGCTACCTGATCGAACGCCCCCAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGAGTTTTATGCACGACTACGCGATGGCCAGCTGAA  >  NC_002947/1966202‑1966493                                                                                                                                                    |                                                                                                                                                 aCGCCCAGGGCTACCTGATCGAACGCCCCCAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGAtt                                                                                                                                                 >  7:72833/1‑149 (MQ=255) aCGCCCAGGGCTACCTGATCGAACGCCCCCAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGAt                                                                                                                                                  <  7:308857/148‑1 (MQ=255)        gggCGACCTGATCGAACGCCCCCAGGTGTTCTCCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCAt                                                                                                                                          <  8:279511/149‑1 (MQ=255)                      aaCGCCCCCAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCgg                                                                                                                             >  5:347592/1‑148 (MQ=255)                      aaCGCCCCCAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCgg                                                                                                                             <  6:347592/148‑1 (MQ=255)                         gCCCCCAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCt                                                                                                                         >  3:72581/1‑149 (MQ=255)                            cccAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTc                                                                                                                                    <  8:276000/135‑1 (MQ=255)                            cccAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTc                                                                                                                                    >  7:276000/1‑135 (MQ=255)                             ccAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGAGCGCCGCCCGGCCGCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGa                                                                                                                     <  6:357312/149‑1 (MQ=255)                               aGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATg                                                                                                                   <  6:225233/149‑1 (MQ=255)                                 gtgtTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATg                                                                                                                   <  1:231782/147‑1 (MQ=255)                                 gtgtTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATg                                                                                                                   >  2:231782/1‑147 (MQ=255)                                       aCCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGc                                                                                                           >  2:355291/1‑149 (MQ=255)                                       aCCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGc                                                                                                           >  2:124296/1‑149 (MQ=255)                                         cTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCa                                                                                                         >  3:65817/1‑149 (MQ=255)                                          tgtgACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAg                                                                                                        >  5:335873/1‑149 (MQ=255)                                                   cTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTAccccc                                                                                                <  2:353382/148‑1 (MQ=255)                                                   cTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTAccccc                                                                                                >  1:353382/1‑148 (MQ=255)                                                    tGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGAATGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTAccccaat                                                                                              >  1:237461/1‑146 (MQ=255)                                                    tGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGAATGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTAccccaat                                                                                              >  1:232351/1‑146 (MQ=255)                                                        gcgcCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATggg                                                                                          <  7:147373/149‑1 (MQ=255)                                                         cgcCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGc                                                                                         >  3:307321/1‑149 (MQ=255)                                                          gcCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGcc                                                                                        >  7:16395/1‑149 (MQ=255)                                                                        gccCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTg                                                                           >  3:326767/1‑148 (MQ=255)                                                                        gccCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGAATGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTA‑CCCCAATGGGCCCAGCAACTGATTGaa                                                                         <  2:232351/149‑1 (MQ=255)                                                                        gccCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGAATGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTA‑CCCCAATGGGCCCAGCAACTGATTGaa                                                                         <  2:237461/149‑1 (MQ=255)                                                                                    ttGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCcgc                                                                >  5:17776/1‑147 (MQ=255)                                                                                       cAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCgaga                                                            <  4:72581/148‑1 (MQ=255)                                                                                              cACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAgcgc                                                    <  4:326767/149‑1 (MQ=255)                                                                                               aCCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGa                                                                          <  3:12688/126‑1 (MQ=255)                                                                                               aCCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGa                                                                          >  4:12688/1‑126 (MQ=255)                                                                                                      acaTTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGggtggt                                            >  4:80941/1‑149 (MQ=255)                                                                                                      acaTTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGggtggt                                            <  6:17776/149‑1 (MQ=255)                                                                                                        aTTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCga                                                              <  1:351542/129‑1 (MQ=255)                                                                                                        aTTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCga                                                              >  2:351542/1‑129 (MQ=255)                                                                                                                 gTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGa                                   <  8:72833/147‑1 (MQ=255)                                                                                                                 gTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGAg                                  <  3:80941/148‑1 (MQ=255)                                                                                                                       aaTCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCgaga                                                            <  7:253263/116‑1 (MQ=255)                                                                                                                       aaTCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCgaga                                                            >  8:253263/1‑116 (MQ=255)                                                                                                                         tCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGAGTTTTATGCa                         <  1:355291/149‑1 (MQ=255)                                                                                                                             aaGGATCAAGGAGCTGAAAGTGATCGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGAGTTTTATGCACGAc                     >  6:24323/1‑149 (MQ=255)                                                                                                                              tggATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCTTGCCTTTGATGGACAACGCCATTTACCCACAAATCGCCCAGCAACTGATAGAAGACTGCCGCGAGAGCAAGCGCCGGGTTGTTGAGCATGAGTTTTATGCAAGAc                     >  4:234094/2‑148 (MQ=255)                                                                                                                                      ggAGCGGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGAGTTTTATGCACGACTACGCGATg            <  6:335873/149‑1 (MQ=255)                                                                                                                                             aaaGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGAGTTTTATGCACGACTACGCGATGGCCAg       <  4:307321/147‑1 (MQ=255)                                                                                                                                             aaaGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGAGTTTTATGCACGACTACGCGATGGCCAGCt     <  1:124296/149‑1 (MQ=255)                                                                                                                                                gggATTGACGCATTAGATCGTATAGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGCGCCCAGCAACTGATTGAAGACTCCCGCGACAGCAAGCGCCGGGTGGTTGAGCATGAGTTTTATGCACGACTACGCGATGGCCAGCTGaa  <  4:65817/147‑1 (MQ=255)                                                                                                                                                gTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCGAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGAGTTTTATGCACGACTACGCGATGGCCAGCTg    >  5:46442/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                 ACGCCCAGGGCTACCTGATCGAACGCCCCCAGGTGTTCACCTGTGACAGTCTGCAGCGCCAACGGATCGCCGCCCGGCCTCTGTTGCAGCGAGCACCTGGCACATTTCGCTAGTTGGAAATCACAAGGATCAAGGAGCTGAAAGTGATTGACGCATTTGTTCGTATCGGGCCTTTGATGGACCCCGCCAGTTACCCCCAATGGGCCCAGCAACTGATTGAAGACTGCCGCGAGAGCAAGCGCCGGGTGGTTGAGCATGAGTTTTATGCACGACTACGCGATGGCCAGCTGAA  >  NC_002947/1966202‑1966493

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