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 APC family amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9511.G98.3% 217.3 / ‑3.8 60coding (65/1386 nt)yifKputative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (32/27);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAA  >  NC_000913/3980827‑3981082
                                                                                                                             |                                                                                                                                    
gATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCAtt                                                                                                                        >  2:312036/1‑139 (MQ=255)
      cccATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTc                                                                                                                  >  1:114851/1‑139 (MQ=255)
       ccATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCg                                                                                                                 <  1:426844/139‑1 (MQ=255)
        cATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCgg                                                                                                                >  1:681411/1‑139 (MQ=255)
         aTAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGc                                                                                                               <  2:205926/139‑1 (MQ=255)
          tAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGcc                                                                                                              >  1:86286/1‑139 (MQ=255)
          tAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGcc                                                                                                              >  1:437730/1‑139 (MQ=255)
              gATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGtt                                                                                                          >  2:280932/1‑139 (MQ=255)
                tAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTa                                                                                                        >  1:738717/1‑139 (MQ=255)
                 aaCCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTAt                                                                                                       >  1:337260/1‑139 (MQ=255)
                 aaCCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTAt                                                                                                       >  2:685991/1‑139 (MQ=255)
                     ggTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                   <  2:681411/139‑1 (MQ=255)
                       tAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                                                                 >  1:260273/1‑139 (MQ=255)
                        aaCACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc                                                                                                >  2:680047/1‑139 (MQ=255)
                          cacCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc                                                                                              >  1:34997/1‑139 (MQ=255)
                          cacCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc                                                                                              <  1:337101/139‑1 (MQ=255)
                          cacCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc                                                                                              >  1:420669/1‑139 (MQ=255)
                            ccGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                                                            <  2:767412/139‑1 (MQ=255)
                             cGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                                                                           >  2:655836/1‑139 (MQ=255)
                                      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCTTGGG‑GGGGCACCATTGGCGTCGGc                                                                                                               >  2:754990/1‑110 (MQ=255)
                                      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGc                                                                                                               <  1:754990/110‑1 (MQ=255)
                                      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaa                                                                                  <  1:586484/139‑1 (MQ=255)
                                                 cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       <  1:238548/139‑1 (MQ=255)
                                                 cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       >  2:231304/1‑139 (MQ=255)
                                                 cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       >  2:579545/1‑139 (MQ=255)
                                                  aggAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGGACCCTGAAATgggccgggcc                                                                      <  1:470223/137‑1 (MQ=255)
                                                       gATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCg                                                                 <  2:260273/139‑1 (MQ=255)
                                                             tGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt                                                           <  2:420669/139‑1 (MQ=255)
                                                             tGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt                                                           <  2:34997/139‑1 (MQ=255)
                                                             tGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt                                                           <  2:738717/139‑1 (MQ=255)
                                                                acagtAACAAACCCGGGCTACAGCGTGGGCTGGAAGCTCGACATAACGAACTCATCGCCCT‑GGAGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATtggtgg                                                        <  1:514900/135‑1 (MQ=255)
                                                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGc                                                                                                               >  2:134803/1‑82 (MQ=255)
                                                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGc                                                                                                               <  1:134803/82‑1 (MQ=255)
                                                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc                                                      >  2:596584/1‑139 (MQ=255)
                                                                             agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGc                                           >  2:675136/1‑139 (MQ=255)
                                                                               agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCg                                         >  2:593452/1‑139 (MQ=255)
                                                                               agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATGGTTGGCCTATatcatcctgt                                         >  2:514900/1‑135 (MQ=255)
                                                                               agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gg‑gggggca                                                                                                                            <  2:134245/55‑4 (MQ=255)
                                                                               agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gg‑gggggca                                                                                                                            >  1:134245/1‑52 (MQ=255)
                                                                                 cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCggg                                       <  1:312036/139‑1 (MQ=255)
                                                                                           ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                                                            <  2:554939/76‑1 (MQ=255)
                                                                                           ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                                                            >  1:554939/1‑76 (MQ=255)
                                                                                           ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCt                             >  2:397825/1‑139 (MQ=255)
                                                                                              tGGAAGCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCtttt                          <  1:596584/139‑1 (MQ=255)
                                                                                                   gCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTtcatc                     <  2:273860/139‑1 (MQ=255)
                                                                                                   gCTCGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTtcatc                     <  2:772122/139‑1 (MQ=255)
                                                                                                      cGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGc                                                                                                               >  2:35584/1‑46 (MQ=255)
                                                                                                      cGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGc                                                                                                               <  1:612885/46‑1 (MQ=255)
                                                                                                      cGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGc                                                                                                               >  2:612885/1‑46 (MQ=255)
                                                                                                      cGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGc                                                                                                               <  1:35584/46‑1 (MQ=255)
                                                                                                      cGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggcc                                                                      <  2:799864/87‑1 (MQ=255)
                                                                                                      cGACATATCGAACTCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggcc                                                                      >  1:799864/1‑87 (MQ=255)
                                                                                                                   tCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                                                                           <  2:701215/53‑1 (MQ=255)
                                                                                                                   tCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                                                                           >  1:701215/1‑53 (MQ=255)
                                                                                                                   tCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgttg                                                         <  2:25624/87‑1 (MQ=255)
                                                                                                                   tCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgttg                                                         >  1:25624/1‑87 (MQ=255)
                                                                                                                   tCATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATTCGTTCAATGGGCg     >  2:682869/1‑139 (MQ=255)
                                                                                                                    cATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCCGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       >  2:125906/1‑72 (MQ=255)
                                                                                                                    cATCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCCGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       <  1:125906/72‑1 (MQ=255)
                                                                                                                     aTCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGaa   <  2:461789/139‑1 (MQ=255)
                                                                                                                     aTCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGaa   >  1:603593/1‑139 (MQ=255)
                                                                                                                      tCGCCCTGGG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGaaa  <  1:231304/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                                    
GATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GG‑GGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAA  >  NC_000913/3980827‑3981082

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