Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,427,024 A→G 7.3% I376T (ATC→ACC)  rhlE‑I ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,427,0240AG 7.3% 156.9 / 2.8 55I376T (ATC→ACC) rhlE‑IATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (22/29);  new base G (2/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

TCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCGGTACGCCCGATGCGGTGCACGTAGT  >  NC_002947/5426875‑5427171
                                                                                                                                                     |                                                                                                                                                   
tCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGa                                                                                                                                                     >  1:1171364/1‑150 (MQ=255)
 cACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGAt                                                                                                                                                    <  1:1158103/150‑1 (MQ=255)
     gccgGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGaaa                                                                                                                                                <  2:569166/150‑1 (MQ=255)
          ttCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCgg                                                                                                                                           >  2:1095148/1‑150 (MQ=255)
          ttCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCgg                                                                                                                                           <  1:32932/150‑1 (MQ=255)
             cGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGctc                                                                                                                                        <  2:571683/150‑1 (MQ=255)
                 cTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCa                                                                                                                                    <  1:95390/150‑1 (MQ=255)
                   ggCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATg                                                                                                                                  <  2:690330/150‑1 (MQ=255)
                     ccGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGgc                                                                                                                                >  1:109242/1‑150 (MQ=255)
                             gTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGGTGAAATCCGGCTCTTCATGGCGCGGCAgg                                                                                                                        >  2:450421/1‑150 (MQ=255)
                             gTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAgg                                                                                                                        >  2:542899/1‑150 (MQ=255)
                              tGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAggg                                                                                                                       >  2:201198/1‑150 (MQ=255)
                                cTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTc                                                                                                                     >  1:758453/1‑150 (MQ=255)
                                    ccatccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGcc                                                                                                                 <  2:49839/150‑1 (MQ=255)
                                     catccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGGTGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCg                                                                                                                >  2:1186812/1‑150 (MQ=255)
                                     catccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCg                                                                                                                >  2:753745/1‑150 (MQ=255)
                                      atccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCgg                                                                                                               >  1:798488/1‑150 (MQ=255)
                                       tccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCggg                                                                                                              <  1:542899/150‑1 (MQ=255)
                                        ccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGt                                                                                                             <  2:109242/150‑1 (MQ=255)
                                        ccaACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGt                                                                                                             >  1:58315/1‑150 (MQ=255)
                                            cGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGc                                                                                                         <  2:798488/150‑1 (MQ=255)
                                             gCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGGTGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCa                                                                                                        <  1:450421/150‑1 (MQ=255)
                                             gCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCa                                                                                                        <  1:753745/150‑1 (MQ=255)
                                                    cccGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCGCTTCAAt                                                                                                 <  1:943258/150‑1 (MQ=255)
                                                             ggCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGCTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCcag                                                                                        <  1:969183/150‑1 (MQ=255)
                                                               cgcTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGGTGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCcagca                                                                                      <  1:1186812/150‑1 (MQ=255)
                                                                gcTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGa                                                                                                                                                  >  1:785691/1‑89 (MQ=255)
                                                                gcTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGa                                                                                                                                                  <  2:785691/89‑1 (MQ=255)
                                                                   gTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTg                                                                                  >  1:67288/1‑150 (MQ=255)
                                                                     tttctttcGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCa                                                                                <  2:1171364/150‑1 (MQ=255)
                                                                      ttctttcGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTTTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCAc                                                                               <  1:201198/150‑1 (MQ=255)
                                                                         tttcGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCCCTTc                                                                            >  2:856525/1‑150 (MQ=255)
                                                                         tttcGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGAATCAAGCAGGCCCATCATGACTCGGCAGGGTATCCAGGGTAGGAACGGCAGTGTCCGCCAGCAGTTGCACTTc                                                                            >  1:21805/1‑150 (MQ=255)
                                                                                ttcttGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCACTTGCCCTTCAGCCGca                                                                     <  2:900917/150‑1 (MQ=255)
                                                                                       gCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCa                                                              <  2:758453/150‑1 (MQ=255)
                                                                                         ttcttcGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGc                                                            <  1:921102/150‑1 (MQ=255)
                                                                                            ttcGTTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGaa                                                         >  2:1172609/1‑150 (MQ=255)
                                                                                             tcGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGTATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCCCCCACCAGCGcaa                                                        >  2:533539/1‑150 (MQ=255)
                                                                                              cGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAAt                                                       <  1:1095148/150‑1 (MQ=255)
                                                                                                    cTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAAtcgcctc                                                 <  2:58315/150‑1 (MQ=255)
                                                                                                        aTAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGccc                                             <  2:700881/150‑1 (MQ=255)
                                                                                                           aCCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCACACACACCAGCGAAATCGCCTCGCCCttg                                          >  1:737907/1‑150 (MQ=255)
                                                                                                             cTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCttgtt                                        >  2:590794/1‑150 (MQ=255)
                                                                                                              tGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCACACACACCAGCGAAATCGCCTCGCCCttgttg                                       <  2:737907/150‑1 (MQ=255)
                                                                                                               gCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTgc                                      <  1:590794/150‑1 (MQ=255)
                                                                                                                 ccACTGGCATCGGTCATCGGCAGCCGGTGGGCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTgccc                                    <  2:881854/150‑1 (MQ=255)
                                                                                                                   aCTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAgg                                                                                                                        <  2:1107017/64‑1 (MQ=255)
                                                                                                                   aCTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAgg                                                                                                                        >  1:1107017/1‑64 (MQ=255)
                                                                                                                   aCTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTgcccgc                                  >  1:1185403/1‑150 (MQ=255)
                                                                                                                             gTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCGGTa                        >  1:540740/1‑150 (MQ=255)
                                                                                                                             gTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCGGTa                        >  2:851271/1‑150 (MQ=255)
                                                                                                                             gTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCGGTa                        <  2:34078/150‑1 (MQ=255)
                                                                                                                              tCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTgccc                                    <  1:947131/137‑1 (MQ=255)
                                                                                                                              tCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTgccc                                    >  2:947131/1‑137 (MQ=255)
                                                                                                                              tCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCGGTAc                       >  2:963810/1‑150 (MQ=255)
                                                                                                                                                   ggATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCGGTACGCCCGATGCGGTGCACGTAGt  <  1:851271/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTCGGCTTCTTGGGCTTCTTCGGTTTCTTGATAACCTGCCCACTGGCATCGGTCATCGGCACCCGGTGCTCGGGGATGAAATCCGGCTCTTCATGGCGCGGCAGGGTCTGCCGGGTCAGCACTTCAATGGCCGCCAGCAGTTGCACTTCATCCGCACACACCAGCGAAATCGCCTCGCCCTTGTTGCCCGCCCGCCCGGTACGCCCGATGCGGTGCACGTAGT  >  NC_002947/5426875‑5427171

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