Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,708,410 T→A intergenic (‑705/+206) dppA ← / ← proK dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,708,4100TA95.2% 195.0 / ‑1.1 62intergenic (‑705/+206)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (30/29);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGG  >  NC_000913/3708313‑3708505
                                                                                                 |                                                                                               
ggCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGc                                                                                               <  2:1648096/100‑1 (MQ=255)
ggCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGc                                                                                               <  1:1771251/100‑1 (MQ=255)
 gCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCa                                                                                              <  2:710143/100‑1 (MQ=255)
   acaTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCaaa                                                                                            >  2:1220603/1‑100 (MQ=255)
    caTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCaaaa                                                                                           >  2:1028505/1‑100 (MQ=255)
        tGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGc                                                                                       <  1:1220600/100‑1 (MQ=255)
         gcgtGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCGCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGcc                                                                                      <  2:564478/97‑1 (MQ=255)
         gCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGcc                                                                                      <  1:689913/100‑1 (MQ=255)
         gCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGcc                                                                                      <  2:570855/100‑1 (MQ=255)
          ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAACCaaaaagg                                                                                        >  1:1914682/1‑95 (MQ=255)
           cTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTg                                                                                    >  1:813390/1‑100 (MQ=255)
           cTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTg                                                                                    >  2:1514750/1‑100 (MQ=255)
               tGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCg                                                                                >  2:779618/1‑100 (MQ=255)
                gCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGt                                                                               >  2:636429/1‑100 (MQ=255)
                 cGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGATGAAAAAAAAAAAAAGCCTGCTCGtt                                                                              <  1:1868587/100‑1 (MQ=18)
                  gACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTg                                                                             <  2:148049/100‑1 (MQ=255)
                     gCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGc                                                                          <  2:854642/100‑1 (MQ=255)
                       tGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAg                                                                        <  1:273522/100‑1 (MQ=255)
                        gACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAgg                                                                       >  2:1066684/1‑100 (MQ=255)
                         aCGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGc                                                                      <  2:283735/100‑1 (MQ=255)
                          cgcgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCAACGAGCAGGCt                                                                     <  2:1921645/100‑1 (MQ=17)
                            cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCttt                                                                   <  2:49796/100‑1 (MQ=255)
                            cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCttt                                                                   >  1:56085/1‑100 (MQ=255)
                            cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCttt                                                                   <  2:792010/100‑1 (MQ=255)
                            cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCttt                                                                   <  1:1196072/100‑1 (MQ=255)
                            cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCAACGAGCAGGCttt                                                                   <  2:1695833/100‑1 (MQ=17)
                                ttATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGa                                                               <  2:1516786/100‑1 (MQ=255)
                                ttATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCAACGAGCAGGCtttttt                                                                <  2:112297/99‑3 (MQ=14)
                                 tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGaa                                                              >  1:1275111/1‑100 (MQ=255)
                                   tCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAtt                                                            >  1:1665596/1‑100 (MQ=255)
                                    cAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtttagatagga                                                           >  1:1072128/1‑92 (MQ=255)
                                      ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTgg                                                         >  1:546036/1‑100 (MQ=255)
                                      ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTgg                                                         >  1:119796/1‑100 (MQ=255)
                                        ccTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTc                                                       <  2:546038/100‑1 (MQ=25)
                                          taggtcgagtcgaACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCgg                                                     <  1:1686485/97‑1 (MQ=11)
                                            cGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTg                                                   <  2:1903537/100‑1 (MQ=14)
                                                        ccTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTAc                                       >  2:1777604/1‑100 (MQ=14)
                                                            ccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGc                                   >  2:1652871/1‑100 (MQ=25)
                                                               gCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCAc                                >  2:1586910/1‑100 (MQ=255)
                                                                  cccGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGCGCAGGCTTTTCGAATTTGGTCGGTGATAGTGGATTACTCGCCACTTc                             >  1:50280/1‑100 (MQ=255)
                                                                    cGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGt                           <  2:1760254/100‑1 (MQ=255)
                                                                      gAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTgg                         >  2:1125970/1‑100 (MQ=255)
                                                                       aGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGc                        <  1:1585286/100‑1 (MQ=255)
                                                                        ggTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCt                       >  1:888257/1‑100 (MQ=255)
                                                                          ttCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCg                     <  2:1334746/100‑1 (MQ=255)
                                                                          ttCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCg                     >  2:71318/1‑100 (MQ=255)
                                                                              cATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCAGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCt                 >  1:1307324/1‑100 (MQ=255)
                                                                                 ccTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGTCCCGCCCTg                >  1:653368/1‑98 (MQ=255)
                                                                                   tCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCggg            <  2:970041/100‑1 (MQ=255)
                                                                                    cATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCGTGCGGGc           >  1:1047700/1‑100 (MQ=255)
                                                                                    cATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGc           >  1:1050730/1‑100 (MQ=255)
                                                                                    cATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGACGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCggg            >  1:1709075/1‑99 (MQ=255)
                                                                                     aTGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGcc          >  2:1345378/1‑100 (MQ=255)
                                                                                       gAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGt        >  2:1862363/1‑100 (MQ=255)
                                                                                        aGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGtt       <  2:736789/100‑1 (MQ=255)
                                                                                         gCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTg      <  1:1153353/100‑1 (MQ=255)
                                                                                            gCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTg   >  2:29263/1‑100 (MQ=255)
                                                                                            gCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTg   >  2:1825969/1‑100 (MQ=255)
                                                                                            gCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTg   <  1:636428/100‑1 (MQ=255)
                                                                                             cAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTGGCCCTGCGGGCCGTTGCt    >  1:1208544/1‑98 (MQ=255)
                                                                                             cAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTgg  <  2:1355806/100‑1 (MQ=255)
                                                                                             cAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTgg  <  2:1303650/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGG  >  NC_000913/3708313‑3708505

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