Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,426,955 A→G 8.1% M383T (ATG→ACG)  rhlE‑I ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,426,9550AG 8.1% 134.6 / 2.4 49M383T (ATG→ACG) rhlE‑IATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base G (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCGG  >  NC_002947/5426809‑5427099
                                                                                                                                                  |                                                                                                                                                
gCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGAATCGGTCATc                                                                                                                                                <  4:311349/149‑1 (MQ=255)
                 ttCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCGTCGGCACCCgg                                                                                                                                       <  8:270889/141‑1 (MQ=255)
                 ttCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCGTCGGCACCCgg                                                                                                                                       >  7:270889/1‑141 (MQ=255)
                 ttCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCggg                                                                                                                               >  7:148763/1‑149 (MQ=255)
                   cACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCggg                                                                                                                               <  5:29437/147‑1 (MQ=255)
                    aCCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGAt                                                                                                                            >  3:94833/1‑149 (MQ=255)
                      cgccgGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATg                                                                                                                           >  5:301301/1‑148 (MQ=255)
                              ccGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCgg                                                                                                                   <  6:301301/148‑1 (MQ=255)
                              ccGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCgg                                                                                                                   <  6:175421/148‑1 (MQ=255)
                                gCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGctc                                                                                                                >  8:208400/1‑149 (MQ=255)
                                      cccGCTTCAGTGCTGTCCATCAAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATg                                                                                                          >  7:148375/2‑149 (MQ=255)
                                        cGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGgc                                                                                                        >  3:326007/1‑149 (MQ=255)
                                                  cTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATc                                                                                                                                                <  6:226175/99‑1 (MQ=255)
                                                  cTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATc                                                                                                                                                >  5:226175/1‑99 (MQ=255)
                                                       catccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCa                                                                                                  >  3:253850/1‑140 (MQ=255)
                                                       catccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCa                                                                                                  <  4:253850/140‑1 (MQ=255)
                                                       catccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGcc                                                                                         >  1:254641/1‑149 (MQ=255)
                                                       catccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGcc                                                                                         >  4:280117/1‑149 (MQ=255)
                                                        atccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCg                                                                                        >  1:186591/1‑149 (MQ=255)
                                                        atccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCg                                                                                        >  5:238603/1‑149 (MQ=255)
                                                        atccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCg                                                                                        <  3:280117/149‑1 (MQ=255)
                                                          ccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCGTCGGCACCCGGTGCTCGGGGATGAAATCCgg                                                                                                                   <  4:183961/120‑1 (MQ=255)
                                                          ccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCGTCGGCACCCGGTGCTCGGGGATGAAATCCgg                                                                                                                   >  3:183961/1‑120 (MQ=255)
                                                           caACGCCGCAGCCGGCGCTTGGCGCTGGTTTCTTTCGGCTTCTTGGGCTTCTTCTGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGt                                                                                     <  8:148763/149‑1 (MQ=255)
                                                           caACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGt                                                                                     <  7:208400/149‑1 (MQ=255)
                                                              cGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAg                                                                                  <  5:118120/149‑1 (MQ=255)
                                                               gCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGc                                                                                 <  2:317496/149‑1 (MQ=255)
                                                                    aGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCgg                                                                                       >  4:267649/1‑138 (MQ=255)
                                                                    aGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCgg                                                                                       <  3:267649/138‑1 (MQ=255)
                                                                       ccGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAAt                                                                         <  4:94833/149‑1 (MQ=255)
                                                                        cgcgCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATg                                                                        >  1:223190/1‑149 (MQ=255)
                                                                               ggCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGgccgc                                                                   <  4:326007/147‑1 (MQ=255)
                                                                                gcgcTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCca                                                                 <  2:186591/148‑1 (MQ=255)
                                                                                                gCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATcc                                                <  2:254641/149‑1 (MQ=255)
                                                                                                        ggCTTCTTCGGTTTCTTGATAACCTGCCCACTGTCATGGGTCATCGGCACCCGGTGCTCGGGGATGAATTCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATg                                                                        <  5:88051/117‑1 (MQ=255)
                                                                                                             cttcGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCa                                      >  4:313603/1‑146 (MQ=255)
                                                                                                             cttcGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCa                                      <  3:313603/146‑1 (MQ=255)
                                                                                                             cttcGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGc                                    <  2:223190/148‑1 (MQ=255)
                                                                                                             cttcGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGc                                    <  8:148375/148‑1 (MQ=255)
                                                                                                                   tttCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAAtcg                             >  4:251052/1‑149 (MQ=255)
                                                                                                                              ccTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTg                                                          >  1:138297/1‑109 (MQ=255)
                                                                                                                              ccTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTg                                                          <  2:138297/109‑1 (MQ=255)
                                                                                                                               cTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCttgt                 >  4:240703/1‑149 (MQ=255)
                                                                                                                                         gCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTgcccgcccg       >  1:247600/1‑149 (MQ=255)
                                                                                                                                           aTCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAAtcgcct                          >  2:322063/1‑128 (MQ=255)
                                                                                                                                            tCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTgcccgcccgccc    >  2:203552/1‑149 (MQ=255)
                                                                                                                                              ggTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGcccgcccgcccg   <  6:238603/148‑1 (MQ=255)
                                                                                                                                              ggTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCgg  >  4:205883/1‑149 (MQ=255)
                                                                                                                                               gTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCgg  >  6:249454/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCGG  >  NC_002947/5426809‑5427099

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