Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,191,853	A→G	15.3%	Q907R (CAG→CGG)	PP_1937 →	helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,191,853	0	A	G	15.3%	87.0 / 9.4	40	Q907R (CAG→CGG)	PP_1937	helicase
Reads supporting (aligned to +/- strand):  ref base A (16/17);  new base G (3/3);  total (19/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 = 8.61e-01

TGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCC  >  NC_002947/2191709‑2192000                                                                                                                                                 |                                                                                                                                                     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                                                         <  5:83027/134‑1 (MQ=255)                                                                                                        gcgcgcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATgg                                                         >  6:83027/1‑134 (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                                              >  7:157502/6‑117 (MQ=255)                                                                                                                                    gagattgatgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCaa                                              <  8:157502/112‑1 (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)                                                                                                                                                 |                                                                                                                                                     TGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGG‑GCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCC  >  NC_002947/2191709‑2192000

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