Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,057,014 A→G 11.4% I121I (ATT→ATC yrhA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,057,0140AG11.4% 109.2 / 7.4 44I121I (ATT→ATCyrhAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/17);  new base G (3/2);  total (25/19)
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.49e-01

CAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTATGAAATCATCATTAT  >  W3110S.gb/4056945‑4057075
                                                                     |                                                             
cAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGaaa                                                              >  1:2701705/1‑71 (MQ=255)
cAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGaaa                                                              >  1:180572/1‑71 (MQ=255)
cAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGaaa                                                              >  1:1754145/1‑71 (MQ=255)
    aaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCt                                                          <  1:2536395/71‑1 (MQ=255)
    aaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCt                                                          <  1:1237476/71‑1 (MQ=255)
       acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCtata                                                       <  1:648959/71‑1 (MQ=255)
       acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCta                                                         <  1:1995717/69‑1 (MQ=255)
            catACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATccc                                                  >  1:1614935/1‑71 (MQ=255)
            catACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATccc                                                  >  1:2341300/1‑71 (MQ=255)
                    aaaTCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAGATGCTATAATCCCCGATCa                                            >  1:1564453/1‑69 (MQ=255)
                            aaGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCg                                     >  1:430889/1‑68 (MQ=255)
                             aGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCt                                 <  1:1416088/71‑1 (MQ=255)
                             aGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCt                                 <  1:1131049/71‑1 (MQ=255)
                               ttGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATaa                                                      >  1:654028/1‑48 (MQ=255)
                                tGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCta                                                         >  1:2394837/1‑44 (MQ=255)
                                tGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCta                                                         >  1:275610/1‑44 (MQ=255)
                                tGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCta                                                         >  1:2266849/1‑44 (MQ=255)
                                tGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGAt                                              >  1:2083787/1‑55 (MQ=255)
                                tGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTg                               >  1:1520742/1‑70 (MQ=255)
                                    agcaTCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTg                               >  1:860401/1‑66 (MQ=255)
                                       aTCACCTTTAATGTCATAAGTAAATATTGAAATGCtat                                                        <  1:2491659/38‑1 (MQ=255)
                                        tCACCTTTAATGTCATAAGTAAATATTGAAATGCtat                                                        <  1:1625690/37‑1 (MQ=255)
                                        tCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCg                                                <  1:228993/45‑1 (MQ=255)
                                        tCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAAcc                                      >  1:557141/1‑55 (MQ=255)
                                        tCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCt                                 <  1:2742071/60‑1 (MQ=255)
                                        tCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCt                                 <  1:414293/60‑1 (MQ=255)
                                        tCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCt                                 <  1:1082014/60‑1 (MQ=255)
                                          aCCTTTAATGTCATAAGTAAATATTGAAATGCtata                                                       >  1:2732240/1‑36 (MQ=255)
                                          aCCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTaa                                        >  1:1088625/1‑51 (MQ=255)
                                          aCCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAAc                                       >  1:479812/1‑52 (MQ=255)
                                               tAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTa                             >  1:1464608/1‑57 (MQ=255)
                                               tAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTAt                >  1:2063362/1‑70 (MQ=255)
                                               tAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTAt                >  1:1246990/1‑70 (MQ=255)
                                                aaTGTCATAAGTAAATATTGAGATGCTATAATCCCCGATCACTAACCGTTCTTg                               >  1:644685/1‑54 (MQ=255)
                                                aaTGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTg                               >  1:354037/1‑54 (MQ=255)
                                                  tGTCATAAGTAAATATTGAGATGCTATAATCCCCGAt                                              <  1:1643048/37‑1 (MQ=255)
                                                  tGTCATAAGTAAATATTGAGATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTATGaa             <  1:2223388/70‑1 (MQ=255)
                                                    tCATAAGTAAATATTGAGATGCTATAATCCCCGATCa                                            >  1:1819948/1‑37 (MQ=255)
                                                    tCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTATGAAatc          <  1:455031/71‑1 (MQ=255)
                                                      aTAAGTAAATATTGAAATGCTATAATCCCCGATCACTACCCGTTCTTGTAGATCAGGGTTTATGAAat           >  1:2414284/1‑68 (MQ=255)
                                                         aGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTATGAAATcatcat     >  1:1998210/1‑71 (MQ=255)
                                                         aGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTATGAAATcatcat     <  1:1830050/71‑1 (MQ=255)
                                                           tAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTATGAAATcatcat     <  1:84899/69‑1 (MQ=255)
                                                            aaaTATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTATGAAATCATCattat  >  1:287195/1‑71 (MQ=255)
                                                             aaTATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATc                         <  1:2321200/47‑1 (MQ=255)
                                                             aaTATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTATGAAATcatcat     >  1:868913/1‑67 (MQ=255)
                                                              atatTGAAATGCTATAATCCCCGATCACTAACCGTTCTTg                               <  1:1548228/40‑1 (MQ=255)
                                                                     |                                                             
CAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTTTAATGTCATAAGTAAATATTGAAATGCTATAATCCCCGATCACTAACCGTTCTTGTAGATCAGGGTTTATGAAATCATCATTAT  >  W3110S.gb/4056945‑4057075

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