Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,085 T→C 50.9% S233S (AGT→AGC araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,0850TC50.9% 0.1 / 74.7 53S233S (AGT→AGCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (10/16);  new base C (13/14);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTT  >  W3110S.gb/71017‑71151
                                                                    |                                                                  
cTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGt                                                                    >  1:1563712/1‑69 (MQ=255)
 tGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAg                                                                 >  1:210406/1‑71 (MQ=255)
  gTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAgg                                                                >  1:1510592/1‑71 (MQ=255)
   tCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAg                                                                 <  1:585516/69‑1 (MQ=255)
     acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCg                                                              <  1:1222541/70‑1 (MQ=255)
      caTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCTAGGACCAACGTATCAGCCAGGCGaa                                                            <  1:1519099/71‑1 (MQ=255)
      caTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                            <  1:1420640/71‑1 (MQ=255)
      caTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGaa                                                            <  1:1291470/71‑1 (MQ=255)
         cTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAgct                                                         >  1:440793/1‑71 (MQ=255)
          ttttCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAgctg                                                        >  1:1601735/1‑71 (MQ=255)
          ttttCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgct                                                         >  1:224049/1‑70 (MQ=255)
              ccGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCttt                                                    >  1:74969/1‑71 (MQ=255)
               cGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtttt                                                   <  1:588608/71‑1 (MQ=255)
                gCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTg                                                  <  1:561416/71‑1 (MQ=255)
                 ccAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGa                                                 <  1:1086873/71‑1 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                <  1:1142815/71‑1 (MQ=255)
                    gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCa                                              >  1:1461745/1‑71 (MQ=255)
                    gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAgga                                              >  1:1293048/1‑69 (MQ=255)
                     cagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAc                                             <  1:190585/71‑1 (MQ=255)
                       gTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAcc                                            >  1:1466933/1‑69 (MQ=255)
                             gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCtt                                                     <  1:1543653/55‑1 (MQ=255)
                              aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgc                                                          <  1:1654700/49‑1 (MQ=255)
                               ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGa                                                             <  1:1566040/45‑1 (MQ=255)
                               ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgct                                                         <  1:1771424/49‑1 (MQ=255)
                                 aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTg                                                  >  1:1484843/1‑54 (MQ=255)
                                 aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTg                                                  >  1:202203/1‑54 (MQ=255)
                                   cGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTa                               >  1:956233/1‑71 (MQ=255)
                                     tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATc                             <  1:1094685/71‑1 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgctgc                                                       <  1:112259/43‑1 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                       <  1:763568/59‑1 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                       <  1:761552/59‑1 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                       <  1:141067/59‑1 (MQ=255)
                                        tAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGc                           >  1:1198161/1‑70 (MQ=255)
                                         aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGTTGCTTTTGAGCACCACCCGGATGCCTATCGc                           <  1:1593218/69‑1 (MQ=255)
                                         aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                                >  1:846008/1‑48 (MQ=255)
                                         aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                                >  1:1725026/1‑48 (MQ=255)
                                         aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCg                            <  1:156642/68‑1 (MQ=255)
                                         aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGcc                          <  1:830762/70‑1 (MQ=255)
                                         aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                         >  1:694422/1‑71 (MQ=255)
                                          aGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCAc                        <  1:1459624/71‑1 (MQ=255)
                                             tGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgt                     >  1:656082/1‑71 (MQ=255)
                                               gcgcGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcg                   <  1:94546/71‑1 (MQ=255)
                                               gcgcGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcg                   <  1:1483129/71‑1 (MQ=255)
                                                 gcgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCgg                                      <  1:1268265/50‑1 (MQ=255)
                                                 gcgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCgtcggt                 >  1:426200/1‑71 (MQ=255)
                                                   gAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAc                                             >  1:1174847/1‑41 (MQ=37)
                                                   gAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCgg                                      <  1:1547279/48‑1 (MQ=255)
                                                    aGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                       <  1:39652/46‑1 (MQ=255)
                                                    aGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGc              <  1:1492563/71‑1 (MQ=255)
                                                      gACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCa                                              >  1:1312121/1‑37 (MQ=25)
                                                         cAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                       <  1:151285/41‑1 (MQ=25)
                                                         cAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGt         >  1:1668924/1‑71 (MQ=255)
                                                            cGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAGCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTgg      >  1:1356362/1‑71 (MQ=255)
                                                                ttAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt  >  1:1244100/1‑71 (MQ=255)
                                                                    |                                                                  
CTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTT  >  W3110S.gb/71017‑71151

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: