Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,436,722 T→C 17.4% intergenic (‑751/+177) ynaE ← / ← uspF predicted DNA‑binding transcriptional regulator/stress‑induced protein, ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,436,7220TC17.4% 108.5 / 15.9 46intergenic (‑751/+177)ynaE/uspFpredicted DNA‑binding transcriptional regulator/stress‑induced protein, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (11/27);  new base C (4/4);  total (15/31)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGT  >  W3110S.gb/1436664‑1436767
                                                          |                                             
tGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACACTTAAAAGCGAc                                    <  1:1888213/70‑1 (MQ=255)
       aCTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAAtt                              >  1:189669/1‑69 (MQ=255)
       aCTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCg                            >  1:961991/1‑71 (MQ=255)
       aCTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCg                            >  1:812206/1‑71 (MQ=255)
       aCTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCg                            >  1:537296/1‑71 (MQ=255)
       aCTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCg                            >  1:479607/1‑71 (MQ=255)
       aCTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCg                            >  1:406647/1‑71 (MQ=255)
       aCTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCg                            >  1:1804262/1‑71 (MQ=255)
       aCTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCg                            >  1:1752166/1‑71 (MQ=255)
        cTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                           <  1:1133486/71‑1 (MQ=255)
        cTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                           <  1:1159826/71‑1 (MQ=255)
        cTTGTTCAGGTTGTATTGTCCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                           <  1:1053724/71‑1 (MQ=255)
           gTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTta                         <  1:996598/70‑1 (MQ=255)
           gTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTta                         <  1:846694/70‑1 (MQ=255)
           gTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTta                         <  1:339963/70‑1 (MQ=255)
           gTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTta                         <  1:1262861/70‑1 (MQ=255)
           gTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTta                         <  1:1294098/70‑1 (MQ=255)
           gTTCAGGTGGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTta                         <  1:1832704/70‑1 (MQ=255)
            ttCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACACTTAAAAGCGACTCAATTCGTtat                        <  1:164650/70‑1 (MQ=255)
            ttCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACACTTAAAAGCGACTCAATTCGTtat                        <  1:985759/70‑1 (MQ=255)
            ttCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACACTTAAAAGCGACTCAATTCGTtat                        <  1:389962/70‑1 (MQ=255)
                   tGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                           >  1:834082/1‑60 (MQ=255)
                   tGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                           >  1:1722237/1‑60 (MQ=255)
                   tGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                           >  1:453384/1‑60 (MQ=255)
                        tGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATgg                    <  1:1332907/62‑1 (MQ=255)
                        tGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATgg                    <  1:1512058/62‑1 (MQ=255)
                          ttctttctTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAagag         <  1:851651/71‑1 (MQ=255)
                          ttctttctTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAagag         <  1:467531/71‑1 (MQ=255)
                          ttctttctTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAagag         <  1:1435937/71‑1 (MQ=255)
                          ttctttctTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAagag         <  1:1458081/71‑1 (MQ=255)
                          ttctttctTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAagag         <  1:1545452/71‑1 (MQ=255)
                           tctttcttACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAagag         <  1:1567116/70‑1 (MQ=255)
                               tcttACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTa       <  1:425588/68‑1 (MQ=255)
                                cttACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGc    <  1:697222/70‑1 (MQ=255)
                                cttACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGc    <  1:128632/70‑1 (MQ=255)
                                  tACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGt  <  1:674965/70‑1 (MQ=255)
                                   aCTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGt  <  1:1736216/69‑1 (MQ=255)
                                    cTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGt  <  1:1546901/68‑1 (MQ=255)
                                     tAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTa       <  1:1549931/62‑1 (MQ=255)
                                       aTTTCTTGATTTTGCGACACTTAAAAGCGACTCAATTCg                            >  1:237030/1‑39 (MQ=255)
                                       aTTTCTTGATTTTGCGACACTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTa       >  1:1133605/1‑60 (MQ=255)
                                       aTTTCTTGATTTTGCGACACTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTa       >  1:1621183/1‑60 (MQ=255)
                                       aTTTCTTGATTTTGCGACACTTAAAAGCGACTCAATTCATTATATGGCATCAGAAGAGTa       >  1:1101117/1‑60 (MQ=255)
                                        tttCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                           <  1:781152/39‑1 (MQ=255)
                                             tGATTTTGCGACATTTAAAAGCGACTCAATTCGTta                         <  1:1576756/36‑1 (MQ=255)
                                                         aTTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGt  <  1:1770425/47‑1 (MQ=255)
                                                          |                                             
TGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGT  >  W3110S.gb/1436664‑1436767

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: