Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,427,343	T→C	10.9%	N270D (AAT→GAT)	rhlE‑I ←	ATP‑dependent RNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,427,343	0	T	C	10.9%	142.9 / 8.3	55	N270D (AAT→GAT)	rhlE‑I	ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base C (4/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

TTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGCT  >  NC_002947/5427198‑5427492                                                                                                                                                   |                                                                                                                                                       ttGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACg                                                                                                                                                    >  2:307718/1‑150 (MQ=255) ttGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACg                                                                                                                                                    >  1:320349/1‑150 (MQ=255)       acCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCg                                                                                                                                             <  2:33254/150‑1 (MQ=255)              gCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCa                                                                                                                                      >  1:579162/1‑150 (MQ=255)                 aGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAgc                                                                                                                                   <  1:303718/150‑1 (MQ=255)                  gTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAgcg                                                                                                                                  >  1:352504/1‑150 (MQ=255)                     aTCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCATCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTc                                                                                                                               >  1:47021/1‑150 (MQ=255)                              cAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTg                                                                                                                      >  2:228169/1‑150 (MQ=255)                                 gCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATc                                                                                                                   >  2:579861/1‑150 (MQ=255)                                            ccACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGc                                                                                                        >  1:111808/1‑150 (MQ=255)                                              acaTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                                      <  2:579162/150‑1 (MQ=255)                                              acaTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                                      <  2:352504/150‑1 (MQ=255)                                                   ggTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCtt                                                                                                 >  2:166910/1‑150 (MQ=255)                                                     tCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTgg                                                                                               <  1:307718/150‑1 (MQ=255)                                                      cGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGc                                                                                              <  2:514085/150‑1 (MQ=255)                                                         caccaGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCa                                                                                                                                         >  2:561974/1‑104 (MQ=255)                                                         caccaGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCa                                                                                                                                         <  1:561974/104‑1 (MQ=255)                                                            caGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGATCACAGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCGGATCGCCGCCATTGCGGGTCTTGGCGAacac                                                                                        >  1:57591/1‑150 (MQ=255)                                                             aGCACTTGCACCTCACGGGCCTTGAAGCTGTACAGCGCACGCTGGCGGGTTGCCTGCGCGGCGTCACCGTGGAAGACGGACCCCAGCGTCGACCTGCGG‑‑ATCAATCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACAcc                                                                                       >  1:293591/1‑150 (MQ=255)                                                               cACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCCGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAg                                                                                     >  1:150214/1‑150 (MQ=255)                                                                    gCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTg                                                                                                <  1:532844/134‑1 (MQ=255)                                                                    gCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTg                                                                                                >  2:532844/1‑134 (MQ=255)                                                                        cTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCtt                                                                            >  1:337972/1‑150 (MQ=255)                                                                        cTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCtt                                                                            >  2:52446/1‑150 (MQ=255)                                                                         tCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTc                                                                           >  1:294979/1‑150 (MQ=255)                                                                          cACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTcc                                                                          >  1:189552/1‑150 (MQ=255)                                                                               gCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAgcgc                                                                     >  1:202924/1‑150 (MQ=255)                                                                                    gAAGCTGTCCAGCGCACGCTGGCGGGTTGACTGTGGCAGGTCACCGTGGATGCC‑GGCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCtt                                                                >  2:6450/1‑150 (MQ=255)                                                                                     aaGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTg                                                               >  1:382518/1‑150 (MQ=255)                                                                                         tGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCgca                                                           <  2:337972/150‑1 (MQ=255)                                                                                            ccAGCGCACGCTGGCGGGTTGCCTGGGTGCCGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCgcagca                                                        <  2:47021/150‑1 (MQ=255)                                                                                             cAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAac                                                                                          >  1:510051/1‑115 (MQ=255)                                                                                             cAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAac                                                                                          <  2:510051/115‑1 (MQ=255)                                                                                                 gcACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTgg                                                   <  1:99048/150‑1 (MQ=255)                                                                                                  tcagCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGa                                                                                                                     >  2:610406/4‑83 (MQ=255)                                                                                                  tcagCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCCGa                                                                                                                     <  1:610406/80‑1 (MQ=255)                                                                                                       tGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAAc                                             <  1:166910/150‑1 (MQ=255)                                                                                                             gTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCgca                                                           <  2:215725/130‑1 (MQ=255)                                                                                                             gTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCgca                                                           >  1:215725/1‑130 (MQ=255)                                                                                                             gTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCa                                       <  2:168089/150‑1 (MQ=255)                                                                                                                  cTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCtt                                  <  2:197450/150‑1 (MQ=255)                                                                                                                       ggCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGct                             <  2:320349/150‑1 (MQ=255)                                                                                                                        gCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGctt                            >  2:336674/1‑150 (MQ=255)                                                                                                                         gggTCACCGTGGATGCC‑GTCGGCACTCACG‑CCTTCGGCCAGCGAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCCAACAGCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTTGCAGAACAGGTCAGCCTTGCGcgtc                           <  2:75904/149‑1 (MQ=255)                                                                                                                         cgctCACCGTCGATGCC‑CACGTCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGACAGCGGCGCGcttc                           <  2:57591/147‑1 (MQ=255)                                                                                                                             cACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGt                       >  1:158401/1‑150 (MQ=255)                                                                                                                              ggcGTGAATGCC‑GACGGCATCCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                      >  1:503391/3‑150 (MQ=255)                                                                                                                              aCCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                      >  2:355208/1‑150 (MQ=255)                                                                                                                               gcGTGAATGCC‑GACGGCATCCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAgc                                                                       <  1:569793/99‑1 (MQ=255)                                                                                                                               gcGTGAATGCC‑GACGGCATCCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAgc                                                                       >  2:569793/2‑100 (MQ=255)                                                                                                                                 gtgaaTGCC‑GACGGCATCCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                      <  2:139074/143‑1 (MQ=255)                                                                                                                                 gTGAATGCC‑GACGGCATCCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                      >  1:139074/1‑147 (MQ=255)                                                                                                                                   ggATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAg                 <  2:293591/150‑1 (MQ=255)                                                                                                                                   ggATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAg                 <  2:150214/150‑1 (MQ=255)                                                                                                                                      tGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCa              <  2:158401/150‑1 (MQ=255)                                                                                                                                                   tCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGCt  >  2:134254/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                       TTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCC‑GTCGGCATTCACG‑CCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGCT  >  NC_002947/5427198‑5427492

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑