Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,290,774 A→G 32.0% intergenic (‑65/‑9) ychJ ← / → rssA conserved hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,290,7740AG32.0% 59.6 / 42.4 50intergenic (‑65/‑9)ychJ/rssAconserved hypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/26);  new base G (8/8);  total (16/34)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAATCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGAATAATGAGAAAGATAAA  >  W3110S.gb/1290711‑1290835
                                                               |                                                             
aataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGc                                                          <  1:3442609/69‑1 (MQ=255)
aataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGc                                                          <  1:970717/69‑1 (MQ=255)
aataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGc                                                          <  1:3583150/69‑1 (MQ=255)
aataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGc                                                          <  1:2405310/69‑1 (MQ=255)
aataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGc                                                          <  1:2319666/69‑1 (MQ=255)
aataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGc                                                          <  1:486117/69‑1 (MQ=255)
aataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGc                                                          <  1:2029009/69‑1 (MQ=255)
 ataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGc                                                          <  1:2180268/68‑1 (MQ=255)
 ataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCg                                                         <  1:120719/69‑1 (MQ=255)
 ataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCg                                                         <  1:2554810/69‑1 (MQ=255)
 ataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCg                                                         <  1:1527156/69‑1 (MQ=255)
 ataatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCg                                                         <  1:165464/69‑1 (MQ=255)
  taatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGa                                                        <  1:2592353/69‑1 (MQ=255)
  taatCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGa                                                        >  1:734992/1‑69 (MQ=255)
      ctcCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGAtgt                                                     <  1:2030906/68‑1 (MQ=255)
           tGGAAACAATACCGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGt                                                             >  1:2587936/1‑55 (MQ=255)
           tGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGt                                                             >  1:2861914/1‑55 (MQ=255)
           tGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGt                                                             >  1:407220/1‑55 (MQ=255)
           tGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTAc                                               <  1:712138/69‑1 (MQ=255)
           tGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTAc                                               <  1:750520/69‑1 (MQ=255)
            tgAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACg                                              <  1:2239057/68‑1 (MQ=255)
               aaCAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGa                                             <  1:1866139/67‑1 (MQ=255)
                   aTAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCAt                                       >  1:797350/1‑69 (MQ=255)
                       cGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCa                                   >  1:1086047/1‑69 (MQ=255)
                              ttAACCGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCATTc                                     >  1:451696/1‑60 (MQ=255)
                              ttAACCGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCATTCCAGGGGAAt                            >  1:3133100/1‑69 (MQ=255)
                              ttAACCGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCATTCCAGGGGAAt                            >  1:139881/1‑69 (MQ=255)
                              ttAACCGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCATTCCAGGGGAAt                            >  1:853960/1‑69 (MQ=255)
                              ttAACCGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCATTCCAGGGGAAt                            >  1:1091163/1‑69 (MQ=255)
                              ttAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTAc                                               <  1:986529/50‑1 (MQ=255)
                              ttAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTAc                                               <  1:2231283/50‑1 (MQ=255)
                                aaCCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGa                                                        >  1:1390712/1‑39 (MQ=255)
                                 aCCGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCtt                         >  1:2910447/1‑69 (MQ=255)
                                  ccGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAAt                            <  1:2084337/65‑1 (MQ=255)
                                   cGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCg                                                         <  1:3084652/35‑1 (MQ=255)
                                   cGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCg                                                         <  1:1240798/35‑1 (MQ=255)
                                   cGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCg                                                         <  1:1360956/35‑1 (MQ=255)
                                   cGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCg                                                         <  1:547611/35‑1 (MQ=255)
                                   cGCCTGAGTAGCACTATGTTAACCGAGCGGTAGCg                                                         <  1:2111712/35‑1 (MQ=255)
                                    gCCTTAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCg                      <  1:118315/69‑1 (MQ=255)
                                    gCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCg                      <  1:550510/69‑1 (MQ=255)
                                           tAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGaata                >  1:1449532/1‑68 (MQ=255)
                                           tAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGaata                <  1:282478/68‑1 (MQ=255)
                                            aGCACTATGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGaat                 <  1:891228/66‑1 (MQ=255)
                                            aGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGaataat              <  1:64362/69‑1 (MQ=255)
                                                   tGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGAATAATg             >  1:687624/1‑63 (MQ=255)
                                                   tGTTAACCGAGCGGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGAATAATg             <  1:3612433/63‑1 (MQ=255)
                                                     ttAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGAATAATGAGAAAGa      >  1:2388041/1‑68 (MQ=255)
                                                         ccGAGCGGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGAATAATGAGAAAGATaaa  <  1:569472/68‑1 (MQ=255)
                                                         ccGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGAATAATGAGAAAGATaaa  <  1:968759/68‑1 (MQ=255)
                                                               |                                                             
AATAATCTCCCTGGAAACAATAACGGCGTATTAACCGCCTGAGTAGCACTATGTTAACCGAGCAGTAGCGATGTGGCTACGATTGCATTCCAGGGGAATCTTGCGGGAATAATGAGAAAGATAAA  >  W3110S.gb/1290711‑1290835

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: