Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 453,565 A→T 21.9% intergenic (+112/+38) phr → / ← ybgH deoxyribodipyrimidine photolyase, FAD‑binding/predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE453,5650AT21.9% 127.3 / 33.7 64intergenic (+112/+38)phr/ybgHdeoxyribodipyrimidine photolyase, FAD‑binding/predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base T (5/9);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGA  >  minE/453499‑453634
                                                                  |                                                                     
cGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGAcgcg                                                                   >  1:569786/1‑71 (MQ=255)
   cTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCa                                                                <  1:5168356/71‑1 (MQ=255)
   cTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCa                                                                <  1:1470684/71‑1 (MQ=255)
    ttATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTc                                                                 >  1:3248053/1‑69 (MQ=255)
      aTCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGc                                                             >  1:1796233/1‑71 (MQ=255)
        ccAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCaa                                                               >  1:4139378/1‑67 (MQ=255)
              tACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGc                                                             <  1:2398575/63‑1 (MQ=255)
              tACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCAt                                                      >  1:3518111/1‑70 (MQ=255)
                cAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATc                                                     <  1:3902978/69‑1 (MQ=255)
                     ttGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGc                                                             >  1:3092361/1‑56 (MQ=255)
                          cAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTg                                          <  1:1882276/70‑1 (MQ=255)
                          cAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTg                                          <  1:469339/70‑1 (MQ=255)
                          cAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGc                                         >  1:3704886/1‑71 (MQ=255)
                              ttCAATAAATTGCAGGAAATACGTAGGCCTGATAAGa                                                                       >  1:2792399/1‑37 (MQ=255)
                              ttCAATAAATTGCAGGAAATACGTAGGCCTGATAAGa                                                                       >  1:1779094/1‑37 (MQ=255)
                              ttCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCgg                                            >  1:3490295/1‑64 (MQ=255)
                               tCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCa                                                    >  1:2564631/1‑55 (MQ=255)
                               tCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCa                                                    >  1:108133/1‑55 (MQ=255)
                                cAATAAATTGCAGGAAATACGTAGGCCTGATAAGa                                                                       >  1:2901243/1‑35 (MQ=255)
                                cAATAAATTGCAGGAAATACGTAGGCCTGATAAGa                                                                       >  1:2407228/1‑35 (MQ=255)
                                   tAAATTGCAGGAAATACGTAGGCCTGATAAGCCGCGTCAAg                                                              >  1:2380836/1‑41 (MQ=255)
                                   tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAg                                                              >  1:2052126/1‑41 (MQ=255)
                                   tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAg                                                              >  1:2242683/1‑41 (MQ=255)
                                   tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAg                                                              >  1:5636674/1‑41 (MQ=255)
                                   tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGc                                                        <  1:1023345/47‑1 (MQ=255)
                                   tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGc                                                        <  1:1872565/47‑1 (MQ=255)
                                   tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCAt                                                      <  1:5446975/49‑1 (MQ=255)
                                   tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTAtt                                <  1:1772202/71‑1 (MQ=255)
                                    aaaTTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCg                                             >  1:5413041/1‑57 (MQ=255)
                                     aaTTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGc                                                        >  1:1411317/1‑45 (MQ=255)
                                      aTTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGc                                                        <  1:1704320/44‑1 (MQ=255)
                                      aTTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTaa                              <  1:1819958/70‑1 (MQ=255)
                                       ttGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTaaa                             >  1:4168872/1‑69 (MQ=255)
                                       ttGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTaa                              >  1:3083005/1‑69 (MQ=255)
                                        tGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATc                                                     >  1:3792552/1‑45 (MQ=255)
                                        tGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATc                                                     >  1:1146138/1‑45 (MQ=255)
                                          cAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTc                         <  1:4068127/71‑1 (MQ=255)
                                           aGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTcc                        <  1:3233539/71‑1 (MQ=255)
                                              aaaTACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagc                     >  1:4814873/1‑71 (MQ=255)
                                                ccccgtaGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtctt                  <  1:5162831‑M1/68‑43 (MQ=255)
                                                ccccgtaGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtct                   >  1:5376657‑M1/4‑29 (MQ=255)
                                                ccccgtaGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtct                   >  1:354102‑M1/4‑29 (MQ=255)
                                                aTACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgc                   >  1:2764064/1‑71 (MQ=255)
                                                aTACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgc                   <  1:1415049/71‑1 (MQ=255)
                                                 cccgtAGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtctta                 <  1:701243‑M1/69‑44 (MQ=255)
                                                   cGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgccag                <  1:4830765/71‑1 (MQ=255)
                                                   cGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgcca                 <  1:297224/70‑1 (MQ=255)
                                                    gTAGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtcttatcc              <  1:3553236‑M1/71‑47 (MQ=255)
                                                    gTAGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtcttatcc              <  1:3546526‑M1/71‑47 (MQ=255)
                                                    gTAGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtcttatcc              <  1:845347‑M1/71‑47 (MQ=255)
                                                    gTAGGCCTGATAAGTCACGTCAGCGgcgcatcaggcatcggttgtc                                       >  1:887445‑M1/1‑25 (MQ=255)
                                                     tAGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtcttatcc              >  1:4465251‑M1/1‑24 (MQ=255)
                                                     tAGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtcttatcc              >  1:4976579‑M1/1‑24 (MQ=255)
                                                        gCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgc                            <  1:987004‑M1/53‑33 (MQ=255)
                                                        gCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgc                            <  1:553258‑M1/53‑33 (MQ=255)
                                                           tGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgc                            <  1:3753859‑M1/50‑33 (MQ=255)
                                                           tGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgc                            <  1:906803‑M1/50‑33 (MQ=255)
                                                           tGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTaa                              <  1:2366748/49‑1 (MQ=255)
                                                           tGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTaa                              <  1:3030704/49‑1 (MQ=255)
                                                            gATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTcca                       <  1:1787791/55‑1 (MQ=255)
                                                            gATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTcca                       <  1:2447993/55‑1 (MQ=255)
                                                             aTAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGc                                         <  1:3682591/36‑1 (MQ=255)
                                                             aTAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGtt       >  1:1308465/1‑70 (MQ=255)
                                                             aTAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGtt       <  1:3756168/70‑1 (MQ=255)
                                                             aTAAGACGCGTAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGtt       <  1:835991/69‑1 (MQ=255)
                                                               aaGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccag                      >  1:5120385/1‑53 (MQ=255)
                                                                aGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTg   <  1:941609/71‑1 (MQ=255)
                                                                 gACGCGTCAAGCGTCGCATCAGGCATAGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGa  <  1:5637452/71‑1 (MQ=255)
                                                                  |                                                                     
CGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGA  >  minE/453499‑453634

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: