Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,681,393	A→T	8.5%	S155C (AGT→TGT)	PP_1476 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,681,393	0	A	T	8.5%	136.6 / 3.6	47	S155C (AGT→TGT)	PP_1476	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/26);  new base T (2/2);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCC  >  NC_002947/1681251‑1681541                                                                                                                                               |                                                                                                                                                     cGGCTGCGTGCGTGTGTGGGTGGGTAAAGTGGTGAGCGCGGTCAAGAAGGTGAAATAGCTGAGCATTGCGGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACTTGGCGCGTAGGATCAAGTGGGACGTGATGATCAGTGAAg                                                                                                                                                <  4:358978/149‑1 (MQ=255) cGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGaa                                                                                                                                                 <  5:317815/148‑1 (MQ=255)        cTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGt                                                                                                                                         >  8:369585/1‑149 (MQ=255)            gAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTc                                                                                                                                     <  4:365974/149‑1 (MQ=255)            gAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTc                                                                                                                                     <  4:230218/149‑1 (MQ=255)                 tgGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGtga                                                                                                                                 <  6:57646/148‑1 (MQ=255)                         tAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGAcc                                                                                                                        >  8:151974/1‑149 (MQ=255)                            aGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTg                                                                                                                     <  4:404165/149‑1 (MQ=255)                             gtggtgAGCGCGGTCAAGAACCTGAAAAACCTGACCATTGCCGCTGTCCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCTAGTGTGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTgg                                                                                                                    <  7:78111/149‑1 (MQ=255)                                gtgAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGc                                                                                                                   >  2:150298/1‑147 (MQ=255)                                       cgGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGaa                                                                                                          <  2:1303/149‑1 (MQ=255)                                        ggTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAg                                                                                                         >  3:405630/1‑149 (MQ=255)                                              gAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAAc                                                                                                   >  6:267896/1‑149 (MQ=255)                                                  cTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCa                                                                                               <  7:141897/149‑1 (MQ=255)                                                     aaaaaCCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTgg                                                                                                                    >  8:147615/1‑125 (MQ=255)                                                     aaaaaCCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTgg                                                                                                                    <  7:147615/125‑1 (MQ=255)                                                      aaaaCCTGAGCATTGCTGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTgg                                                                                           <  7:383714/149‑1 (MQ=255)                                                            tGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTc                                                                                     >  5:373533/1‑149 (MQ=255)                                                               gCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGt                                                                                  <  6:373533/149‑1 (MQ=255)                                                                 aTTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGt                                                                                  >  1:394425/1‑147 (MQ=255)                                                                  ttGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAg                                                                               >  7:178160/1‑149 (MQ=255)                                                                    gCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGc                                                                              <  4:218879/148‑1 (MQ=255)                                                                       gCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTgg                                                                          <  8:178160/149‑1 (MQ=255)                                                                           tgCCGCAAGAGGTGCTGGAAACCCTGGCGACGGACATGCCGCGTATCATCATGGGTGACTTGATGATAAGTGAATGTACGGTATTCTTGACTGACGACCTTTTCCAGCATGAAGTGCAACTTCAGTGGCTGCTCGGTGAGCgttgtt                                                                        <  8:24194/147‑1 (MQ=255)                                                                             ccGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCTTAGCATCAAGTGTGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGTCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGaa                                                                     <  5:224950/148‑1 (MQ=255)                                                                             ccGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGaa                                                                     <  7:8376/148‑1 (MQ=255)                                                                                cAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCaa                                                                 <  1:150298/149‑1 (MQ=255)                                                                                            gAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTAt                                                      <  7:369585/148‑1 (MQ=255)                                                                                                 ccTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCTGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGaa                                                                     <  3:120605/128‑1 (MQ=255)                                                                                                 ccTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCTGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGaa                                                                     <  3:126298/128‑1 (MQ=255)                                                                                                 ccTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCTGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGaa                                                                     >  4:126298/1‑128 (MQ=255)                                                                                                 ccTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCTGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGaa                                                                     >  4:120605/1‑128 (MQ=255)                                                                                                      cgaccgacATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGcc                                           <  7:151974/149‑1 (MQ=255)                                                                                                                     gTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGAccc                            >  7:65405/1‑149 (MQ=255)                                                                                                                      tAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCt                           <  7:173348/149‑1 (MQ=255)                                                                                                                       aGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTg                          >  2:107049/1‑149 (MQ=255)                                                                                                                        gCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTg                          >  7:57357/1‑148 (MQ=255)                                                                                                                           tCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGt                                                                                  >  6:154235/1‑89 (MQ=255)                                                                                                                           tCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAaca                                                               <  3:378552/108‑1 (MQ=255)                                                                                                                           tCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAaca                                                               >  4:378552/1‑108 (MQ=255)                                                                                                                                 gggACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAacac                                                              <  1:205668/103‑1 (MQ=255)                                                                                                                                 gggACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAacac                                                              >  2:205668/1‑103 (MQ=255)                                                                                                                                  ggACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTgg                <  8:65405/148‑1 (MQ=255)                                                                                                                                    aCGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCg                                         <  2:196576/121‑1 (MQ=255)                                                                                                                                    aCGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCg                                         >  1:196576/1‑121 (MQ=255)                                                                                                                                     cGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTaa            <  2:394425/149‑1 (MQ=255)                                                                                                                                        gatgatCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCt         >  8:19376/1‑149 (MQ=255)                                                                                                                                               aGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCAccc  >  5:70405/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     CGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCC  >  NC_002947/1681251‑1681541

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