Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,225,816 A→T 9.2% intergenic (+332/+40) ycgI → / ← minE hypothetical protein/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,225,8160AT 9.2% 171.5 / 7.4 64intergenic (+332/+40)ycgI/minEhypothetical protein/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base T (4/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGT  >  AP009048/1225748‑1225885
                                                                    |                                                                     
agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                      <  1:5156204/70‑1 (MQ=255)
       gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                              <  1:1445424/71‑1 (MQ=255)
          aGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                             >  1:7360424/1‑69 (MQ=255)
          aGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                           <  1:6113811/71‑1 (MQ=255)
           gAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                           >  1:1258811/1‑70 (MQ=255)
             aTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                          >  1:3670320/1‑69 (MQ=255)
              tACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                       <  1:1044483/71‑1 (MQ=255)
                cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAACCTCAAGGCAGAGATAACTCTGCCTTGAAGa                                                      <  1:2978467/70‑1 (MQ=255)
                cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                     <  1:743882/71‑1 (MQ=255)
                cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                     <  1:4589808/71‑1 (MQ=255)
                   gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaaa                                                  >  1:2149230/1‑71 (MQ=255)
                    ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAct                                                                  >  1:638864/1‑54 (MQ=255)
                    ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAct                                                                  >  1:53215/1‑54 (MQ=255)
                    ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                 <  1:381470/71‑1 (MQ=255)
                     cgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                                >  1:68299/1‑71 (MQ=255)
                     cgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                                >  1:1889594/1‑71 (MQ=255)
                      gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActct                                                                >  1:3973679/1‑54 (MQ=255)
                      gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                               >  1:4348233/1‑71 (MQ=255)
                      gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                               >  1:7106393/1‑71 (MQ=255)
                       cTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                              >  1:3761247/1‑71 (MQ=255)
                        tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATaa                                                                    >  1:2266672/1‑48 (MQ=255)
                        tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                       <  1:125694/61‑1 (MQ=255)
                        tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcgc                                             >  1:4742371/1‑71 (MQ=255)
                         ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt                                            >  1:7038584/1‑71 (MQ=255)
                         ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt                                            >  1:4558819/1‑71 (MQ=255)
                             cacaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                        <  1:5500267/71‑1 (MQ=255)
                               caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt                                            >  1:6238972/1‑65 (MQ=255)
                               caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                        >  1:2550394/1‑69 (MQ=255)
                               caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                        >  1:6065406/1‑69 (MQ=255)
                                aaCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                     >  1:920132/1‑71 (MQ=255)
                                  cGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                                   >  1:1653794/1‑71 (MQ=255)
                                   gTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                  <  1:6767975/71‑1 (MQ=255)
                                   gTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                  <  1:7521284/71‑1 (MQ=255)
                                      gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTGTAAAAGCgcat                                <  1:2121482/70‑4 (MQ=255)
                                       aTGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                  <  1:5288757/67‑1 (MQ=255)
                                          agagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                            <  1:827758/70‑1 (MQ=255)
                                           gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAttt                          <  1:5369527/71‑1 (MQ=255)
                                             ggtaTGAAAAACTCAAGGCAGTGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                         <  1:3412021/67‑1 (MQ=255)
                                                aTGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                                    >  1:3203705/1‑56 (MQ=255)
                                                   aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAttt                          <  1:613258/63‑1 (MQ=255)
                                                   aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGctct                   <  1:5911357/70‑1 (MQ=255)
                                                   aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtctt                  <  1:7809269/71‑1 (MQ=255)
                                                    aaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                 <  1:3069975/71‑1 (MQ=255)
                                                      aaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct                >  1:3259003/1‑70 (MQ=255)
                                                       aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                  >  1:4951468/1‑51 (MQ=255)
                                                       aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg               >  1:1444189/1‑70 (MQ=255)
                                                           aaGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                              <  1:5975/35‑1 (MQ=255)
                                                           aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTGTGCTTcc          >  1:6997066/1‑71 (MQ=255)
                                                           aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc          >  1:5124286/1‑71 (MQ=255)
                                                            aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTAc                                       <  1:273651/41‑1 (MQ=21)
                                                            aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTAc                                       <  1:2004083/41‑1 (MQ=21)
                                                            aGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCg         >  1:662957/1‑71 (MQ=255)
                                                              gCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAg                                     >  1:5641261/1‑41 (MQ=25)
                                                              gCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAg                                     >  1:4130513/1‑41 (MQ=25)
                                                              gCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCgg        >  1:4177815/1‑70 (MQ=255)
                                                              gCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCgg        >  1:3557465/1‑70 (MQ=255)
                                                              gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                     <  1:7240866/41‑1 (MQ=255)
                                                              gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                     <  1:2119099/41‑1 (MQ=255)
                                                              gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                     <  1:2011373/41‑1 (MQ=255)
                                                              gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt       <  1:809162/71‑1 (MQ=255)
                                                              gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt       <  1:237033/71‑1 (MQ=255)
                                                                agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAtt                           >  1:968508/1‑49 (MQ=255)
                                                                agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAtt                           >  1:4836855/1‑49 (MQ=255)
                                                                agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt       >  1:3179636/1‑69 (MQ=255)
                                                                 gagaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg    <  1:361630/71‑1 (MQ=255)
                                                                  agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg   >  1:1165348/1‑71 (MQ=255)
                                                                    aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt       <  1:2711735/65‑1 (MQ=255)
                                                                    aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGt  <  1:1279513/70‑1 (MQ=255)
                                                                    |                                                                     
AGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGT  >  AP009048/1225748‑1225885

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