Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,191,856	C→T	14.6%	P908L (CCG→CTG)	PP_1937 →	helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,191,856	0	C	T	14.6%	94.7 / 7.9	41	P908L (CCG→CTG)	PP_1937	helicase
Reads supporting (aligned to +/- strand):  ref base C (17/18);  new base T (3/3);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-02

TGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCCCAGG  >  NC_002947/2191709‑2192004                                                                                                                                                    |                                                                                                                                                      tgtATGAGCTTTATCGAAACGACACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACATTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCg                                                                                                                                                      >  7:55330/1‑149 (MQ=255)     tGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCAGGCTCTGATGACCGGCTGGtca                                                                                                                                                  >  6:209784/1‑149 (MQ=255)     tGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGACGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGtca                                                                                                                                                  >  6:336853/1‑149 (MQ=255)        gCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCGGCTGGtcatca                                                                                                                                               <  5:209784/149‑1 (MQ=255)            tATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCGGCTGGtcatcaga                                                                                                                                             <  1:223930/147‑3 (MQ=255)                    cgaaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCcgc                                                                                                                                       >  1:204866/1‑145 (MQ=255)                       aacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTg                                                                                                                                <  2:66927/149‑1 (MQ=255)                        acTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTTGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCGCGGAGATTga                                                                                                                               <  5:202036/149‑2 (MQ=255)                          tAACAAGGTCATGACGCGTGTGAGCGTTGTCCCTTCCAGTGAACCATGGGTGTTACTCTGCTAGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGaaa                                                                                                                             >  2:82680/1‑149 (MQ=255)                            aCAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGAtt                                                                                                                                 >  6:304349/1‑143 (MQ=255)                            aCAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGAtt                                                                                                                                 <  5:304349/143‑1 (MQ=255)                            aCAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCaa                                                                                                                             <  2:70965/147‑1 (MQ=255)                              aaGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACg                                                                                                                           >  1:213016/1‑147 (MQ=255)                                 gTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCt                                                                                                                      <  2:204866/149‑1 (MQ=255)                                    aTGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAg                                                                                                                   >  5:147718/1‑149 (MQ=255)                                          cgTGTGAGCGTTGTCCGTTCCATTGAACCTTGGTTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTTATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACg                                                                                                             <  4:241960/149‑1 (MQ=255)                                                         ccTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGAttt                                                                                              >  3:81487/1‑149 (MQ=255)                                                           ttCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTg                                                                                             >  5:48298/1‑148 (MQ=255)                                                                                    tGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑Gc                                                                     <  4:81487/146‑1 (MQ=255)                                                                                    tGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTg                                                                   >  4:289537/1‑148 (MQ=255)                                                                                          gCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATgg                                                             <  2:213016/148‑1 (MQ=255)                                                                                           ccGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGc                                                           <  6:48298/149‑1 (MQ=255)                                                                                                  gaATTCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCAtgc                                                       >  5:272532/1‑146 (MQ=255)                                                                                                  gaATGCGCGCTCACCATGGTGGTTTTTCCTTTACGGCTCTTATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACTGTTATCAGACGATTTGCAGAAGCTCAACCCGACAAGG‑GCTGCAATGGGCAtgc                                                       <  6:272532/146‑1 (MQ=255)                                                                                                        gcgcgcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATgg                                                             >  6:83027/1‑134 (MQ=255)                                                                                                        gcgcgcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATgg                                                             <  5:83027/134‑1 (MQ=255)                                                                                                               ccAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCAtgctg                                                     >  2:316473/1‑135 (MQ=255)                                                                                                               ccAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCAtgctg                                                     <  1:316473/135‑1 (MQ=255)                                                                                                                         ttttCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGc                             <  6:147718/149‑1 (MQ=255)                                                                                                                                 tCCGGCTCTGATGACCCGCCGGGCATCAATCTCCGCGCGATTGCAACGCAGCTGAGGGGACGGGTATCAGACGTTTTGCAGAACCTCAACCCGACAAGGTGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGAt                        <  1:82680/147‑1 (MQ=255)                                                                                                                                   cGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGgaga                   <  3:289537/149‑1 (MQ=255)                                                                                                                                    gagattgatgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCaa                                                  <  8:157502/112‑1 (MQ=255)                                                                                                                                    gagattgatgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCaa                                                  >  7:157502/6‑117 (MQ=255)                                                                                                                                         tgatgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGc             >  3:44857/1‑149 (MQ=255)                                                                                                                                          gatgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGc             >  7:41512/1‑148 (MQ=255)                                                                                                                                          gatgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTg                                                                   >  6:24698/1‑94 (MQ=255)                                                                                                                                          gatgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTg                                                                   <  5:24698/94‑1 (MQ=255)                                                                                                                                           atgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCa                                     >  3:108290/1‑123 (MQ=255)                                                                                                                                           atgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCa                                     <  4:108290/123‑1 (MQ=255)                                                                                                                                                cAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCACACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGcc      <  8:55330/149‑1 (MQ=255)                                                                                                                                                  gCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGGACGGTTAACACACGATTTGCAGAACCTCAACCCGACAAGG‑GCAGCAATGGGCATGCTGCAATCGTTTCAGGCCACGCTGGGCGACATGGAGAAAATGATACTccc      >  5:310425/1‑147 (MQ=255)                                                                                                                                                    cGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCCCAgg  <  4:44857/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCCCAGG  >  NC_002947/2191709‑2192004

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