Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,979,546 Δ1 bp coding (1047/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,5450G.99.3% 634.4 / ‑3.3 142coding (1046/1386 nt)yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (70/71);  minor base C (1/0);  total (71/71)
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.75e-01

Only 100 of 144 total aligned reads displayed.

AAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGCTCATGCCTGGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCGGGGATGG  >  NC_000913/3979472‑3979616
                                                                          |                                                                       
 aaGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTgaa                                                                       <  2:4271468/76‑3 (MQ=255)
    tttCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTa                                                                   >  2:2763196/1‑76 (MQ=255)
    tttCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTa                                                                   <  2:4437683/76‑1 (MQ=255)
      tCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTAca                                                                 >  1:4510258/1‑76 (MQ=255)
      tCCCGTCACGGCGTACCGGTTGCGGGGGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTAca                                                                 <  2:4496268/76‑1 (MQ=255)
        ccGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGACCTAcatc                                                               >  2:148566/1‑76 (MQ=255)
        ccGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTAcatc                                                               <  1:244422/76‑1 (MQ=255)
         cGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTAcatca                                                              <  1:422728/76‑1 (MQ=255)
          gTCGCGGCGTACCGGATGCGGGTGTGGCAGTATCTATTGCTCTTCTGCTAATTGGCTCATGCCTGAACTAcatcat                                                             <  1:4848711/76‑1 (MQ=255)
            cACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTc                                                           >  1:1736373/1‑76 (MQ=255)
             aCGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTcc                                                          <  1:4831065/76‑1 (MQ=255)
             aCGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTcc                                                          <  1:784568/76‑1 (MQ=255)
              cGGCGTGCCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTccc                                                         <  2:203952/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: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                                                     <  2:469123/76‑1 (MQ=255)
                  gTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATc                                                     >  1:1179430/1‑76 (MQ=255)
                   tgccGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATcc                                                    <  2:1273983/74‑1 (MQ=255)
                   tACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATcc                                                    >  1:4860689/1‑76 (MQ=255)
                   tACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATcc                                                    <  1:4337559/76‑1 (MQ=255)
                   tACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATcc                                                    <  2:809298/76‑1 (MQ=255)
                    aCCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCCTTCCCAAcccc                                                   >  2:1747612/1‑75 (MQ=255)
                     ccGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGc                                                  >  1:3144895/1‑76 (MQ=255)
                     ccGGTTGCGGGTGTGGCAGTATATATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCAGTCCCCATcctc                                                  >  1:3436547/1‑74 (MQ=255)
                      cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCCTGCCTGAACTACATCATTCCCACTCCGcc                                                 >  2:1530662/1‑75 (MQ=255)
                      cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGACATACATCATTCCCCATCTGCa                                                 >  1:2609499/1‑76 (MQ=255)
                      cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATccccc                                                 >  2:3214992/1‑73 (MQ=255)
                      cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCa                                                 <  1:4318084/76‑1 (MQ=255)
                      cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCa                                                 >  1:2601785/1‑76 (MQ=255)
                      cGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCa                                                 <  1:3328470/76‑1 (MQ=255)
                      cGGTTGCGCGTGTGCCAGTATCTATTGCTATTCTTCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCa                                                 <  1:3343407/76‑1 (MQ=255)
                         ttGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCcgcg                                              >  2:4329093/1‑76 (MQ=255)
                         ttGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCg                                              <  1:870035/76‑1 (MQ=255)
                          tGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgt                                             >  2:2825553/1‑76 (MQ=255)
                           gCGGGTGTGGCAGTCTCTATTGCTATTCTGCTAATTGGCTCATGGCTGAACTCCATCATTCCCAATCCGCAGCgtg                                            <  2:4487685/76‑1 (MQ=255)
                            cGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgtgt                                           <  2:4371251/76‑1 (MQ=255)
                            cGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCACCgtgt                                           >  2:4391472/1‑76 (MQ=255)
                            cGGGTGTCGCAGTATCTATTGCTATTCTGCTAATTGGCTCATTCCTGAACTACATCATGCCCAATCCGCAGCgtgt                                           <  1:3468633/76‑1 (MQ=255)
                             gggTGTGGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgtgtg                                          >  1:1937187/1‑76 (MQ=255)
                                tgtgGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAg                                                <  2:3129668/67‑1 (MQ=255)
                                tgtgGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAg                                                >  1:3129668/1‑67 (MQ=255)
                                tgtgGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGttt                                       >  2:1080039/1‑76 (MQ=255)
                                tgtgGCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGttt                                       >  2:1901361/1‑76 (MQ=255)
                                ggtgGCAGGATCTATTGCTATTCTGCGAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGttt                                       <  2:699248/75‑1 (MQ=255)
                                   ggCAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtc                                    >  1:2048408/1‑76 (MQ=255)
                                     cAGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtcta                                  <  1:2590892/76‑1 (MQ=255)
                                     cAGTAACTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtcta                                  >  1:3586323/1‑76 (MQ=255)
                                      aGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctac                                 <  1:1561019/76‑1 (MQ=255)
                                      aGTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctac                                 <  1:1942693/76‑1 (MQ=255)
                                       gTATCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacg                                <  1:2793398/76‑1 (MQ=255)
                                         aTCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtc                              >  1:4243793/1‑76 (MQ=255)
                                          tCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtct                             <  2:2242837/76‑1 (MQ=255)
                                          tCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtct                             >  2:4482314/1‑76 (MQ=255)
                                          tCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtct                             >  1:3031488/1‑76 (MQ=255)
                                          tCTATTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTgtctacgtct                             >  2:413217/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)
                                             aTTGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACa                          <  2:4064196/76‑1 (MQ=255)
                                              ttGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAg                         >  2:2837102/1‑76 (MQ=255)
                                               tGCTATTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGt                        >  1:2343321/1‑76 (MQ=255)
                                                gCTATTCTGCTACTTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTg                       >  2:1333882/1‑76 (MQ=255)
                                                  tATTCTGCTAATTGGCGCATGCCTGAACTACATCACTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGcc                     <  1:1230088/76‑1 (MQ=255)
                                                   aTTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCa                    <  2:3764672/76‑1 (MQ=255)
                                                   aTTCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGCCTCCCGTGCCa                    >  2:111943/1‑76 (MQ=255)
                                                   aTTCTGATAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCa                    <  2:949843/76‑1 (MQ=255)
                                                    ttCTGCTAATTGGCTCATGCCTGACCTACATCATCCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAg                   <  1:4568945/76‑1 (MQ=255)
                                                     tCTGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGc                  <  1:639053/76‑1 (MQ=255)
                                                       tGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt                >  2:2137396/1‑76 (MQ=255)
                                                       tGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt                <  1:69807/76‑1 (MQ=255)
                                                       tGCTAATTGGCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt                <  1:893355/76‑1 (MQ=255)
                                                       tGCTAATTGGCTCAGGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGt                <  2:3698322/76‑1 (MQ=255)
                                                             ttggCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTcc          <  1:2918019/73‑1 (MQ=255)
                                                             ttggCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTcc          >  2:3879918/4‑76 (MQ=255)
                                                              tggCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCg         <  1:1158083/74‑1 (MQ=255)
                                                              tggCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCg         >  1:1576550/3‑76 (MQ=255)
                                                               ggCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgt                                             <  1:4440296/38‑1 (MQ=38)
                                                               ggCTCATGCCTGAACTACATCATTCCCAATCCGCAGCgt                                             >  2:4440296/2‑39 (MQ=38)
                                                               ggCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCgg        <  1:3932975/75‑1 (MQ=255)
                                                               ggCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCgg        >  1:3549840/2‑76 (MQ=255)
                                                                gCTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCggg       <  2:1003383/76‑1 (MQ=255)
                                                                 cTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCgggg      <  2:2046254/76‑1 (MQ=255)
                                                                 cTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCgggg      <  1:4326136/76‑1 (MQ=255)
                                                                 cTCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCgggg      <  2:2404995/76‑1 (MQ=255)
                                                                  tCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCGGGGa     <  1:3568563/76‑1 (MQ=255)
                                                                  tCATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCATGCTTCCGGGGa     >  2:2144377/1‑76 (MQ=255)
                                                                   cATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCGGGGAt    >  1:4277918/1‑76 (MQ=255)
                                                                   cATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCGGGGAt    >  1:2755923/1‑76 (MQ=255)
                                                                   cATGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCGGGGAt    >  1:2754557/1‑76 (MQ=255)
                                                                    aTGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCGGGGATg   >  2:3512140/1‑76 (MQ=255)
                                                                     tGCCTGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCGGGGATgg  <  2:2732326/76‑1 (MQ=255)
                                                                          |                                                                       
AAAGTTTCCCGTCACGGCGTACCGGTTGCGGGTGTGGCAGTATCTATTGCTATTCTGCTAATTGCTCATGCCTGGAACTACATCATTCCCAATCCGCAGCGTGTGTTTGTCTACGTCTACAGTGCCAGCGTGCTTCCGGGGATGG  >  NC_000913/3979472‑3979616

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