Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,427,343 T→C 12.1% 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_0029475,427,3430TC12.1% 104.6 / 5.6 42N270D (AAT→GAT) rhlE‑IATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (17/19);  new base C (3/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

AAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGC  >  NC_002947/5427213‑5427491
                                                                                                                                  |                                                                                                                                                    
aaGTCATCGATATCCACGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATAGCGGCTGCATTCACGCCTTCGGCCAGCAAg                                                                                                                                   <  1:222645/150‑1 (MQ=255)
            tCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGc                                                                                                                       <  1:457891/150‑1 (MQ=255)
             ccAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCt                                                                                                                      >  1:591023/1‑150 (MQ=255)
               aGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCCGCAGGAGGGCCTGTGGGCGGTCACCGTGAATGTGGTCGGCATTCACTCGTTCGGCCAGCAGGCGCTCCACCACCTGa                                                                                                                    >  1:376614/1‑150 (MQ=255)
                    aCGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGAc                                                                                                               >  1:557742/1‑150 (MQ=255)
                         ccgccACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCAt                                                                                                          <  2:557742/150‑1 (MQ=255)
                              acaTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                                     <  1:534582/150‑1 (MQ=255)
                              acaTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGATGGTATGCTCGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                                     <  2:793200/150‑1 (MQ=255)
                                    gTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGc                                                                                                                                          >  2:325765/1‑107 (MQ=255)
                                    gTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGc                                                                                                                                          <  1:325765/107‑1 (MQ=255)
                                                        cTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCCCCTGCtt                                                                           >  2:360843/1‑150 (MQ=255)
                                                             gggCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCa                                                                                                           >  1:327799/1‑113 (MQ=255)
                                                             gggCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCa                                                                                                           <  2:327799/113‑1 (MQ=255)
                                                                ccTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCgct                                                                   <  1:248893/150‑1 (MQ=255)
                                                                  ttGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAgc                                                                      <  2:584509/145‑1 (MQ=255)
                                                                  ttGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAgc                                                                      >  1:584509/1‑145 (MQ=255)
                                                                           tCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCg                                                                                                      >  2:59247/1‑104 (MQ=255)
                                                                           tCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCg                                                                                                      <  1:59247/104‑1 (MQ=255)
                                                                              aGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAgg                                                     >  2:654527/1‑150 (MQ=255)
                                                                                  cACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGc                                                 <  2:591023/150‑1 (MQ=255)
                                                                                   aCGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAgcgc                                                                                                                                <  1:119191/70‑1 (MQ=255)
                                                                                   aCGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAgcgc                                                                                                                                >  2:119191/1‑70 (MQ=255)
                                                                                      cTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGaa                                             <  2:376614/150‑1 (MQ=255)
                                                                                       tGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAAc                                            <  1:436063/150‑1 (MQ=255)
                                                                                           gggTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAggg                                        >  1:131790/1‑149 (MQ=255)
                                                                                                   tgtgGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGc                                                                                                       >  2:400723/1‑79 (MQ=255)
                                                                                                   tgtgGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGc                                                                                                       <  1:400723/79‑1 (MQ=255)
                                                                                                     tgGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTgcg                              >  1:260431/1‑150 (MQ=255)
                                                                                                        gCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCa                                                <  2:508987/129‑1 (MQ=255)
                                                                                                        gCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCa                                                >  1:508987/1‑129 (MQ=255)
                                                                                                           gTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCa                                                                                  >  1:33404/1‑92 (MQ=255)
                                                                                                           gTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCa                                                                                  <  2:33404/92‑1 (MQ=255)
                                                                                                           gTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttctt                        <  2:260431/150‑1 (MQ=255)
                                                                                                              ggcGTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                     >  2:80211/3‑150 (MQ=255)
                                                                                                                cGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCca                   >  1:537559/1‑150 (MQ=255)
                                                                                                                attGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCca                   <  1:273221/148‑1 (MQ=255)
                                                                                                                 gtgaaTGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                     <  1:194323/143‑1 (MQ=255)
                                                                                                                 gTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                     >  2:194323/1‑147 (MQ=255)
                                                                                                                 gTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCcac                  >  1:294738/1‑150 (MQ=255)
                                                                                                                   gaaTGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAg                <  2:294738/148‑1 (MQ=255)
                                                                                                                          gTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttctt                        <  2:249473/135‑1 (MQ=255)
                                                                                                                          gTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttctt                        >  1:249473/1‑135 (MQ=255)
                                                                                                                                 ttCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACATTCACCAGGCACTGc  <  2:131790/150‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                                    
AAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGC  >  NC_002947/5427213‑5427491

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