Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,620,710 C→T 87.3% intergenic (‑88/‑124) marC ← / → marR predicted transporter/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,620,7100CT87.3% 145.5 / 13.0 55intergenic (‑88/‑124)marC/marRpredicted transporter/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base T (35/13);  total (37/18)
Fisher's exact test for biased strand distribution p-value = 3.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TATCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGCATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAATTCATTCATTTGACT  >  W3110S.gb/1620642‑1620776
                                                                    |                                                                  
tATCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTAt                                                                  >  1:2809321/1‑71 (MQ=255)
tATCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTAt                                                                  >  1:1860592/1‑71 (MQ=255)
tATCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTAt                                                                  >  1:1398310/1‑71 (MQ=255)
tATCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTAt                                                                  >  1:873481/1‑71 (MQ=255)
tATCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGCAt                                                                  >  1:111018/1‑71 (MQ=255)
  tCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATct                                                                >  1:1841489/1‑70 (MQ=255)
  tCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCg                                                                >  1:2687191/1‑71 (MQ=255)
  tCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCg                                                                >  1:927426/1‑71 (MQ=255)
        ttATCAGCGAAAAGTATAAGGGGTAAACAATGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCGCATTg                                                           <  1:2304413/70‑1 (MQ=255)
        ttATCAGCGAAAAGTATAAGGGGTAAAAAATGATAAAGTGTCACTCTTTAGCTAGCCTTGTATc                                                                 >  1:1064207/1‑64 (MQ=255)
                      tataAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGCATCGCATTGAACAAAACTTGAAcc                                            <  1:609724/71‑1 (MQ=255)
                      tataAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGCATCGCATTGAACAAAACTTGAAcc                                            <  1:1411893/71‑1 (MQ=255)
                      tataAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGCATCGCATTGAACAAAACTTGAAcc                                            <  1:2047310/71‑1 (MQ=255)
                          aGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGCATCGCATTGAACAAAACTTGAACCGAtt                                        >  1:2162513/1‑71 (MQ=255)
                            gggTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTg                                                <  1:125918/61‑1 (MQ=255)
                               tAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        <  1:2980475/66‑1 (MQ=255)
                               tAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        <  1:2609094/66‑1 (MQ=255)
                               tAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        <  1:3436617/66‑1 (MQ=255)
                               tAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        <  1:788686/66‑1 (MQ=255)
                               tAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        <  1:1541791/66‑1 (MQ=255)
                                 aaCAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCGCATTGAACaa                                                      <  1:2227196/50‑1 (MQ=255)
                                 aaCAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGTATCGCATTGAACaa                                                      <  1:3332574/50‑1 (MQ=255)
                                              tgtCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGt                    >  1:2001613/1‑71 (MQ=255)
                                              tgtCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGt                    >  1:969064/1‑71 (MQ=255)
                                              tgtCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGt                    >  1:1128618/1‑71 (MQ=255)
                                               gtCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGt                    >  1:3249205/1‑70 (MQ=255)
                                               gtCACTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGt                    >  1:1526698/1‑70 (MQ=255)
                                                tCACTCTTTAGCTAGCCTTGCATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTc                   <  1:2894340/70‑1 (MQ=255)
                                                  aCTCTTTAGCTAGCCTTGTATCGCATTGAACAAAAc                                                   >  1:3514254/1‑36 (MQ=255)
                                                  aCTCTTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAat                >  1:397491/1‑71 (MQ=255)
                                                  aCTCTTTAGCTAGCCTTGTATAGCATTGAACAAAAc                                                   >  1:1074532/1‑36 (MQ=38)
                                                    tctTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAatt               >  1:813872/1‑70 (MQ=255)
                                                    tctTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAatt               >  1:1549193/1‑70 (MQ=255)
                                                    tctTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAatt               >  1:1952347/1‑70 (MQ=255)
                                                     ctTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAattca             >  1:857210/1‑71 (MQ=255)
                                                     ctTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAattca             >  1:834587/1‑71 (MQ=255)
                                                     ctTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAattca             >  1:2016086/1‑71 (MQ=255)
                                                     ctTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAattca             >  1:367010/1‑71 (MQ=255)
                                                     ctTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAattca             >  1:54897/1‑71 (MQ=255)
                                                     ctTTAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAattc              >  1:3357126/1‑70 (MQ=255)
                                                      tttAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        >  1:3148608/1‑43 (MQ=255)
                                                      tttAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        >  1:9502/1‑43 (MQ=255)
                                                      tttAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        >  1:606875/1‑43 (MQ=255)
                                                      tttAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        >  1:2458060/1‑43 (MQ=255)
                                                      tttAGCTAGCCTTGTATCGCATTGAACAAAACTTGAACCGAtt                                        >  1:320381/1‑43 (MQ=255)
                                                        tagctagcCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAattcatt           >  1:2749698/1‑70 (MQ=255)
                                                        tagctagcCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAattcatt           >  1:1467577/1‑70 (MQ=255)
                                                         agctagcCTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAattcattc          <  1:3070376/70‑1 (MQ=255)
                                                          gctagcCTTGCATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAAttcattca         <  1:1306009/70‑1 (MQ=255)
                                                                cTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTgg                           >  1:1869965/1‑46 (MQ=255)
                                                                cTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTgg                           >  1:1747536/1‑46 (MQ=255)
                                                                cTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTgg                           >  1:865037/1‑46 (MQ=255)
                                                                cTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAATTCATTCATTTGACt  <  1:3667635/71‑1 (MQ=255)
                                                                cTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAATTCATTCATTTGACt  <  1:3211833/71‑1 (MQ=255)
                                                                cTTGTATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAATTCATTCATTTGACt  <  1:1287370/71‑1 (MQ=255)
                                                                    |                                                                  
TATCTGGGTTATCAGCGAAAAGTATAAGGGGTAAACAAGGATAAAGTGTCACTCTTTAGCTAGCCTTGCATCGCATTGAACAAAACTTGAACCGATTTAGCAAAACGTGGCATCGGTCAATTCATTCATTTGACT  >  W3110S.gb/1620642‑1620776

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: