Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,764,566 T→A 8.3% intergenic (‑134/‑8) yfjM ← / → yfjN hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,764,5660TA 8.3% 160.0 / 5.1 59intergenic (‑134/‑8)yfjM/yfjNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base A (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTGAATC  >  W3110S.gb/2764499‑2764632
                                                                   |                                                                  
aaaTTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc                                                                   >  1:864360/1‑69 (MQ=255)
  aTTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAg                                                                >  1:348441/1‑70 (MQ=255)
  aTTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGAACCAgg                                                               <  1:1152222/71‑1 (MQ=255)
    ttAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGa                                                              <  1:3932157/70‑1 (MQ=255)
        cAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGa                                                         <  1:1330101/71‑1 (MQ=255)
         aaaTCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGAc                                                        <  1:349625/71‑1 (MQ=255)
               aTACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAg                                                                <  1:3324485/57‑1 (MQ=255)
               aTACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAg                                                  <  1:3947942/71‑1 (MQ=255)
                 acacTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAg                                                  >  1:3369940/1‑69 (MQ=255)
                   acTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGt                                              <  1:547738/71‑1 (MQ=255)
                           ggggTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaaaa                                      >  1:3947698/1‑71 (MQ=255)
                               tattatTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTa                                  <  1:2618766/71‑1 (MQ=255)
                                 ttattGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAg                                                  >  1:2469553/1‑53 (MQ=255)
                                 ttattGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGTAg                                               <  1:3319068/56‑1 (MQ=255)
                                 ttattGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaa                                        <  1:12627/63‑1 (MQ=255)
                                 ttattGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAAc                                     >  1:1433198/1‑66 (MQ=255)
                                 ttattGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAAc                                     >  1:4174494/1‑66 (MQ=255)
                                 ttattGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTaa                                 <  1:1164771/70‑1 (MQ=255)
                                   attGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGaaa                                                            >  1:2313313/1‑41 (MQ=255)
                                   attGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAAt                                                     >  1:1580071/1‑48 (MQ=255)
                                   attGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAAt                                                     >  1:3974638/1‑48 (MQ=255)
                                   attGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAg                                                  >  1:979994/1‑51 (MQ=255)
                                   attGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAg                                                  <  1:347621/51‑1 (MQ=255)
                                   attGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAg                                                  >  1:2035024/1‑51 (MQ=255)
                                   attGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGa                                                >  1:1222450/1‑53 (MQ=255)
                                   attGTAGAGTTTCCCCATATCTTTCTATGGGATCCAGGAAATGAc                                                        >  1:1104552/1‑45 (MQ=255)
                                      gTAGAGTTTCCCCATATGTTTCTATGGGAACCAGGAAATGACAATCAGGAGTTACAAAAACTTaaa                                >  1:578290/1‑66 (MQ=255)
                                       tAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGt                          <  1:4513440/71‑1 (MQ=255)
                                       tAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGt                          <  1:4408794/71‑1 (MQ=255)
                                        agagTTTCCCCATATGTTTCTATGGGAACCAGGAAATGACa                                                       >  1:2412382/1‑41 (MQ=255)
                                         gagTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGt                          >  1:1952995/1‑69 (MQ=255)
                                         gagTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGt                          >  1:812205/1‑69 (MQ=255)
                                            tttCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTa                                            <  1:787460/48‑1 (MQ=255)
                                            tttCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaa                                        <  1:4130907/52‑1 (MQ=255)
                                            tttCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTaa                                 <  1:3795336/59‑1 (MQ=255)
                                              tCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGc                    >  1:4544275/1‑70 (MQ=255)
                                              tCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGc                    >  1:2293240/1‑70 (MQ=255)
                                                 ccATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCaaa                 >  1:3385572/1‑70 (MQ=255)
                                                  cATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAg                                               >  1:1813652/1‑39 (MQ=255)
                                                   atatGTTTCTATGTGATCCAGGAAATGACAATCAGGAGTTACaaa                                        <  1:368846/45‑1 (MQ=255)
                                                   atatGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAatat              >  1:2851014/1‑71 (MQ=255)
                                                    tatGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaa                                        <  1:4608556/44‑1 (MQ=255)
                                                    tatGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCa                   <  1:2737262/65‑1 (MQ=255)
                                                    tatGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATc             <  1:2831346/71‑1 (MQ=255)
                                                     atGTTTCTATGGGAACCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATc                              <  1:2924340/53‑1 (MQ=255)
                                                      tGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt                                    <  1:3816455/46‑1 (MQ=255)
                                                      tGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATc             >  1:2959119/1‑69 (MQ=255)
                                                      tGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATc             >  1:2447009/1‑69 (MQ=255)
                                                      tGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATc             >  1:2306919/1‑69 (MQ=255)
                                                       gTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt                                    <  1:4062275/45‑1 (MQ=255)
                                                        tttCTATGGGAACCAGGAAATGACAATCAGGAGTTACAAAAACTTaa                                 <  1:2406130/47‑1 (MQ=255)
                                                         ttCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAAt                               >  1:3896845/1‑48 (MQ=255)
                                                          tCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTaa                                 <  1:3753236/45‑1 (MQ=255)
                                                          tCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCa                        >  1:1208576/1‑54 (MQ=255)
                                                          tCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCa                        >  1:2447735/1‑54 (MQ=255)
                                                          tCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACt       >  1:2159394/1‑71 (MQ=255)
                                                              tGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTg      <  1:674810/68‑1 (MQ=255)
                                                                ggATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTGAATc  <  1:1673269/70‑1 (MQ=255)
                                                                 gATCCAGGAAATGACAATCAGGAGTTACAAAAACTTaa                                 <  1:2672/38‑1 (MQ=255)
                                                                   |                                                                  
AAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTGAATC  >  W3110S.gb/2764499‑2764632

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: