Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,811,559 A→G 5.4% intergenic (‑89/‑191) yicE ← / → gltS predicted transporter/glutamate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,811,5590AG 5.4% 318.6 / 7.0 111intergenic (‑89/‑191)yicE/gltSpredicted transporter/glutamate transporter
Reads supporting (aligned to +/- strand):  ref base A (43/62);  new base G (3/3);  total (46/65)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 117 total aligned reads displayed.

GCAAAGCGGGCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGTTCA  >  AP009048/3811490‑3811627
                                                                     |                                                                    
gCAAAGCGGGCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTaa                                                                      <  1:5634923/70‑1 (MQ=255)
  aaaGCGGGCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAAt                                                                    <  1:85224/70‑1 (MQ=255)
  aaaGCGGGCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAAt                                                                    <  1:5180893/70‑1 (MQ=255)
   aaGCGGGCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATTTAAATATTTAGATga                                                                  <  1:5690791/71‑1 (MQ=255)
    aGCGGGCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATTTAAATATTTAGATgag                                                                 <  1:3844468/71‑1 (MQ=255)
    aGCGGGCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAAt                                                                    >  1:872571/1‑68 (MQ=255)
     gCGGGCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAAt                                                                    >  1:4441646/1‑67 (MQ=255)
        ggCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATaa                                                             <  1:794122/71‑1 (MQ=255)
           aTTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGAt                                                               <  1:1313405/66‑1 (MQ=255)
            ttttAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAAt                                                           >  1:4588579/1‑69 (MQ=255)
              ttAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTaaa                                                                     >  1:6850614/1‑57 (MQ=255)
              ttAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCAc                                                       >  1:27716/1‑71 (MQ=255)
              ttAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCAc                                                       >  1:2392714/1‑71 (MQ=255)
               tAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACt                                                      >  1:3267369/1‑71 (MQ=255)
                  cTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACtttt                                                   <  1:29647/71‑1 (MQ=255)
                  atGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACtttt                                                   <  1:4437269/70‑1 (MQ=255)
                    gACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATTTAAATATTTAGATgagg                                                                >  1:349819/1‑55 (MQ=37)
                     aCCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACtttttt                                                 >  1:4786131/1‑70 (MQ=255)
                     aCCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACtttttt                                                 <  1:1637742/70‑1 (MQ=255)
                     aCCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACtttttt                                                 <  1:6855041/70‑1 (MQ=255)
                     aCCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACttttt                                                  <  1:2084217/69‑1 (MQ=255)
                     aCCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTta                                                >  1:3583205/1‑71 (MQ=255)
                      ccGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTtat                                               <  1:682212/71‑1 (MQ=255)
                       cGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTtata                                              >  1:7590425/1‑71 (MQ=255)
                       cGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTtata                                              >  1:1369232/1‑71 (MQ=255)
                           agagAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAAc                                          >  1:40938/1‑71 (MQ=255)
                            gagaAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACt                                         <  1:6739864/71‑1 (MQ=255)
                            gagaAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACt                                         <  1:450284/71‑1 (MQ=255)
                             agaAAAGGAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGAt                                                               <  1:5584399/48‑1 (MQ=255)
                             agaAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAAt                                                                    <  1:7373077/43‑1 (MQ=255)
                             agaAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAAt                                                                    <  1:3766900/43‑1 (MQ=255)
                             agaAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGAt                                                               <  1:3410964/48‑1 (MQ=255)
                               aaaaGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATaa                                                             >  1:4586311/1‑48 (MQ=255)
                               aaaaGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACtttt                                      <  1:1608759/71‑1 (MQ=255)
                                 aaGCAAACGTTTGCCACATCCTCATCTAAATATTTaaa                                                                     <  1:1623639/38‑1 (MQ=255)
                                 aaGCAAACGTTTGCCACATCCTCATCTAAATATTTaaa                                                                     <  1:364691/38‑1 (MQ=255)
                                  aGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACttttt                                                  <  1:4123662/56‑1 (MQ=255)
                                   gCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTa                                  >  1:2462995/1‑71 (MQ=255)
                                    cAAACGTTTGCCACATCCTCATCTAAATATTTAAAt                                                                    <  1:3688609/36‑1 (MQ=255)
                                      aaCGTTTGCCACATCCTCATCTAAATATTTAAATGCGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:7090281/71‑1 (MQ=255)
                                      aaCGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:1931544/71‑1 (MQ=255)
                                      aaCGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:5183817/71‑1 (MQ=255)
                                      aaCGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:4962323/71‑1 (MQ=255)
                                      aaCGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:4367173/71‑1 (MQ=255)
                                       aCGTTTGCCACATCCTCATCTAAATATTTAAATGAGATa                                                              <  1:7456184/39‑1 (MQ=255)
                                       aCGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAtttt                              >  1:6464803/1‑71 (MQ=255)
                                       aCGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAtttt                              >  1:7631917/1‑71 (MQ=255)
                                        cGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTta                                                >  1:1141670/1‑52 (MQ=255)
                                        cGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACttt                                       <  1:2364790/61‑1 (MQ=255)
                                        cGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACttt                                       <  1:5808229/61‑1 (MQ=255)
                                          tttGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCt                           >  1:5947568/1‑71 (MQ=255)
                                            tGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGc                         >  1:2227630/1‑71 (MQ=255)
                                              ccACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:6236464/63‑1 (MQ=255)
                                              ccACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:7744287/63‑1 (MQ=255)
                                               cacaTCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACttt                                       <  1:5580297/54‑1 (MQ=255)
                                               cacaTCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCtt                       >  1:2561906/1‑70 (MQ=255)
                                                acctCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTa                      >  1:6723657/4‑70 (MQ=255)
                                                  aTCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTaaaa                   <  1:2213675/71‑1 (MQ=255)
                                                  aTCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTaaaa                   <  1:6216874/71‑1 (MQ=255)
                                                   tCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTta                                                >  1:1789673/1‑41 (MQ=255)
                                                   tCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTta                                                >  1:5992712/1‑41 (MQ=255)
                                                   tCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTta                                                >  1:3156231/1‑41 (MQ=255)
                                                   tCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACtttt                                      >  1:3525991/1‑51 (MQ=255)
                                                   tCCTCATCTAAATATTTAAATGAGATAAATGCAATTTTTta                                                >  1:6573971/1‑41 (MQ=255)
                                                     cTCATTTAAATATTTAGATGAGATAAATGCACTTTTTtatat                                             >  1:1520324/1‑42 (MQ=25)
                                                     cTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATc                 <  1:7094620/70‑1 (MQ=255)
                                                      tCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTa                                  <  1:7502930/52‑1 (MQ=255)
                                                        aTCTAAATATTTAAATGAGATAAATGCACTTTTTtat                                               <  1:3543388/37‑1 (MQ=255)
                                                        aTCTAAATATTTAAATGAGATAAATGCACTTTTTtat                                               >  1:2568558/1‑37 (MQ=255)
                                                        aTCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAtt                                >  1:1593414/1‑52 (MQ=255)
                                                        aTCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCa               <  1:494231/69‑1 (MQ=255)
                                                        aTCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCa               <  1:2604623/69‑1 (MQ=255)
                                                        aTCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCAt              <  1:816902/70‑1 (MQ=255)
                                                        aTCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATg             >  1:2183599/1‑71 (MQ=255)
                                                          cTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTaaaa                   <  1:5690549/63‑1 (MQ=255)
                                                          cTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGc            >  1:333324/1‑70 (MQ=255)
                                                           tAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCa          >  1:7383380/1‑71 (MQ=255)
                                                            aaaTATTTAAATGAGATAAATGCACTTTTTTATATAACttt                                       <  1:1061193/41‑1 (MQ=255)
                                                            aaaTATTTAAATGAGATAAATGCACTTTTTTATATAACttt                                       <  1:1047953/41‑1 (MQ=255)
                                                            aaaTATTTAAATGAGATAAATGCACTTTTTTATATAACttt                                       <  1:6538357/41‑1 (MQ=255)
                                                            aaaTATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCAt         >  1:3644596/1‑71 (MQ=255)
                                                            aaaTATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCAt         >  1:7538883/1‑71 (MQ=255)
                                                             aaTATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:2364927/48‑1 (MQ=255)
                                                             aaTATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCa          <  1:309846/69‑1 (MQ=255)
                                                             aaTATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCAtt        >  1:6740363/1‑71 (MQ=255)
                                                             aaTATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCAt         <  1:6829654/70‑1 (MQ=255)
                                                              atatTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:1578998/47‑1 (MQ=255)
                                                              atatTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttt                               <  1:2809548/47‑1 (MQ=255)
                                                               tatTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTg      <  1:5120981/71‑1 (MQ=255)
                                                                 tttAGATGAGATAAATGCACTTTTTTATATAACTTTTGGTAttttt                             >  1:6674992/1‑46 (MQ=255)
                                                                 tttAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGt     <  1:3611920/70‑1 (MQ=255)
                                                                 tttAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGt     <  1:3779574/70‑1 (MQ=255)
                                                                 tttAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGt     <  1:3183542/70‑1 (MQ=255)
                                                                 tttAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGt     <  1:1311097/70‑1 (MQ=255)
                                                                 tttAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGt     <  1:5678002/70‑1 (MQ=255)
                                                                 tttAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGt     <  1:6351045/70‑1 (MQ=255)
                                                                   tAAATGAGATAAATGCACTTTTTTATATAACTTTTg                                     >  1:7061704/1‑36 (MQ=255)
                                                                   tAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGTTc   <  1:5892911/70‑1 (MQ=255)
                                                                    aaaTGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGTTc   <  1:7818299/69‑1 (MQ=255)
                                                                     aaTGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGTTCa  <  1:7673729/69‑1 (MQ=255)
                                                                     |                                                                    
GCAAAGCGGGCATTTTAGCTGACCGAAAGAGAAAAGCAAACGTTTGCCACATCCTCATCTAAATATTTAAATGAGATAAATGCACTTTTTTATATAACTTTTGGTATTTTTCTGCTTAAAATCCATGCCATTTGTTCA  >  AP009048/3811490‑3811627

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