Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,561,090 G→A 55.6% intergenic (‑88/‑585) yjiC ← / → yjiD hypothetical protein/DNA replication/recombination/repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,561,0900GA55.6% 4.2 / 33.0 27intergenic (‑88/‑585)yjiC/yjiDhypothetical protein/DNA replication/recombination/repair protein
Reads supporting (aligned to +/- strand):  ref base G (7/5);  new base A (11/4);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAGGGTAGGAAACACTCTA  >  W3110S.gb/4561022‑4561152
                                                                    |                                                              
tATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatgca                                                              >  1:74824/1‑71 (MQ=255)
    tCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGacatacat                                                             >  1:759344/1‑68 (MQ=255)
    tCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGtt                                                          >  1:721856/1‑71 (MQ=255)
    tCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATACATGtt                                                          >  1:1050426/1‑71 (MQ=255)
       cTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATACATGttt                                                         <  1:186686/69‑1 (MQ=255)
        tGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATACATGTTTCtg                                                      >  1:1001368/1‑71 (MQ=255)
          aTGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTTCtg                                                      <  1:1345779/69‑1 (MQ=255)
           tGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATACATGTTTCtgttg                                                   <  1:960084/71‑1 (MQ=255)
                aGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATACATGTTTCTGTTGTAGt                                               <  1:1568637/70‑1 (MQ=255)
                         gATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATACATGTTTCTGTTGTAGTGATATCAAct                                     >  1:155896/1‑71 (MQ=255)
                              aTTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTTCtgtt                                                    >  1:1497135/1‑51 (MQ=255)
                                       tagtaAAAAGAGAAGTAAACAAATGACATACATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGatat                       >  1:194602/1‑71 (MQ=255)
                                       tagtaAAAAGAGAAGTAAACAAATGACATACATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGatat                       >  1:58653/1‑71 (MQ=255)
                                       tagtaAAAAGAGAAGTAAACAAATGACATACATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGatat                       >  1:375595/1‑71 (MQ=255)
                                          taAAAAGAGAAGTAAACAAATGACATACATGTTTCtgttg                                                   <  1:264167/40‑1 (MQ=255)
                                          taAAAAGAGAAGTAAACAAATGACATACATGTTTCtgtt                                                    >  1:917602/1‑39 (MQ=255)
                                               agagAAGTAAACAAATGACATGCATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTa                     <  1:677223/65‑1 (MQ=255)
                                                 agaAGTAAACAAATGACATACATGTTTCTGTTGTa                                                 >  1:1732428/1‑35 (MQ=255)
                                                      tAAACAAATGACATGCATGTTTCTGTTGTAGTGATATCAActc                                    <  1:1168992/43‑1 (MQ=255)
                                                        aaCAAATGACATGCATGTTTCTGTTGTAGTGGTATCAACTCTACACGGTGATATTAAAGGGTAGGAAaca       <  1:316373/70‑1 (MQ=255)
                                                        aaCAAATGACATACATGTTTCTGTTGTAGTGATATCAACTCTACACgg                             >  1:1474406/1‑48 (MQ=255)
                                                        aaCAAATGACATACATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGa                          >  1:1631972/1‑51 (MQ=255)
                                                            aaTGACATGCATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAGGGTAGGAAACACTCTa  <  1:1086334/71‑1 (MQ=255)
                                                            aaTGACATACATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAGGGTAGGAGACActct   >  1:618233/1‑70 (MQ=255)
                                                               gACATGCATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAgg                 >  1:1381696/1‑53 (MQ=255)
                                                               gACATGCATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAgg                 >  1:674525/1‑53 (MQ=255)
                                                               gACATGCATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAgg                 >  1:1735292/1‑53 (MQ=255)
                                                               gACATGCATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAgg                 >  1:1683905/1‑53 (MQ=255)
                                                               gACATGCATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAgg                 >  1:749963/1‑53 (MQ=255)
                                                                  atgcatgTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAGGGt               >  1:1597621/1‑52 (MQ=255)
                                                                   tgcatgTTTCTGTTGTAGTGATATCAACTCTACACGGTGata                        >  1:770884/1‑42 (MQ=255)
                                                                    |                                                              
TATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTTCTGTTGTAGTGATATCAACTCTACACGGTGATATTAAAGGGTAGGAAACACTCTA  >  W3110S.gb/4561022‑4561152

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: