Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,715,660 (G)5→6 coding (238/588 nt) eamB → cysteine/O‑acetylserine exporter EamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,715,6551.G75.2% 314.4 / 111.5 109coding (233/588 nt)eamBcysteine/O‑acetylserine exporter EamB
Reads supporting (aligned to +/- strand):  ref base . (12/15);  new base G (28/54);  total (40/69)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 141 total aligned reads displayed.

TTTTCACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCGCCACCAGCCC  >  NC_000913/2715606‑2715703
                                                   |                                                 
ttttCACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggc                                                    <  1:2906844/50‑1 (MQ=255)
ttttCACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggc                                                    <  2:182307/50‑1 (MQ=255)
ttttCACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggc                                                    >  2:1959748/1‑50 (MQ=255)
 tttCACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcg                                                  >  1:71773/1‑50 (MQ=255)
  ttCACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcgg                                                 >  2:2839663/1‑50 (MQ=255)
  ttCACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcg                                                  <  2:4173955/49‑1 (MQ=255)
   tCACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcggg                                                >  1:4252413/1‑50 (MQ=255)
    cACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcgggg                                              >  2:110545/1‑50 (MQ=255)
    cACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcgggg                                              >  2:4478087/1‑50 (MQ=255)
    cACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcgggg                                              >  1:50873/1‑50 (MQ=255)
    cACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcgggg                                              >  1:2421950/1‑50 (MQ=255)
    cACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcggg                                                <  2:181062/49‑1 (MQ=255)
     aCTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcggggg                                             >  1:2898846/1‑50 (MQ=255)
     aCTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcgggg                                              <  1:1716650/49‑1 (MQ=255)
      cTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgtgcgggggg                                            <  2:1352734/50‑2 (MQ=255)
      cTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcgggggg                                            >  1:1871516/1‑49 (MQ=255)
       tGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcggggggc                                           >  1:138329/1‑48 (MQ=255)
       tGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcgggggcg                                           >  1:44107/1‑50 (MQ=255)
       tGGCAGTGATTGACCCGGCAGCGGTACACCTTTTCAGTTgggcggggggc                                           >  1:4596542/1‑48 (MQ=39)
         gCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcggggggc                                           <  2:2619754/48‑3 (MQ=255)
         gCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcggggggc                                           <  1:3209113/48‑3 (MQ=255)
         gCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTgggcggggggc                                           <  1:349396/48‑3 (MQ=255)
          cAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCa                                        >  1:213087/1‑50 (MQ=255)
           aGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCa                                        >  1:2030716/1‑50 (MQ=255)
           aGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCa                                        <  1:1206986/50‑1 (MQ=255)
           aGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCa                                        >  2:3998089/1‑50 (MQ=255)
            gTGATTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCat                                       <  2:5553801/50‑1 (MQ=255)
             tGATTGACCCGGCGGCGGTACACCTTTGGAGTTGGGCGGGGGGCGGCata                                      <  2:5432216/50‑1 (MQ=255)
               aTTGACCCGGCAGCGGTACGCCTTTTGAGTTGGGCGGGGGGCGGCatata                                    <  2:4852694/50‑1 (MQ=255)
               aTTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCatata                                    >  1:1290243/1‑50 (MQ=255)
               aTTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCatata                                    <  2:5303555/50‑1 (MQ=255)
               aTTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCatat                                     <  1:2476499/49‑1 (MQ=255)
               aTTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCatat                                     <  2:1567928/49‑1 (MQ=255)
               aTTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCatatat                                   <  1:942455/50‑1 (MQ=255)
               aTTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCatata                                    <  2:4178000/49‑1 (MQ=255)
               aTTGACCCGGCAGCGGTACACCTTTTGAGGTGGGCGGGGGGCGGCatata                                    >  2:1624247/1‑50 (MQ=255)
                ttGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCatatat                                   <  1:5628771/50‑1 (MQ=255)
                ttGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATAtt                                  <  1:1567228/50‑1 (MQ=255)
                 tGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATAtt                                  >  1:771277/1‑50 (MQ=255)
                  gACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTg                                 >  2:4366303/1‑50 (MQ=255)
                  gACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTg                                 <  2:4458215/49‑1 (MQ=255)
                   tcccGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGt                                <  1:3170415/49‑1 (MQ=255)
                   aCCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGt                                <  1:3047687/50‑1 (MQ=255)
                   aCCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGt                                <  1:4485274/50‑1 (MQ=255)
                   aCCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTGt                                <  2:355064/49‑1 (MQ=255)
                   aCCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTGTc                               >  2:779177/1‑50 (MQ=255)
                    cccGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTc                               <  2:656427/50‑1 (MQ=255)
                     ccGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTct                              >  2:476147/1‑50 (MQ=255)
                     ccGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTct                              <  2:2601210/50‑1 (MQ=255)
                     ccGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTGTctg                             <  2:4809160/50‑1 (MQ=255)
                      cGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctg                             <  2:1086228/50‑1 (MQ=255)
                      cGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctg                             >  1:3571613/1‑50 (MQ=255)
                      cGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctg                             <  1:2328970/50‑1 (MQ=255)
                      cGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctg                             <  2:3148084/50‑1 (MQ=255)
                      cGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctg                             <  1:3667817/50‑1 (MQ=255)
                      cGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTGTctgg                            <  2:3117604/50‑1 (MQ=255)
                       ggCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctgg                            <  1:5132429/50‑1 (MQ=255)
                       ggCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctgg                            >  2:3115839/1‑50 (MQ=255)
                        gCAGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctggc                           >  1:2403103/1‑50 (MQ=255)
                         cAGCGGTACGCCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctggct                          <  2:933030/50‑1 (MQ=255)
                         cAGCGGTACACCTGTTGAGTTGGGCGGGGGGCGGCATGTATTGTgtggct                          <  2:900558/50‑1 (MQ=255)
                          aGCGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTctggctg                         <  1:5305718/50‑1 (MQ=255)
                          aGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTGTctggctgg                        >  1:5485427/1‑50 (MQ=255)
                          aGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTGTctggctgg                        >  2:421344/1‑50 (MQ=255)
                            cGGTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCtggctggc                       >  2:5069667/1‑50 (MQ=255)
                             ggTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCg                      >  2:3692365/1‑50 (MQ=255)
                             ggTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTGTCTGGCTGGCGt                     <  1:2596115/50‑1 (MQ=255)
                              gTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGt                     <  1:1749685/50‑1 (MQ=255)
                              gTACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGt                     >  1:207499/1‑50 (MQ=255)
                              gTACACCTTTTGAGTTGGGCGGGTGGCGGCATATATTGTCTGGCTGGCGt                     <  1:3267778/50‑1 (MQ=255)
                              gTACACCTTTTGAGTGGGCGGGGGGGCGGCGGATATTGCCTGGCTGCCGt                     >  2:377786/1‑50 (MQ=255)
                               tACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTg                    >  1:3639807/1‑50 (MQ=255)
                               tACACCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTg                    <  1:5113190/50‑1 (MQ=255)
                                acacCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTg                    >  2:317731/1‑49 (MQ=39)
                                 cacCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTgg                   >  1:4642126/1‑49 (MQ=39)
                                 cacCTTTTGAGTTGGGCGGGGGCGGCATATATTGTCTGGCTGGCGTGGa                  <  1:3431097/49‑1 (MQ=255)
                                  acCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGaa                 <  2:3535441/50‑1 (MQ=255)
                                  acCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGaa                 <  2:5425517/50‑1 (MQ=255)
                                  acCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGaa                 <  2:4873331/50‑1 (MQ=255)
                                  acCTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGaa                 <  1:180358/50‑1 (MQ=255)
                                   ccTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGaaa                <  1:3921130/50‑1 (MQ=255)
                                    cTTTTGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGaaaa               <  1:2263556/50‑1 (MQ=255)
                                    cTTTTGAGTTGGGCGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAAt              <  1:5598826/50‑1 (MQ=255)
                                    cTTTTGAGTTGGGCGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAAt              >  1:5345236/1‑50 (MQ=255)
                                      tttGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATc             >  1:3420566/1‑50 (MQ=255)
                                      tttGAGTTGGGCGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCg            >  1:4590590/1‑50 (MQ=255)
                                       ttGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAAt              <  1:439950/48‑1 (MQ=39)
                                       ttGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAAt              <  1:4860594/48‑1 (MQ=39)
                                       ttGAGTTGGGCGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATc             <  1:2065012/48‑1 (MQ=255)
                                        tGAGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCg            <  1:2733227/49‑1 (MQ=39)
                                          aGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCGcca         >  1:2743695/1‑50 (MQ=255)
                                          aGTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCGcc          <  1:3740100/49‑1 (MQ=39)
                                           gTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCGccac        >  2:4239116/1‑50 (MQ=255)
                                           gTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCGccac        <  2:1638258/50‑1 (MQ=255)
                                           gTTGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCGccac        <  1:1756317/50‑1 (MQ=255)
                                           gTTGGGCGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCGccacc       <  1:1465262/50‑1 (MQ=255)
                                             tGGGCGGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCGccac        >  2:3719844/1‑48 (MQ=39)
                                              gggcgggggcggCATATATTGTCTGGCTAGCGTGGAAAATCGCCACCAg     <  2:3458425/49‑1 (MQ=255)
                                               ggcgggggcggCATATATTGTCTGGCTGGCGTGGAAAATCGCCACCAGcc   <  2:3776551/50‑1 (MQ=255)
                                                 gcgggggcggCATATATTGTCTGGCTGGCGTGGAAAATCGCCACCAGccc  <  2:4938063/50‑1 (MQ=255)
                                                   |                                                 
TTTTCACTGGCAGTGATTGACCCGGCAGCGGTACACCTTTTGAGTTGGGCGGGGGCGGCATATATTGTCTGGCTGGCGTGGAAAATCGCCACCAGCCC  >  NC_000913/2715606‑2715703

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