Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,427,892 A→T 9.8% intergenic (+198/‑276) ydaY → / → ynaA hypothetical protein/predicted tail protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,427,8920AT 9.8% 188.8 / 9.3 71intergenic (+198/‑276)ydaY/ynaAhypothetical protein/predicted tail protein
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base T (4/3);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCT  >  AP009048/1427825‑1427961
                                                                   |                                                                     
cTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACa                                                                     <  1:5450799/70‑1 (MQ=255)
   ttCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACa                                                                     >  1:1827768/1‑67 (MQ=255)
    tCGGTAATATCCGGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAgaga                                                                 >  1:7103914/1‑70 (MQ=255)
    tCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAgaga                                                                 >  1:351150/1‑70 (MQ=255)
    tCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAgaga                                                                 >  1:2439889/1‑70 (MQ=255)
     cGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGAt                                                                >  1:2751175/1‑70 (MQ=255)
      ggTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaa                                                              <  1:568896/71‑1 (MQ=255)
      ggTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaa                                                              <  1:4571126/71‑1 (MQ=255)
      ggTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaa                                                              <  1:5294115/71‑1 (MQ=255)
      ggTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaa                                                              <  1:2214469/71‑1 (MQ=255)
        tAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAt                                                            <  1:5175586/71‑1 (MQ=255)
        tAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAt                                                            <  1:4651089/71‑1 (MQ=255)
        tAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAt                                                            >  1:1256091/1‑71 (MQ=255)
         aaTATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAt                                                            >  1:7352648/1‑70 (MQ=255)
          atatCCCGTGCCCGACGTTATGCTTGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATg                                                           >  1:7523764/1‑70 (MQ=255)
          atatCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaa                                                              <  1:3695499/67‑1 (MQ=255)
             tCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGaa                                                         >  1:320798/1‑69 (MQ=255)
             tCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAat                                                        <  1:2205598/70‑1 (MQ=255)
                cGTGCCTGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAGTGAATAttt                                                    <  1:7312190/71‑1 (MQ=255)
                  tGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTa                                                   <  1:607103/70‑1 (MQ=255)
                    cacGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGcc                                                <  1:4339248/69‑1 (MQ=255)
                     ccGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGcc                                                <  1:4864369/70‑1 (MQ=255)
                      cGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAAAGAGATAAATGAATATTTAGCCAc                                              >  1:7258281/1‑71 (MQ=255)
                       gACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                             >  1:2753790/1‑55 (MQ=255)
                       gACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAtt                                                     >  1:1109619/1‑63 (MQ=255)
                       gACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAtt                                                     >  1:3864689/1‑63 (MQ=255)
                       gACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAg                                                  >  1:3887715/1‑66 (MQ=255)
                       gACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACt                                             >  1:5575264/1‑71 (MQ=255)
                        aCGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAc                                                                      <  1:4122113/45‑1 (MQ=255)
                        aCGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACt                                             >  1:2564239/1‑70 (MQ=255)
                         cGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAg                                                  >  1:1376680/1‑64 (MQ=255)
                             aTGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCAtgag                                       >  1:6447463/1‑68 (MQ=255)
                               gCTGGAATGGCTGGGGTGCCTTTGCATTTATCTCTGTCAGAGATa                                                               >  1:7160863/1‑45 (MQ=25)
                               gCTGGAATGGCTGGGGTGCCTTTGCATTTATCTCTGTCAGAGATa                                                               >  1:4655619/1‑45 (MQ=25)
                                 tGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCgg                                      <  1:1233705/68‑1 (MQ=255)
                                     aTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgat                               >  1:6941771/1‑71 (MQ=255)
                                      tGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgatt                              <  1:6336350/71‑1 (MQ=255)
                                      tGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgatt                              <  1:7613150/71‑1 (MQ=255)
                                        gCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTAtttattga                            <  1:4876748/71‑1 (MQ=255)
                                            gggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAgcg                         <  1:6373869/70‑1 (MQ=255)
                                             ggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGcgcg                       <  1:6474881/71‑1 (MQ=255)
                                               tgcctttgccttTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattg                                 <  1:7722678/59‑1 (MQ=255)
                                               tgcctttgccttTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattg                                 <  1:3578179/59‑1 (MQ=255)
                                                 cctttgcctttATCTCTGACAGAGATAAATGAATAttt                                                    >  1:4361293/1‑38 (MQ=255)
                                                 cctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATg                    <  1:1705239/70‑1 (MQ=255)
                                                  ctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGaa                  <  1:48703/71‑1 (MQ=255)
                                                  ctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGaa                  <  1:6536552/71‑1 (MQ=255)
                                                   tttgcctttATCTCTGACAGAGATAAATGAATATTTAGc                                                 <  1:3993025/39‑1 (MQ=255)
                                                    ttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAtt                >  1:7854137/1‑71 (MQ=255)
                                                     tgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGaa                  <  1:2378749/68‑1 (MQ=255)
                                                     tgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAttt               >  1:724205/1‑71 (MQ=255)
                                                      gcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGaa                  <  1:1004649/67‑1 (MQ=255)
                                                      gcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAttt               >  1:6284217/1‑70 (MQ=255)
                                                      gcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAttt               <  1:6028357/70‑1 (MQ=255)
                                                      gcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAttt               <  1:2270465/70‑1 (MQ=255)
                                                      gcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGAGATGAAttt               <  1:1255943/70‑1 (MQ=255)
                                                        ctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGcgcg                       <  1:5762720/60‑1 (MQ=255)
                                                        atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCCTCCGGTATTGATTGAg                           >  1:6764980/2‑56 (MQ=37)
                                                        atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTg              >  1:2937126/2‑69 (MQ=255)
                                                         tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGt                                     <  1:3074417/45‑1 (MQ=255)
                                                          ttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCAt                                          <  1:2383875/39‑1 (MQ=255)
                                                          ttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTat                                   >  1:7750368/1‑46 (MQ=255)
                                                           tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAt                 >  1:1035618/1‑63 (MQ=255)
                                                           tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAt                 >  1:709347/1‑63 (MQ=255)
                                                            aTCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgatt                              <  1:4165019/49‑1 (MQ=255)
                                                            aTCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgatt                              <  1:2816280/49‑1 (MQ=255)
                                                               tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAt                 <  1:3731584/59‑1 (MQ=255)
                                                                 tGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAttt               <  1:811630/59‑1 (MQ=255)
                                                                 tGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATc   >  1:7232358/1‑71 (MQ=255)
                                                                 tGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATc   >  1:6243174/1‑71 (MQ=255)
                                                                  gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCt  >  1:4574316/1‑71 (MQ=255)
                                                                   aCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCt  <  1:3245365/70‑1 (MQ=255)
                                                                   |                                                                     
CTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCT  >  AP009048/1427825‑1427961

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: