Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,209,427	A→G	60.8%	R236R (CGT→CGC)	rnd ←	ribonuclease D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,209,427	0	A	G	60.8%	34.8 / 50.8	51	R236R (CGT→CGC)	rnd	ribonuclease D
Reads supporting (aligned to +/- strand):  ref base A (9/11);  new base G (18/13);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TGGCGGCGCGCTTGATCAGCTGGATGAGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAAC  >  NC_002947/5209279‑5209565                                                                                                                                                      |                                                                                                                                           tGGCGGCGCGCTTGATCAGCTGGATGAGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTa                                                                                                                                            >  8:221751/1‑149 (MQ=255)    cggcgCGCTTGATCAGCTGGATGAGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTgcgc                                                                                                                                         <  4:151981/149‑1 (MQ=255)          gcTTGATCAGCTGGATGAGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGcaca                                                                                                                                   <  5:285214/149‑1 (MQ=255)               aTCAGCTGGATGAGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATcgcg                                                                                                                              <  4:189054/149‑1 (MQ=255)               aTCAGCTGGATGAGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATcgcg                                                                                                                              <  6:101294/149‑1 (MQ=255)                   gCTGGATGAGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTg                                                                                                                          >  2:215242/1‑149 (MQ=255)                       gATGAGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTgcgcg                                                                                                                      >  8:187125/1‑149 (MQ=255)                           aGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCttc                                                                                                                  <  1:215242/149‑1 (MQ=255)                           aGGAACTCACCGTCCTGACGAATGGTGCGCGGGAGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGGTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCtg                                                                                                                   <  7:105419/148‑2 (MQ=255)                            ggAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCttcg                                                                                                                 >  3:327775/1‑149 (MQ=255)                             gAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCttcgc                                                                                                                <  7:221751/149‑1 (MQ=255)                               aCTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCttcgcgt                                                                                                              >  8:137662/1‑149 (MQ=255)                                 tCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCttcgcgttc                                                                                                            >  7:207091/1‑149 (MQ=255)                                 tCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCttcgcgttc                                                                                                            >  1:253653/1‑149 (MQ=255)                                   aCCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCttcgcgttcgc                                                                                                          >  1:140794/1‑149 (MQ=255)                                      gTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGcc                                                                                                       <  4:327775/149‑1 (MQ=255)                                               aaTGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGcaca                                                                                              <  2:253653/149‑1 (MQ=255)                                                       gcgGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGt                                                                                      <  7:187125/149‑1 (MQ=255)                                                           gTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCa                                                                                  >  4:246558/1‑149 (MQ=255)                                                                tctcGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGc                                                                              >  5:140879/1‑148 (MQ=255)                                                                    gtcgATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCATCGGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAgt                                                                          >  1:266599/1‑149 (MQ=255)                                                                    gtcgATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCa                                                                            <  3:246558/146‑1 (MQ=255)                                                                    gtcgATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAgtt                                                                         >  4:102650/1‑149 (MQ=255)                                                                             ggCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGa                                                                    <  3:154440/145‑1 (MQ=255)                                                                                cAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGgcgc                                                               >  5:205839/1‑147 (MQ=255)                                                                                            gttgttCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGc                                                   >  5:234046/1‑147 (MQ=255)                                                                                            gttgttCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGc                                                   >  3:105302/1‑147 (MQ=255)                                                                                              tgttCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGc                                               >  3:353147/1‑149 (MQ=255)                                                                                                ttCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCtt                                             <  2:266599/149‑1 (MQ=255)                                                                                                  cGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGa                                           >  5:134113/1‑149 (MQ=255)                                                                                                         ctTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTcgcg                                    <  6:134113/149‑1 (MQ=255)                                                                                                          ttttGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCgg                                   >  3:179790/1‑149 (MQ=255)                                                                                                            ttGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCgg                                   <  4:179790/147‑1 (MQ=255)                                                                                                               gCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAg                              <  6:140879/149‑1 (MQ=255)                                                                                                                ccAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAgg                             >  5:256634/1‑149 (MQ=255)                                                                                                                 cAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGc                            >  6:335501/1‑149 (MQ=255)                                                                                                                   t‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCtt                          >  5:300930/1‑149 (MQ=255)                                                                                                                       gCCACAGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCg                       <  4:105302/149‑1 (MQ=255)                                                                                                                        ccACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCgg                      >  4:183880/1‑149 (MQ=255)                                                                                                                            aGCGAGTGCTCCTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTc                  <  7:137662/149‑1 (MQ=255)                                                                                                                            aGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTc                  <  5:335501/149‑1 (MQ=255)                                                                                                                                    cTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCAc          >  6:137556/1‑149 (MQ=255)                                                                                                                                     tCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGgcgc                                                               <  6:299469/95‑1 (MQ=255)                                                                                                                                     tCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGgcgc                                                               >  5:299469/1‑95 (MQ=255)                                                                                                                                      ccTTGAGGATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACgg        >  5:352089/1‑149 (MQ=255)                                                                                                                                       cTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTc                                                        <  4:129111/100‑1 (MQ=255)                                                                                                                                       cTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTc                                                        >  3:129111/1‑100 (MQ=255)                                                                                                                                          gAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCt                                                          <  6:210046/95‑1 (MQ=255)                                                                                                                                          gAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCt                                                          >  5:210046/1‑95 (MQ=255)                                                                                                                                           aGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCa    <  5:137556/148‑1 (MQ=255)                                                                                                                                            ggATACGATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAAc  <  8:207091/149‑1 (MQ=255)                                                                                                                                                  gATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACa                               >  6:184137/1‑114 (MQ=255)                                                                                                                                                  gATTGCGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACa                               <  5:184137/114‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                           TGGCGGCGCGCTTGATCAGCTGGATGAGGAACTCACCGTCCTGACGAATGGTGCGCGGGTGCATCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCAT‑GGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAAC  >  NC_002947/5209279‑5209565

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