Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,979,535 +G coding (1036/1386 nt) yifK → putative amino acid/amine transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,979,5341.G98.4% 339.2 / ‑1.0 128coding (1035/1386 nt)yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (62/64);  total (64/64)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 137 total aligned reads displayed.

GCGGCAATGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGCTCATGCCTGGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGT  >  NC_000913/3979460‑3979602
                                                                            |                                                                    
gcggcAATGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTg                                                                     <  2:1856994/76‑1 (MQ=255)
gaggcAATCGAGAAAGGGTCCAGTCACGGCGTGACGGTGGAGGGTGTCGCAGTAGAGATTGCTATTCTGCTGAtta                                                                      <  2:1886808/74‑2 (MQ=255)
 cggcAATGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTgg                                                                    >  1:1928039/1‑75 (MQ=255)
 cggcAATGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTgg                                                                    <  2:3117500/76‑2 (MQ=255)
  ggcAATGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTggc                                                                   <  2:2340059/76‑3 (MQ=255)
   gcAATGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTggct                                                                  <  2:56389/76‑4 (MQ=255)
      aTGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCAt                                                                <  1:4748352/76‑1 (MQ=255)
       tGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATg                                                               <  2:3812561/76‑1 (MQ=255)
       tGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATg                                                               >  1:4859675/1‑76 (MQ=255)
        ggCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGCTGTGGCAGTATCTATTGCTATTCTGCTAATTTGCTCACGc                                                              >  2:1497510/1‑76 (MQ=255)
             aaGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTgaa                                                         <  2:4271468/76‑3 (MQ=255)
                tttCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTa                                                     <  1:3664168/76‑1 (MQ=255)
                  tCCCGTCACGGCGTACCGGTTGCGGGGGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTAca                                                   <  2:4496268/76‑1 (MQ=255)
                    ccGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGACCTAcatc                                                 >  2:148566/1‑76 (MQ=255)
                    ccGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTAcatc                                                 <  1:244422/76‑1 (MQ=255)
                    ccGTCACGGCGTACCGGTTGCGGGCGTGGCAGCATCTATGGCTATTCTTCTCATTGGCTCATGCCTGAACTTcatc                                                 >  2:2165470/1‑76 (MQ=255)
                     cGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTTCTAATTGGCTCATGCCTGAACCAcatca                                                >  1:2605568/1‑76 (MQ=255)
                     cGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTAcatca                                                <  1:422728/76‑1 (MQ=255)
                      gTCGCGGCGTACCGGATGCGGGTGTGGCAGTATCTATTGCTCTTCTGCTAATTGGCTCATGCCTGAACTAcatcat                                               <  1:4848711/76‑1 (MQ=255)
                        cACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTc                                             >  1:1736373/1‑76 (MQ=255)
                         aCGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCAtttc                                            >  1:847008/1‑74 (MQ=255)
                         aCGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTcc                                            <  1:784568/76‑1 (MQ=255)
                          cGGCGTGCCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTccc                                           <  2:203952/76‑1 (MQ=255)
                          cGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTccc                                           <  1:1672321/76‑1 (MQ=255)
                          cGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTccc                                           >  2:3464608/1‑76 (MQ=255)
                           ggggTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCa                                          <  2:4570222/73‑1 (MQ=255)
                           ggCGTGCCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCa                                          <  2:1531361/76‑1 (MQ=255)
                           ggCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCa                                          >  1:3109249/1‑76 (MQ=255)
                           ggCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCa                                          >  1:3549889/1‑76 (MQ=255)
                           ggCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCa                                          >  1:2672523/1‑76 (MQ=255)
                           ggCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCa                                          >  2:637114/1‑76 (MQ=255)
                            gCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCaa                                         <  1:565691/76‑1 (MQ=255)
                             tgTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCaac                                        >  2:374388/2‑75 (MQ=255)
                             cGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAAt                                        >  1:3723829/1‑76 (MQ=255)
                              gTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCCATc                                       >  2:2367813/1‑76 (MQ=255)
                              gTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATc                                       >  1:1179430/1‑76 (MQ=255)
                              gTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATc                                       <  2:1761527/76‑1 (MQ=255)
                              gTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCCGAACTACATCATTCCCAATc                                       >  2:2614326/1‑76 (MQ=255)
                               tgccGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATcc                                      <  2:1273983/74‑1 (MQ=255)
                               tACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATcc                                      <  2:809298/76‑1 (MQ=255)
                               tACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATcc                                      >  1:4860689/1‑76 (MQ=255)
                               tACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATcc                                      <  1:4337559/76‑1 (MQ=255)
                               tACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATcc                                      <  1:4283187/76‑1 (MQ=255)
                                 ccGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGc                                    >  1:3144895/1‑76 (MQ=255)
                                  cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCCTGCCTGAACTACATCATTCCCACTCCGcc                                   >  2:1530662/1‑75 (MQ=255)
                                  cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGACATACATCATTCCCCATCTGCa                                   >  1:2609499/1‑76 (MQ=255)
                                  cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATccccc                                   >  2:3214992/1‑73 (MQ=255)
                                  cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCa                                   <  2:4470675/76‑1 (MQ=255)
                                  cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCa                                   <  1:3328470/76‑1 (MQ=255)
                                  cGGTTGCGCGTGTGCCAGTATCTATTGCTATTCTTCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCa                                   <  1:3343407/76‑1 (MQ=255)
                                   ggTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCCTTCCCCATCCGCAg                                  >  2:3115430/1‑76 (MQ=255)
                                     ttGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCcgcg                                >  2:4329093/1‑76 (MQ=255)
                                      tGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgt                               >  2:2825553/1‑76 (MQ=255)
                                        cGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgtgt                             <  2:4371251/76‑1 (MQ=255)
                                        cGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCACCgtgt                             >  2:4391472/1‑76 (MQ=255)
                                        cGGGTGTCGCAGTATCTATTGCTATTCTGCTAATTGGCTCATTCCTGAACTACATCATGCCCAATCCGCAGCgtgt                             <  1:3468633/76‑1 (MQ=255)
                                         gggTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgtgtg                            >  2:2774526/1‑76 (MQ=255)
                                         gggTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgtgtg                            >  1:1937187/1‑76 (MQ=255)
                                         gggTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgtgtg                            <  2:4522105/76‑1 (MQ=255)
                                           gtgtGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGtt                          >  1:3353027/1‑76 (MQ=255)
                                           gtgtGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGtt                          <  2:3729110/76‑1 (MQ=255)
                                            tgtgGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAg                                  >  1:3129668/1‑67 (MQ=255)
                                            tgtgGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGttt                         >  2:1080039/1‑76 (MQ=255)
                                             gtgGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTg                        >  1:2274171/1‑76 (MQ=255)
                                               ggCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtc                      >  1:2048408/1‑76 (MQ=255)
                                               ggCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtc                      >  1:2810191/1‑76 (MQ=255)
                                                 cAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtcta                    <  1:2590892/76‑1 (MQ=255)
                                                 cAGTAACTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtcta                    >  1:3586323/1‑76 (MQ=255)
                                                  aGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctac                   <  2:3977932/76‑1 (MQ=255)
                                                   gTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacg                  >  2:2024954/1‑76 (MQ=255)
                                                   gTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacg                  <  1:2793398/76‑1 (MQ=255)
                                                     aTCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtc                >  1:4243793/1‑76 (MQ=255)
                                                     aTCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtc                >  1:4497154/1‑76 (MQ=255)
                                                      tCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtct               >  2:4400499/1‑76 (MQ=255)
                                                      tCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtct               >  2:4482314/1‑76 (MQ=255)
                                                      tCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtct               >  1:1765714/1‑76 (MQ=255)
                                                      tCTATTGCTATTCTGCTAATTGGCGCATGCCTGACCTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtct               <  2:3784235/76‑1 (MQ=255)
                                                       cTATTGCTCTTGTGCTAATTGGCTCATGCCTGAACTACATCATGCCCAATCCGCAGCGTGTGTTTgtctacgtcta              <  2:4045369/76‑1 (MQ=255)
                                                       cTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtcta              <  2:1992241/76‑1 (MQ=255)
                                                       cTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCACCGTGTGTTTctctacctcta              >  2:1442204/1‑76 (MQ=255)
                                                        tATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtctac             <  1:4332655/76‑1 (MQ=255)
                                                         aTTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACa            <  2:4064196/76‑1 (MQ=255)
                                                          ttGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAg           >  2:2837102/1‑76 (MQ=255)
                                                            gCTATTCTGCTACTTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTg         >  2:1333882/1‑76 (MQ=255)
                                                            gCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTg         >  1:2573154/1‑76 (MQ=255)
                                                            gCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTg         >  1:1375436/1‑76 (MQ=255)
                                                             cTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGc        >  2:605212/1‑76 (MQ=255)
                                                              tATTCTGGTAATTGGCGCATGCCTGAATTACAGCGTTCCCAATGCGCAGCGTGTGTTTGTCTACGTCTCCAGTGcc       <  2:715068/76‑1 (MQ=255)
                                                              tATTCTGCTAATTGGCGCATGCCTGAACTACATCACTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGcc       <  1:1230088/76‑1 (MQ=255)
                                                               aTTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCa      <  2:3764672/76‑1 (MQ=255)
                                                                ttCTGCTAAGTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAg     >  2:1343286/1‑76 (MQ=255)
                                                                ttCTGCTAATTGGCTCATGCCTGACCTACATCATCCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAg     <  1:4568945/76‑1 (MQ=255)
                                                                 tCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGc    <  2:4229225/76‑1 (MQ=255)
                                                                 tCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGc    <  1:639053/76‑1 (MQ=255)
                                                                   tGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt  >  2:3457360/1‑76 (MQ=255)
                                                                   tGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt  >  2:2137396/1‑76 (MQ=255)
                                                                   tGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt  <  1:930914/76‑1 (MQ=255)
                                                                   tGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt  <  1:893355/76‑1 (MQ=255)
                                                                   tGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt  <  1:69807/76‑1 (MQ=255)
                                                                   tGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt  <  1:1900712/76‑1 (MQ=255)
                                                                            |                                                                    
GCGGCAATGGCGAAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGCTCATGCCTGGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGT  >  NC_000913/3979460‑3979602

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: