Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,980,957 (G)6→7 coding (71/1386 nt) yifK → putative transporter YifK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9511.G93.2% 264.7 / 12.6 73coding (65/1386 nt)yifKputative transporter YifK
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base G (29/39);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

Only 100 of 102 total aligned reads displayed.

CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995
                                                  |                                               
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                  <  2:31492/50‑1 (MQ=255)
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                  <  1:3488811/50‑1 (MQ=255)
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                  <  1:3316704/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                <  1:4796889/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                >  1:5047199/1‑50 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                <  2:2150792/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                <  1:2659910/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                >  1:2469122/1‑50 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgg                                               <  1:4054998/50‑1 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgg                                               <  1:190281/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                            <  1:4319658/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                            <  2:5230450/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                            >  2:3857314/1‑50 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                            <  2:1503010/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                            <  2:3452596/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                            >  1:259286/1‑50 (MQ=255)
     gCTACAGCGTGGGCTGGAAGCTCGGCATATCGAACTCAACGCCCCgggg‑‑g                                         <  2:5843511/50‑1 (MQ=38)
     gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑g                                         >  1:4680032/1‑50 (MQ=255)
     gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑g                                         <  2:2908885/50‑1 (MQ=255)
     gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑g                                         >  1:1665342/1‑50 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑gg                                        >  1:1589877/1‑50 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑gg                                        <  1:3430063/50‑1 (MQ=255)
       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑ggg                                       >  2:5074987/1‑49 (MQ=255)
       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑ggg                                       >  1:5144699/1‑49 (MQ=255)
        aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑ggg                                       >  2:3362877/1‑48 (MQ=255)
         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑gggc                                      <  1:2353164/49‑3 (MQ=255)
         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑ggg                                       <  2:2522617/48‑2 (MQ=255)
         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg‑‑ggg                                       <  1:3205806/48‑2 (MQ=255)
          aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGG‑‑GGCaac                                    <  2:4406052/50‑3 (MQ=255)
            cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCa                                   <  2:3245626/50‑1 (MQ=255)
            cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCa                                   >  1:131049/1‑50 (MQ=255)
            cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCa                                   <  2:2003452/50‑1 (MQ=255)
            cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCa                                   <  1:4542300/50‑1 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCAt                                  >  1:5713703/1‑50 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCAt                                  >  1:4285588/1‑50 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCAt                                  >  1:4275231/1‑50 (MQ=255)
              tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCAtt                                 <  1:2753696/50‑1 (MQ=255)
               gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTg                                >  2:5352750/1‑50 (MQ=255)
               gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTg                                >  2:5237248/1‑50 (MQ=255)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTgg                               <  2:607451/50‑1 (MQ=255)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTg                                <  2:4584017/49‑1 (MQ=39)
                ggCTGGAAGCGCGACATATCGAACTCACCGCCCTGGGGG‑‑GGCACCATTg                                <  2:1734011/49‑1 (MQ=25)
                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGc                              <  2:3951169/50‑1 (MQ=255)
                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGc                              <  2:4937033/50‑1 (MQ=255)
                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGc                              >  2:5101660/1‑50 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCg                             >  1:3526523/1‑50 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCg                             >  2:2879584/1‑50 (MQ=255)
                   tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGt                            >  1:2233269/1‑50 (MQ=255)
                   tGGAAGCTCGACATATCGAACTCATCGCCCTGGGG‑‑GGCACCATTGGCGTc                           >  2:682435/1‑50 (MQ=255)
                    ggAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTc                           >  1:2127362/1‑50 (MQ=255)
                    ggAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTc                           >  1:4864326/1‑50 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTc                           <  2:4193185/49‑1 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCg                          <  2:4182654/50‑1 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCg                          >  1:590377/1‑50 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCg                          <  1:1583496/50‑1 (MQ=255)
                      aaGCTCGACATATCGAACTCATCGCCCTGGTGG‑‑GGCACCATTGGCGTCgg                         >  2:3569922/1‑50 (MQ=255)
                       aGCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGc                        >  2:2439133/1‑50 (MQ=255)
                        gCTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGcc                       <  2:1548391/50‑1 (MQ=255)
                         cTCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCt                      <  1:544731/50‑1 (MQ=255)
                         cTCGACATATCGAACTCATCGCCCTGGGG‑‑GGCACCATTGGCGTCGGCCTg                     <  2:1589877/50‑1 (MQ=255)
                         cTCGACATAACGAACTAATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCt                      <  2:2103636/50‑1 (MQ=255)
                          tCGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTg                     >  2:3060029/1‑50 (MQ=255)
                           cGACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGt                    >  2:1167167/1‑50 (MQ=255)
                           cGACATATCGAACTCATCGCCCTGGGGGTGGCACCATTGGCGTCGGCCTg                     >  2:5228835/1‑50 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGtt                   <  2:639576/50‑1 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGt                    <  1:3208023/49‑1 (MQ=39)
                            gACATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGt                    <  2:4799592/49‑1 (MQ=39)
                             aCATATCGAACTCATCGCCCTGGGGG‑‑GGCGCCATTGGCGTCGGCCTGttt                  <  1:2873105/50‑1 (MQ=255)
                             aCATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGttt                  <  1:800507/50‑1 (MQ=255)
                             aCATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGtt                   >  1:1041462/1‑49 (MQ=39)
                              cATATCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTa                 >  1:3301346/1‑50 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTAt                >  2:3061483/1‑50 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTAt                <  1:1676954/50‑1 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTAt                <  2:5232886/50‑1 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTAt                >  1:4159434/1‑50 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTAt                >  2:1593989/1‑50 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTAt                >  2:3616987/1‑50 (MQ=255)
                                tatCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATg               <  2:817000/50‑1 (MQ=255)
                                tatCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATg               >  1:631916/1‑50 (MQ=255)
                                 atCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATgg              >  1:1627314/1‑50 (MQ=255)
                                  tCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATggg             <  2:3017893/50‑1 (MQ=255)
                                  tCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATggg             <  1:2025956/50‑1 (MQ=255)
                                  tCGAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATggg             <  1:3556418/50‑1 (MQ=255)
                                    gAACTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATggggg           <  1:3492035/50‑1 (MQ=255)
                                     aaCTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGgc          <  2:5165795/50‑1 (MQ=255)
                                     aaCTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGCGGgc          <  2:2283593/50‑1 (MQ=255)
                                      aCTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:3604918/50‑1 (MQ=255)
                                      aCTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:4632910/50‑1 (MQ=255)
                                      aCTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:4703055/50‑1 (MQ=255)
                                      aCTCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:1399039/50‑1 (MQ=255)
                                        tCATCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:5621537/1‑50 (MQ=255)
                                        tCATCGCCCTGGGGG‑‑GGCACCATTGCCGTCGGCCTGCTTATGGGcgccgc       >  2:224844/1‑50 (MQ=255)
                                          aTCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc      <  2:5445070/49‑1 (MQ=39)
                                          aTCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       <  2:3485792/48‑1 (MQ=39)
                                          aTCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:2629266/50‑1 (MQ=255)
                                           tCGCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg    <  2:1986115/50‑1 (MQ=255)
                                             gCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt   <  1:3738736/49‑1 (MQ=39)
                                             gCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  >  2:3266220/1‑50 (MQ=255)
                                             gCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  <  1:2171797/50‑1 (MQ=255)
                                             gCCCTGGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  >  1:593528/1‑50 (MQ=255)
                                                  |                                               
CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGG‑‑GGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995

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: