Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,281,279 A→T A216A (GCA→GCT)  yejH → putative ATP‑dependent DNA or RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,281,2790AT94.4% 243.2 / 5.0 72A216A (GCA→GCT) yejHputative ATP‑dependent DNA or RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (30/38);  total (32/40)
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.01e-01

AGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGATTATGGAGTT  >  NC_000913/2281182‑2281374
                                                                                                 |                                                                                               
aGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGctca                                                                                               <  2:1173425/100‑1 (MQ=255)
 gCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCag                                                                                              <  2:863730/100‑1 (MQ=255)
 gCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCag                                                                                              <  1:1490671/100‑1 (MQ=255)
 gCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACag                                                                                              >  2:1665686/1‑100 (MQ=255)
   tGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCagag                                                                                            <  1:437798/100‑1 (MQ=255)
       gCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCGGAGTAAc                                                                                        >  1:258537/1‑100 (MQ=255)
       gCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCGGAGTAAc                                                                                        >  1:254972/1‑100 (MQ=255)
       gCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAAc                                                                                        >  1:135341/1‑100 (MQ=255)
        cTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACg                                                                                       <  2:1732680/100‑1 (MQ=255)
           cGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGc                                                                                    <  2:1426033/100‑1 (MQ=255)
                tatGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTc                                                                               <  1:847333/100‑1 (MQ=255)
                tatGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTc                                                                               >  1:2124503/1‑100 (MQ=255)
                 atGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGGTTTCAGCCGCTTGCAGGCTCGGAGTAACGGGCTGTTCa                                                                              <  2:254970/100‑1 (MQ=255)
                 atGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCGGAGTAACGGGCTGTTCa                                                                              <  2:258535/100‑1 (MQ=255)
                  tGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAg                                                                             <  2:2212829/100‑1 (MQ=255)
                   gATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGc                                                                            <  2:662503/100‑1 (MQ=255)
                    aTTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCg                                                                           <  2:546655/100‑1 (MQ=255)
                    aTTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCg                                                                           <  2:260030/100‑1 (MQ=255)
                       aaaCACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAGCGGGCTGTTCAGCGAAg                                                                        >  1:196125/1‑100 (MQ=255)
                        aaCACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGc                                                                       >  1:1376176/1‑100 (MQ=255)
                          cacGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCg                                                                     >  1:264071/1‑100 (MQ=255)
                          cacGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCg                                                                     <  2:1582416/100‑1 (MQ=255)
                          cacGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCg                                                                     <  2:1582487/100‑1 (MQ=255)
                            cGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGAt                                                                   >  2:1116632/1‑100 (MQ=255)
                             ggCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGAtc                                                                  <  2:2124496/100‑1 (MQ=255)
                             ggCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGAtc                                                                  <  1:96195/100‑1 (MQ=255)
                               cTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGAtctc                                                                <  2:1468390/100‑1 (MQ=255)
                                tATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCa                                                               >  1:847574/1‑100 (MQ=255)
                                   cTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAAcc                                                            >  1:322310/1‑100 (MQ=255)
                                     gACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGt                                                          >  2:904662/1‑100 (MQ=255)
                                      aCGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTg                                                         <  1:1665691/100‑1 (MQ=255)
                                      aCGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTg                                                         <  1:2106016/100‑1 (MQ=255)
                                        tctgccAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAg                                                       <  1:2199886/97‑1 (MQ=255)
                                        gccgccAGAACGACTGGATATGCCAGTAGTACAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAg                                                       >  2:398550/1‑100 (MQ=255)
                                         ccgccAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGc                                                      >  1:789569/1‑100 (MQ=255)
                                           gccAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAgaga                                                    >  2:2199878/1‑97 (MQ=255)
                                           gccAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTg                                                    >  1:637201/1‑100 (MQ=255)
                                            ccAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGa                                                   <  2:287189/100‑1 (MQ=255)
                                             cAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGaa                                                  >  1:277100/1‑100 (MQ=255)
                                               gAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGaaaa                                                <  1:904664/100‑1 (MQ=255)
                                               gAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGACGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGaaaa                                                <  2:1274814/100‑1 (MQ=255)
                                                aaCGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGaaaaa                                               <  2:2016503/100‑1 (MQ=255)
                                                 aCGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGaaaaaa                                              >  2:1132786/1‑100 (MQ=255)
                                                   gACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaaca                                            <  1:253398/100‑1 (MQ=255)
                                                    aCTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaacaa                                           >  1:1935143/1‑100 (MQ=255)
                                                      tGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaacaaca                                         <  2:1226144/100‑1 (MQ=255)
                                                        gATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaacaacaac                                       <  1:1116634/100‑1 (MQ=255)
                                                         atatGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAATAAACAACAACg                                      <  2:1235035/100‑1 (MQ=255)
                                                         atatGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACg                                      <  1:1072531/100‑1 (MQ=255)
                                                            tGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTAt                                   >  2:1057509/1‑100 (MQ=255)
                                                             gCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTAtt                                  >  2:1918312/1‑100 (MQ=255)
                                                              ccAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTa                                 >  2:1776530/1‑100 (MQ=255)
                                                               cAGTAGCGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTAc                                >  2:1902441/1‑100 (MQ=255)
                                                                 gtagtGCAATACGATTTCAGCCGCTTGCAGGCTCCGAGTAACGGCCTCTTCAGCGAAGCCGATCTCAACCGCGAGCTGAAAAAACAACAACGTCTTAccc                              <  2:1052600/100‑1 (MQ=255)
                                                                  tagtGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTAcccc                             >  1:655475/1‑100 (MQ=255)
                                                                   agtGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCg                            >  1:1794114/1‑100 (MQ=255)
                                                                   agtGCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCg                            >  1:2076870/1‑100 (MQ=255)
                                                                      gCAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGcac                         <  1:1902448/100‑1 (MQ=255)
                                                                       cAATACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGcaca                        <  1:1188948/100‑1 (MQ=255)
                                                                        aaTACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCAcat                       >  1:1008050/1‑100 (MQ=255)
                                                                        aaTACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCAcat                       >  1:1410851/1‑100 (MQ=255)
                                                                         aTACGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCAcatc                      >  2:800114/1‑100 (MQ=255)
                                                                           aCGATTTCAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCAcatcat                    <  1:398552/100‑1 (MQ=255)
                                                                                  cAGCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAg             >  2:877978/1‑100 (MQ=255)
                                                                                    gCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGAt           <  2:1410848/100‑1 (MQ=255)
                                                                                    gCCGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGAt           >  1:57387/1‑100 (MQ=255)
                                                                                     ccGCTTGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGAtt          <  1:1826760/100‑1 (MQ=255)
                                                                                        cTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGATTATg       >  2:625967/1‑100 (MQ=255)
                                                                                         ttGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGATTATgg      <  2:2076863/100‑1 (MQ=255)
                                                                                          tGCAGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGATTATGGa     <  2:1556751/100‑1 (MQ=255)
                                                                                             aGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGATTATGGAGtt  <  2:322308/100‑1 (MQ=255)
                                                                                             aGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGATTATGGAGtt  <  2:2050372/100‑1 (MQ=255)
                                                                                             aGGCTCAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGATTATGGAGtt  <  2:637200/100‑1 (MQ=255)
                                                                                                 |                                                                                               
AGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAACGTATTACCCCGCACATCATCAGCCAGATTATGGAGTT  >  NC_000913/2281182‑2281374

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