Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,275,296 A→G 12.5% intergenic (‑120/‑206) ychN ← / → ychP conserved hypothetical protein/predicted invasin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,275,2960AG12.5% 77.0 / 5.0 32intergenic (‑120/‑206)ychN/ychPconserved hypothetical protein/predicted invasin
Reads supporting (aligned to +/- strand):  ref base A (16/12);  new base G (2/2);  total (18/14)
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.91e-01

ACCAATCGGTCACATTTTTAAGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAAGTTACCGACGTTTT  >  W3110S.gb/1275241‑1275362
                                                       |                                                                  
aCCAATCGGTCACATTTTTAAGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAGGTTGCCAGTAATCGt                                                     <  1:635997/71‑1 (MQ=255)
aCCAATCGGTCACATTTTTAAGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGt                                                     <  1:2510942/71‑1 (MQ=255)
     tCGGTCACATTTTTAAGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTc                                                >  1:3877469/1‑71 (MQ=255)
      cGGTCACATTTTTAAGGATTTTCCTGAACGCGCGAGAAAATACGACaaaa                                                                    >  1:2879008/1‑50 (MQ=255)
       ggTCACATTTTTAAGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCtt                                              >  1:3849790/1‑71 (MQ=255)
                tttAAGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTAt                                     >  1:1515115/1‑71 (MQ=255)
                tttAAGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTAt                                     >  1:3914557/1‑71 (MQ=255)
                    aGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCttt                                             <  1:511302/59‑1 (MQ=255)
                    aGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGt                                  <  1:189079/70‑1 (MQ=255)
                          ttCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTc                                                >  1:2808889/1‑50 (MQ=255)
                          ttCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCt                                               >  1:856996/1‑51 (MQ=255)
                          ttCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCt                                               >  1:3399192/1‑51 (MQ=255)
                          ttCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCt                                               >  1:182871/1‑51 (MQ=255)
                           tCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCat                           <  1:1964772/70‑1 (MQ=255)
                           tCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCat                           <  1:3481493/70‑1 (MQ=255)
                           tCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCat                           <  1:3445619/70‑1 (MQ=255)
                            ccTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCatt                          <  1:3887065/70‑1 (MQ=255)
                            ccTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCatt                          <  1:1370913/70‑1 (MQ=255)
                                aaaGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTa                                                   >  1:4011125/1‑41 (MQ=255)
                                aaaGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTa                                                   >  1:3035144/1‑41 (MQ=255)
                                  aGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGt                                                     <  1:1112240/37‑1 (MQ=255)
                                     gcgAGAAAATACGACAAAAGTTGCCAGTAATCGTTAtt                                                 >  1:2676659/1‑38 (MQ=255)
                                     gcgAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTa                                            >  1:3496119/1‑43 (MQ=255)
                                     gcgAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTa                                            >  1:1707994/1‑43 (MQ=255)
                                          aaaaTACGACAAAGGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGttt                                <  1:2254223/50‑1 (MQ=255)
                                             aTACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAAGTTACCGa        <  1:2262439/71‑1 (MQ=255)
                                             aTACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAAGTTACCGa        <  1:1650566/71‑1 (MQ=255)
                                               aCGACAAAGGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCatta                         >  1:1246233/1‑52 (MQ=255)
                                               aCGACAAAGGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCatt                          >  1:419156/1‑51 (MQ=255)
                                               aCGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAAGTTACCGACg      >  1:2227762/1‑71 (MQ=255)
                                               aCGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAAGTTACCGACg      >  1:1032849/1‑71 (MQ=255)
                                               aCGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAACTTACCGACg      >  1:1737605/1‑71 (MQ=255)
                                                   cAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAAGTTACCTACGtttt  <  1:2978309/71‑1 (MQ=255)
                                                      aaGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAAGTTAc           >  1:1970171/1‑59 (MQ=255)
                                                      aaGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAAGTTAc           >  1:1393741/1‑59 (MQ=255)
                                                       |                                                                  
ACCAATCGGTCACATTTTTAAGGATTTTCCTGAAAGCGCGAGAAAATACGACAAAAGTTGCCAGTAATCGTTATTCTTTAAGGCTATGGTTTTTCATTATTACCGGAAGTTACCGACGTTTT  >  W3110S.gb/1275241‑1275362

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: