Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,271,667	A→G	8.7%	intergenic (+51/‑49)	PP_t32 → / → PP_t33	tRNA‑Asp/tRNA‑Asp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,271,667	0	A	G	8.7%	130.6 / 3.1	46	intergenic (+51/‑49)	PP_t32/PP_t33	tRNA‑Asp/tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base G (2/2);  total (23/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 = 8.16e-01

CCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCT  >  NC_002947/2271521‑2271811                                                                                                                                                   |                                                                                                                                                  ccAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATg                                                                                                                                                 >  6:128212/1‑149 (MQ=255)  cAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACAGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGc                                                                                                                                                >  3:222915/1‑149 (MQ=255)        ggAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTg‑gg                                                                                                                                         <  5:128212/149‑1 (MQ=255)         gAACGCAGTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGTCTGGCACGGGTCTGAAGCCACAAGGCACATGCCTTG‑GGc                                                                                                                                        <  4:170866/149‑1 (MQ=255)         gAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGc                                                                                                                                        >  7:155066/1‑149 (MQ=255)         gAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCC‑CAAGGCATGTGCCTTgtggc                                                                                                                                        <  1:224205/149‑5 (MQ=21)         gAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCC‑CAAGGCATGTGCCTTgtggc                                                                                                                                        >  2:224205/1‑145 (MQ=21)               aCTGCGCAGCGGTAGTTCAGTCTGTTAGAAAACCGGCCTGTCACGACGGGGGTCGCGGGTTCGAGACCAGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCCCCGGTCTGAAGCCACAAGGCACAGGCGGTG‑GGCGGATcc                                                                                                                                  >  8:174740/1‑149 (MQ=21)                         ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAAt                                                                                                                        <  7:111968/149‑1 (MQ=255)                                 aGTCGGTTAGAATACCGGCCTGTCACGCCGGGTGTCGCGGGTTCGAGTCCAGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTATGAAGCCACAAGGCACATGCCATG‑CGCTGATCCCAGTTTCAATCGGgagc                                                                                                                 >  8:175163/1‑148 (MQ=37)                                      gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCgg                                                                                                                     >  2:86433/1‑139 (MQ=255)                                      gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTATTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCgg                                                                                                                     <  1:86433/139‑1 (MQ=255)                                          gAATACCGGCCTGTCACGCCGGGGTTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                        <  7:175163/148‑1 (MQ=255)                                                   ccTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTTAACGATACGCAg                                                                                              <  1:269233/149‑1 (MQ=255)                                                   ccTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAg                                                                                              <  7:153172/149‑1 (MQ=255)                                                       tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                            <  5:83884/147‑1 (MQ=255)                                                         aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCt                                                                                                                                              <  1:5717/96‑1 (MQ=255)                                                         aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCt                                                                                                                                              >  2:5717/1‑96 (MQ=255)                                                           gCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCTGGCGCAAGCCTGAACTATACGCAGCGGTAGtt                                                                                      <  8:155066/149‑1 (MQ=255)                                                             cGGGGGTCGCGGGTTCGAGGCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCa                                                                                    <  7:174740/149‑1 (MQ=255)                                                                  gTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTAATTAACAGATCGGTGCCTGTCACCGTTCTGAAGCCACAAGTCACATGCCTTG‑GGCTTATCCCATTTTCATTCGTGCTCATTGCTTATCTATTCGCAGCGGTAGTTCAGTCgg                                                                               <  4:222915/149‑1 (MQ=255)                                                                            cGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGtt                                                                                                                             <  8:247780/93‑1 (MQ=255)                                                                            cGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGtt                                                                                                                             >  7:247780/1‑93 (MQ=255)                                                                                     tCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                        >  3:138023/1‑105 (MQ=255)                                                                                     tCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                        <  4:138023/105‑1 (MQ=255)                                                                                       cGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACg                                                          >  5:226396/1‑149 (MQ=255)                                                                                          tGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTg                                                               >  5:108569/1‑141 (MQ=255)                                                                                          tGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTg                                                               <  6:108569/141‑1 (MQ=255)                                                                                                aTTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGc                                                                                             >  6:169280/1‑105 (MQ=255)                                                                                                aTTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGc                                                                                             <  5:169280/105‑1 (MQ=255)                                                                                                              gATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTcc                                   >  5:124086/1‑149 (MQ=255)                                                                                                                   gTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTcc                              >  3:85214/1‑149 (MQ=255)                                                                                                                          gCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTgcgc                       >  5:67868/1‑149 (MQ=255)                                                                                                                              cGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATc                   >  8:62906/1‑149 (MQ=255)                                                                                                                                    gAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGtt             >  6:93114/1‑149 (MQ=255)                                                                                                                                         cccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAg                                      <  6:103302/118‑1 (MQ=17)                                                                                                                                         cccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAg                                      >  5:103302/3‑120 (MQ=17)                                                                                                                                         cacaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGt              <  1:138324/143‑1 (MQ=255)                                                                                                                                         cacaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGt              >  2:138324/1‑143 (MQ=255)                                                                                                                                            aaGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCTTTTCAAg        <  7:208507/146‑1 (MQ=255)                                                                                                                                            aaGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCCTTTCAAg        >  8:208507/1‑146 (MQ=255)                                                                                                                                              ggCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc           >  7:198522/1‑141 (MQ=255)                                                                                                                                              ggCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc           <  8:198522/141‑1 (MQ=255)                                                                                                                                              ggCACATGCCTTG‑CGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGTGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc   <  6:226396/149‑1 (MQ=255)                                                                                                                                               gCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc   >  4:240706/1‑148 (MQ=255)                                                                                                                                               gCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCt  <  7:62906/149‑1 (MQ=255)                                                                                                                                                cacaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGa                                       >  5:169815/1‑111 (MQ=255)                                                                                                                                                cacaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGa                                       <  6:169815/111‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCT  >  NC_002947/2271521‑2271811

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