Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,628,826 C→T 8.4% *215* (TGA→TAA)  ytjB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,628,8260CT 8.4% 182.0 / 5.3 59*215* (TGA→TAA) ytjBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/31);  new base T (2/3);  total (25/34)
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.72e-01

CGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCT  >  AP009048/4628761‑4628895
                                                                 |                                                                     
tggaaGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCActct                                                                  <  1:1323987/68‑1 (MQ=255)
 gTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCActct                                                                  <  1:3602343/70‑1 (MQ=255)
 gTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCActct                                                                  >  1:2115380/1‑70 (MQ=255)
   aaGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTTACTCtttt                                                               >  1:1061739/1‑71 (MQ=255)
   aaGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttt                                                               >  1:2254113/1‑71 (MQ=255)
    agagTCAGAGATGAGCTGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttt                                                              <  1:2234133/71‑1 (MQ=255)
    agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTACTGGGTTACTCttttt                                                              >  1:2666783/1‑71 (MQ=255)
    agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttt                                                              <  1:3438335/71‑1 (MQ=255)
    agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttt                                                              <  1:1885416/71‑1 (MQ=255)
     gagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttt                                                              >  1:3498373/1‑70 (MQ=255)
       gTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtc                                                            <  1:109322/70‑1 (MQ=255)
       gTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtc                                                            >  1:4573244/1‑70 (MQ=255)
        tCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTc                                                                       <  1:4319793/58‑1 (MQ=255)
            agaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTTACTCTTTTTTCTCGCtt                                                      <  1:2579012/71‑1 (MQ=255)
            agaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt                                                      >  1:2286332/1‑71 (MQ=255)
             gaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCttt                                                     >  1:3268295/1‑71 (MQ=255)
             gaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt                                                      <  1:1702333/70‑1 (MQ=255)
              aTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTTACTCTTTTTTCTCGCttt                                                     <  1:1463777/70‑1 (MQ=255)
                    aGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTtcttcctc                                              >  1:1353986/1‑71 (MQ=255)
                         tAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCa                                                                      <  1:2312151/42‑1 (MQ=255)
                         tAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCg                                         <  1:2435347/71‑1 (MQ=255)
                         tAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCg                                         >  1:2405887/1‑71 (MQ=255)
                           aTGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGc                                       <  1:1611602/71‑1 (MQ=255)
                              tgGACATAACGATTTCCTTTCTTGTAGTAATGGGTTACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCAc                                     <  1:901840/70‑1 (MQ=255)
                               ggACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCAc                                     <  1:4083031/69‑1 (MQ=255)
                                     aaCGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGttt                               >  1:1074513/1‑69 (MQ=255)
                                     aaCGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTg                              >  1:14258/1‑70 (MQ=255)
                                     aaCGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGc                             >  1:2302650/1‑71 (MQ=255)
                                     aaCGATTTCCTTTCTTGTAGTAATGGGGCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGtt                                <  1:1287529/68‑1 (MQ=255)
                                      aCGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTg                              >  1:1492215/1‑69 (MQ=255)
                                       cGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtc                                                            >  1:4391227/1‑38 (MQ=255)
                                          tttCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCa                                      >  1:1850471/1‑57 (MQ=255)
                                           ttCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCgg                       >  1:1974938/1‑71 (MQ=255)
                                           ttCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCgg                       >  1:3725544/1‑71 (MQ=255)
                                            tCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTtcttcct                                               >  1:2629472/1‑46 (MQ=255)
                                            tCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCgg                       <  1:164324/70‑1 (MQ=255)
                                              cTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGtt                     <  1:2086382/70‑1 (MQ=255)
                                               tttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTttcttc                                                 <  1:3030870/41‑1 (MQ=255)
                                               tttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTttcttc                                                 <  1:624688/41‑1 (MQ=255)
                                               tttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGttt                               >  1:4036221/1‑59 (MQ=255)
                                                ttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTATTCCTCTTCCGGCACCGGTTTGCTGGCGGTTaa                   >  1:4300417/1‑70 (MQ=255)
                                                   ctttagtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAg                <  1:4376547/67‑1 (MQ=255)
                                                    tGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAgg               >  1:2647609/1‑70 (MQ=255)
                                                     gtagtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGc                                       >  1:4254644/1‑45 (MQ=255)
                                                     gtagtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGa              <  1:4011778/70‑1 (MQ=255)
                                                      tagtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGaa             <  1:4010166/70‑1 (MQ=255)
                                                        gtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAggg          <  1:3667081/71‑1 (MQ=255)
                                                        gtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAgg           <  1:4604924/70‑1 (MQ=255)
                                                         taATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCAc                                     <  1:4075858/43‑1 (MQ=255)
                                                         taATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGc         <  1:2810517/71‑1 (MQ=255)
                                                         taATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGc         <  1:2976132/71‑1 (MQ=255)
                                                         taATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGc         <  1:156571/71‑1 (MQ=255)
                                                           aTGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTgg                          <  1:3639678/52‑1 (MQ=255)
                                                              ggTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGtt                                >  1:1141158/1‑43 (MQ=255)
                                                              ggTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATtg    <  1:2518446/71‑1 (MQ=255)
                                                              ggTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATtg    >  1:3244787/1‑71 (MQ=255)
                                                               gTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTa                    <  1:3775445/54‑1 (MQ=255)
                                                               gTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATtgc   <  1:1726146/71‑1 (MQ=255)
                                                                tCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTaa                   <  1:2346432/54‑1 (MQ=255)
                                                                tCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATtgct  <  1:2816373/71‑1 (MQ=255)
                                                                 cACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATtgct  <  1:1355246/70‑1 (MQ=255)
                                                                 |                                                                     
CGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCT  >  AP009048/4628761‑4628895

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: