Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 453,572 Δ1 bp 22.2% intergenic (+119/+31) phr → / ← ybgH deoxyribodipyrimidine photolyase, FAD‑binding/predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported453,5720A.22.2% 142.5 / 48.3 63intergenic (+119/+31)phr/ybgHdeoxyribodipyrimidine photolyase, FAD‑binding/predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (21/28);  new base . (5/9);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCA  >  Exported/453502‑453640
                                                                      |                                                                    
cTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCa                                                                      <  1:5168356/71‑1 (MQ=255)
cTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCa                                                                      <  1:1470684/71‑1 (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:469339/70‑1 (MQ=255)
                       cAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTg                                                <  1:1882276/70‑1 (MQ=255)
                       cAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGc                                               >  1:3704886/1‑71 (MQ=255)
                           ttCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCgg                                                  >  1:3490295/1‑64 (MQ=255)
                            tCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCa                                                          >  1:2564631/1‑55 (MQ=255)
                            tCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCa                                                          >  1:108133/1‑55 (MQ=255)
                                tAAATTGCAGGAAATACGTAGGCCTGATAAGCCGCGTCAAg                                                                    >  1:2380836/1‑41 (MQ=255)
                                tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAg                                                                    >  1:5636674/1‑41 (MQ=255)
                                tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAg                                                                    >  1:2242683/1‑41 (MQ=255)
                                tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAg                                                                    >  1:2052126/1‑41 (MQ=255)
                                tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGc                                                              <  1:1872565/47‑1 (MQ=255)
                                tAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGc                                                              <  1:1023345/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:1146138/1‑45 (MQ=255)
                                     tGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATc                                                           >  1:3792552/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:1415049/71‑1 (MQ=255)
                                             aTACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgc                         >  1:2764064/1‑71 (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:845347‑M1/71‑47 (MQ=255)
                                                 gTAGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtcttatcc                    <  1:3546526‑M1/71‑47 (MQ=255)
                                                 gTAGGCCTGATAAGTCGCGTCAGCGgcgcatcaggcatcggttgtcggatgcgacgctatgcgtcttatcc                    <  1:3553236‑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:3030704/49‑1 (MQ=255)
                                                        tGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTaa                                    <  1:2366748/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:3756168/70‑1 (MQ=255)
                                                          aTAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGtt             >  1:1308465/1‑70 (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)
                                                                 gcgTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAAtt     <  1:1500171/71‑1 (MQ=255)
                                                                   gTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAgcgcg                  <  1:4447817/56‑1 (MQ=255)
                                                                    tCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTg         <  1:1160827/64‑1 (MQ=255)
                                                                    tCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCa  <  1:1189997/71‑1 (MQ=255)
                                                                      |                                                                    
CTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCA  >  Exported/453502‑453640

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: