Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,273,734 G→A 7.3% intergenic (+15/‑96) aphA → / → yjbQ acid phosphatase/phosphotransferase, class B, non‑specific/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,273,7340GA 7.3% 214.2 / 5.4 68intergenic (+15/‑96)aphA/yjbQacid phosphatase/phosphotransferase, class B, non‑specific/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base A (3/2);  total (37/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 = 4.69e-01

CCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAT  >  W3110S.gb/4273664‑4273802
                                                                      |                                                                    
ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCg                                                                      <  1:7397882/71‑1 (MQ=255)
ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCg                                                                      <  1:1089344/71‑1 (MQ=255)
ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCg                                                                      <  1:6488976/71‑1 (MQ=255)
 cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                     <  1:4652089/71‑1 (MQ=255)
   tGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGa                                                                   >  1:974744/1‑71 (MQ=255)
    gCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTg                                                                    <  1:1169190/69‑1 (MQ=255)
    gCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt                                                                  >  1:1012264/1‑71 (MQ=255)
        caAGCGTGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGctc                                                              >  1:7273453/1‑71 (MQ=255)
        caAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGctc                                                              >  1:5630849/1‑71 (MQ=255)
        caAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGctc                                                              >  1:5315665/1‑71 (MQ=255)
        caAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGctc                                                              >  1:5275783/1‑71 (MQ=255)
        caAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGctc                                                              >  1:4370238/1‑71 (MQ=255)
         aaGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg                                                                 >  1:2509903/1‑67 (MQ=255)
         aaGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGctct                                                             <  1:3376027/71‑1 (MQ=255)
         aaGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGctct                                                             <  1:2342196/71‑1 (MQ=255)
           gCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGa                                                                   >  1:7234337/1‑63 (MQ=255)
             gggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCg                                                         <  1:6268709/71‑1 (MQ=255)
             gggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCg                                                         <  1:3224419/71‑1 (MQ=255)
              ggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCg                                                         <  1:3215250/70‑1 (MQ=255)
              ggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGa                                                        <  1:3611735/71‑1 (MQ=255)
               gTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                     <  1:4427304/57‑1 (MQ=255)
               gTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                     <  1:4759210/57‑1 (MQ=255)
               gTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATGATGCTCTCCCGaa                                                       <  1:5242993/71‑1 (MQ=255)
                  cGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGc                                                    >  1:2968568/1‑71 (MQ=255)
                   gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCt                                                   >  1:5107854/1‑71 (MQ=255)
                    tttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCt                                                   <  1:3183617/70‑1 (MQ=255)
                            agagGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttttttt                                           <  1:2503720/70‑1 (MQ=255)
                             gagGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTcc                                                           >  1:6984305/1‑53 (MQ=255)
                             gagGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTcc                                                           >  1:1081376/1‑53 (MQ=255)
                             gagGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTcc                                                           >  1:2078344/1‑53 (MQ=255)
                                gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                     <  1:1187125/40‑1 (MQ=255)
                                gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                     <  1:1644037/40‑1 (MQ=255)
                                gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                     <  1:4878263/40‑1 (MQ=255)
                                gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                     <  1:4251263/40‑1 (MQ=255)
                                gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAta                                       >  1:3869364/1‑69 (MQ=255)
                                gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                      >  1:4341210/1‑71 (MQ=255)
                                   aTCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTg                                                                    >  1:7475515/1‑38 (MQ=255)
                                   aTCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTCTTAATcagtt                                   <  1:4909571/71‑4 (MQ=255)
                                      gTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttt                                            <  1:809715/59‑1 (MQ=255)
                                       tCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCttt                               <  1:4011850/71‑1 (MQ=255)
                                        cAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttt                                               <  1:2894169/54‑1 (MQ=255)
                                        cAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTa                              >  1:2478256/1‑71 (MQ=255)
                                        cAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTccc                                                          >  1:381165/1‑43 (MQ=21)
                                            tCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTccc                                                          >  1:6401643/1‑39 (MQ=255)
                                             cAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAt                                        >  1:177871/1‑56 (MQ=255)
                                             cAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCtt                         >  1:7734626/1‑71 (MQ=255)
                                              aGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCttt                        <  1:4022454/71‑1 (MQ=255)
                                              aGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCtt                         <  1:4135758/70‑1 (MQ=255)
                                              aGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCtt                         <  1:6754993/70‑1 (MQ=255)
                                               gAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTc                       >  1:1053500/1‑71 (MQ=255)
                                                aaTACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCg                      >  1:7574915/1‑71 (MQ=255)
                                                       aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCtt               >  1:5052712/1‑71 (MQ=255)
                                                       aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCt                <  1:6357582/70‑1 (MQ=255)
                                                       aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCt                <  1:7527913/70‑1 (MQ=255)
                                                       aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCt                <  1:2957550/70‑1 (MQ=255)
                                                        cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCg                      >  1:954632/1‑63 (MQ=255)
                                                        cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtg              <  1:6929056/71‑1 (MQ=255)
                                                         agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACAcc                                  >  1:7555007/1‑50 (MQ=255)
                                                         agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgc             >  1:91390/1‑71 (MQ=255)
                                                         agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgc             >  1:7163036/1‑71 (MQ=255)
                                                           agCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctg           >  1:4520437/1‑71 (MQ=255)
                                                             cGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctg           >  1:1871767/1‑69 (MQ=255)
                                                                  gagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                   >  1:3474469/1‑40 (MQ=255)
                                                                  gagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                   >  1:6547985/1‑40 (MQ=255)
                                                                  gagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACt     >  1:3979077/1‑70 (MQ=255)
                                                                  gagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACt     >  1:3609283/1‑70 (MQ=255)
                                                                   agCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGa   <  1:3244398/71‑1 (MQ=255)
                                                                     cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAt  >  1:7282971/1‑70 (MQ=255)
                                                                     cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAt  <  1:1159188/70‑1 (MQ=255)
                                                                     cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAt  >  1:3192568/1‑70 (MQ=255)
                                                                     cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAt  <  1:3196488/70‑1 (MQ=255)
                                                                      gTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCTCTGTCTTGCTGCAAACt     <  1:6149732/66‑1 (MQ=255)
                                                                      |                                                                    
CCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAT  >  W3110S.gb/4273664‑4273802

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: