Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,311,503	Δ1 bp	8.9%	coding (223/1152 nt)	PP_3784 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,311,503	0	A	.	8.9%	136.2 / 7.6	45	coding (223/1152 nt)	PP_3784	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base . (2/2);  total (23/22)
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.07e-01

CCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCGTGATCGGCTCGCCCGCGGCGTTACAATGGC  >  NC_002947/4311357‑4311650                                                                                                                                                   |                                                                                                                                                    ccACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAg                                                                                                                                                   >  2:126912/1‑149 (MQ=255)            ccTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGc                                                                                                                                        <  3:239317/149‑1 (MQ=255)                                ccTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAgg                                                                                                                    >  8:82250/1‑149 (MQ=255)                                 cTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAgg                                                                                                                    >  5:87974/1‑148 (MQ=255)                                         aTCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTc                                                                                                           <  7:202465/149‑1 (MQ=255)                                            aTTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAgc                                                                                                        >  8:87969/1‑149 (MQ=255)                                                ccgccATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAgccgc                                                                                                     <  2:265973/148‑1 (MQ=255)                                                      tCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTtctc                                                                                              <  1:126912/149‑1 (MQ=255)                                                        cAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCGACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCgg                                                                                            <  6:236941/149‑1 (MQ=255)                                                         aGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGAACTGCTTTCAAGGGCTGAGGTATTCGTTCAGCCGACTTTCTCggg                                                                                           >  7:288248/1‑149 (MQ=255)                                                              cGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAg                                                                                      >  8:325897/1‑149 (MQ=255)                                                               gAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTAGGGGCAAGt                                                                                     >  3:36125/1‑149 (MQ=255)                                                               gAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACG‑AGGATTTCCAGGCACTGCTTTCACGGGCTGAg                                                                                                                     >  7:29897/1‑116 (MQ=255)                                                               gAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACG‑AGGATTTCCAGGCACTGCTTTCACGGGCTGAg                                                                                                                     >  7:24492/1‑116 (MQ=255)                                                               gAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACG‑AGGATTTCCAGGCACTGCTTTCACGGGCTGAg                                                                                                                     <  8:29897/116‑1 (MQ=255)                                                               gAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACG‑AGGATTTCCAGGCACTGCTTTCACGGGCTGAg                                                                                                                     <  8:24492/116‑1 (MQ=255)                                                                     gCACCACTTGGCGTGCTGCTTCTTGTCTTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAACGATTACCAGGCACTGATTTCACGGGCTGAGGTAGTCGTTAAGCCGCCTTTCTAGGGGCAAGTTCaaa                                                                                >  5:220508/1‑148 (MQ=255)                                                                      caccacTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTACAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGCTCCAAGTTCAAAt                                                                               >  6:304433/1‑148 (MQ=255)                                                                        ccacTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAg                                                                            <  8:243410/149‑1 (MQ=255)                                                                                 gTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCctggct                                                                   <  7:325897/149‑1 (MQ=255)                                                                                 gTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCctggct                                                                   <  6:87974/149‑1 (MQ=255)                                                                                     tgcGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCctggctgg                                                                 >  7:230687/1‑147 (MQ=255)                                                                                       cGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTa                                                               <  8:230687/147‑1 (MQ=255)                                                                                        gTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGa                                                             >  7:273375/1‑148 (MQ=255)                                                                                           ttGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGAcaca                                                         <  4:133214/149‑1 (MQ=255)                                                                                           ttGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTGCAAATCAGTGCCTGGCTGGTAGaca                                                           >  6:156195/1‑147 (MQ=255)                                                                                             gTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGacac                                                          <  7:82250/146‑1 (MQ=255)                                                                                              tCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTtctc                                                                                              >  7:194907/1‑109 (MQ=255)                                                                                              tCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTtctc                                                                                              <  8:194907/109‑1 (MQ=255)                                                                                                gTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGc                                                    <  7:87969/149‑1 (MQ=255)                                                                                                      gTGATCAACCATAAACAGATCGGCAACGATCTGGAGCAACTGACGAAGGATTTCCAGGCAATGCTTTCACGGGCTTAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGc                                              <  8:288248/149‑1 (MQ=255)                                                                                                      gTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGc                                              >  5:187610/1‑149 (MQ=255)                                                                                                      gTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTAGTAGACACAGCTGCCATTGc                                              >  3:354738/1‑149 (MQ=255)                                                                                                                    aCAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGATCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCg                                >  3:159929/1‑149 (MQ=255)                                                                                                                             cACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGaca                                                           <  4:195161/113‑1 (MQ=255)                                                                                                                             cACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGaca                                                           >  3:195161/1‑113 (MQ=255)                                                                                                                              aCCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGATCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCGTGATCGGCTc                      <  4:159929/149‑1 (MQ=255)                                                                                                                               acGTTCTGGATCAACTTACGAAGGATTTCCAGGCACTGCTTTCTCTGGATGATGTAGTCGTTCAGCCGCATTTCTCGGTGCAAGTTAATATCATTGCCTGGCTGGTAGAAACAGATGACATTGACGATGCACAAGCCGTGATCGGCta                      <  3:82768/147‑2 (MQ=255)                                                                                                                                cGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCGTGATCGGCTCGc                    <  8:294714/149‑1 (MQ=255)                                                                                                                                  ttCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCGTGATCGGCTCGccc                  >  8:22511/1‑149 (MQ=255)                                                                                                                                         gCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCGTGATCGGCTCGCCCGCGGCGt           <  4:36125/149‑1 (MQ=255)                                                                                                                                          cAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCGTGATCGGCTCGCCCGCGGCGtt          <  6:187610/149‑1 (MQ=255)                                                                                                                                                aCGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCGTGATCGGCTCGCCCGCGGCGTTACAATg    >  6:307599/1‑149 (MQ=255)                                                                                                                                                  gAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTTCCGATGCACCAGCCGTGATCGGCTCGCCCGCGGCGTTACAATGGc  >  8:120242/1‑149 (MQ=255)                                                                                                                                                  gAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCGTGATCGGCTCGCCCGCGGCGTTACAATGGc  >  6:247681/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGCAACCAGAAACAGATCGGCACCGTTCTGGAGCAACTGACGAAGGATTTCCAGGCACTGCTTTCACGGGCTGAGGTAGTCGTTCAGCCGCCTTTCTCGGGGCAAGTTCAAATCAGTGCCTGGCTGGTAGACACAGCTGCCATTGCCGATGCACCAGCCGTGATCGGCTCGCCCGCGGCGTTACAATGGC  >  NC_002947/4311357‑4311650

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