Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,847,052 A→G V10A (GTT→GCT)  recD ← exonuclease V (RecBCD complex), alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,847,0520AG98.5% 267.7 / NA 67V10A (GTT→GCT) recDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (25/41);  minor base T (1/0);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAACAGGGCTGCCGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAACTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCATTGGGTC  >  REL606/2846954‑2847147
                                                                                                  |                                                                                               
ttAACAGGGCTGCCGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCt                                                                                               <  2:77559/101‑1 (MQ=255)
   aCAGGGCTGCCGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTcc                                                                                            <  1:436651/101‑1 (MQ=255)
      gggCTGCCGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGt                                                                                         <  1:37977/101‑1 (MQ=255)
         cTGCCGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAAt                                                                                      <  2:181800/101‑1 (MQ=255)
         cTGCCGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAAt                                                                                      <  2:986339/101‑1 (MQ=255)
          tGCCGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAAtt                                                                                     >  2:838436/1‑101 (MQ=255)
            ccGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTCATTGc                                                                                   <  2:1424383/101‑1 (MQ=255)
            ccGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGc                                                                                   <  2:1499344/101‑1 (MQ=255)
               cgAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCttt                                                                                <  1:125436/101‑1 (MQ=255)
                gAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCtttt                                                                               <  1:785483/101‑1 (MQ=255)
                   ggTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCa                                                                            <  1:1310648/101‑1 (MQ=255)
                   ggTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCa                                                                            >  2:765381/1‑101 (MQ=255)
                   ggTGACGGCAGGATGTTCATCCCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCa                                                                            <  2:243131/101‑1 (MQ=255)
                    gTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCaa                                                                           >  2:169160/1‑101 (MQ=255)
                       aCGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAAttt                                                                        <  2:1104471/101‑1 (MQ=255)
                         ggCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCa                                                                      >  2:598583/1‑101 (MQ=255)
                           cAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCAtt                                                                    <  2:1374089/101‑1 (MQ=255)
                            aGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTa                                                                   <  2:1551339/101‑1 (MQ=255)
                            aGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTa                                                                   <  1:185043/101‑1 (MQ=255)
                             ggATGTTCATCTCCCGCCACGGTTAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTAc                                                                  <  1:805110/101‑1 (MQ=255)
                             ggATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTAc                                                                  <  1:866239/101‑1 (MQ=255)
                             ggATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTAc                                                                  <  1:693419/101‑1 (MQ=255)
                             ggATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTAc                                                                  <  2:8569/101‑1 (MQ=255)
                                tGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGGTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCAttcactc                                                               >  2:334314/1‑96 (MQ=255)
                                 tttCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGcct                                                              <  2:40608/100‑1 (MQ=255)
                                 gTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGcct                                                              <  1:405572/101‑1 (MQ=255)
                                 gTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGcct                                                              >  2:879111/1‑101 (MQ=255)
                                 gTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGcct                                                              <  1:785348/101‑1 (MQ=255)
                                   tCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGcctcc                                                            <  2:514349/101‑1 (MQ=255)
                                      tctcCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCtcctcc                                                         <  1:1024614/101‑1 (MQ=255)
                                         cccGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAgg                                                      >  1:432221/1‑101 (MQ=255)
                                             ccACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCa                                                  <  2:999873/101‑1 (MQ=255)
                                             ccACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCa                                                  >  2:629017/1‑101 (MQ=255)
                                              cACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCAt                                                 <  1:838441/101‑1 (MQ=255)
                                              cACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTCCATTACGCCTCCTCCAGGGTCAt                                                 >  1:247862/1‑101 (MQ=255)
                                               aCGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATa                                                <  1:1463760/101‑1 (MQ=255)
                                                cGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATAc                                               <  2:1330029/101‑1 (MQ=255)
                                                 ggTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATAcc                                              <  1:60357/101‑1 (MQ=255)
                                                    cAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGc                                           >  1:885653/1‑101 (MQ=255)
                                                     aGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCa                                          <  2:318288/101‑1 (MQ=255)
                                                        gCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAAc                                       <  1:1271183/101‑1 (MQ=255)
                                                          aaaTTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACAt                                     >  1:789861/1‑101 (MQ=255)
                                                          aaaTTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACAt                                     >  1:876577/1‑101 (MQ=255)
                                                            aTTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACAtct                                   <  1:307316/101‑1 (MQ=255)
                                                             ttGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACAtctc                                  >  2:45758/1‑101 (MQ=255)
                                                              tGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTca                                 <  2:490400/101‑1 (MQ=255)
                                                               gCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTcat                                <  1:1340482/101‑1 (MQ=255)
                                                                 acaTCCAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTcatcc                              <  2:1166002/101‑1 (MQ=255)
                                                                     ccAGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCa                          >  1:639802/1‑101 (MQ=255)
                                                                       aGCGGGCGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAgg                        >  1:1499704/1‑101 (MQ=255)
                                                                         cGGGCGTAGCTGTTTGTGCTCCACATCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGtc                      >  2:811823/1‑101 (MQ=255)
                                                                             cGTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGtcaatc                  <  1:598588/101‑1 (MQ=255)
                                                                              gTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGtcaatca                 >  1:352185/1‑101 (MQ=255)
                                                                              gTAGCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGtcaatca                 >  2:1256572/1‑101 (MQ=255)
                                                                                 gCTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAAcc              >  2:863776/1‑101 (MQ=255)
                                                                                  cTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAAccc             >  1:1317875/1‑101 (MQ=255)
                                                                                  cTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAAccc             >  2:1402302/1‑101 (MQ=255)
                                                                                  cTGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAAccc             <  1:747432/101‑1 (MQ=255)
                                                                                   tGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCt            <  1:832836/101‑1 (MQ=255)
                                                                                   tGTTTGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCt            >  2:1060979/1‑101 (MQ=255)
                                                                                      ttGTGCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCa         <  1:169160/101‑1 (MQ=255)
                                                                                       tgtgCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCAt        >  1:1256329/1‑101 (MQ=255)
                                                                                       cttgCTCCACAGCCTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCAt        <  2:1647062/99‑1 (MQ=255)
                                                                                         tgCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCATTg      <  1:961317/101‑1 (MQ=255)
                                                                                          gCTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCATTgg     <  1:663640/101‑1 (MQ=255)
                                                                                           cTCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCATTggg    >  1:738786/1‑101 (MQ=255)
                                                                                            tCCACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCATTGGGt   >  1:1287712/1‑101 (MQ=255)
                                                                                             ccACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCATTGGGTc  >  1:587676/1‑101 (MQ=255)
                                                                                             ccACAGCTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCATTGGGTc  >  1:371208/1‑101 (MQ=255)
                                                                                                  |                                                                                               
TTAACAGGGCTGCCGCGAGGGTGACGGCAGGATGTTCATCTCCCGCCACGGTCAGGGCAAATTGCACATCCAGCGGGCGTAGCTGTTTGTGCTCCACAACTTCCAGTAATTGCTTTTGCAATTTCATTACGCCTCCTCCAGGGTCATACCGGCAAACATCTCATCCATCAGGTCAATCAACCCTGCATTGGGTC  >  REL606/2846954‑2847147

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: