Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,400,364 T→C D21G (GAT→GGT)  abgB ← predicted peptidase, aminobenzoyl‑glutamate utilization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,400,3640TC98.0% 197.6 / ‑3.9 50D21G (GAT→GGT) abgBpredicted peptidase, aminobenzoyl‑glutamate utilization protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/26);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTC  >  REL606/1400291‑1400435
                                                                         |                                                                       
cGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCg                                                                       <  1:835655/76‑1 (MQ=255)
   aGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGta                                                                    >  2:431304/1‑76 (MQ=255)
      cAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATaa                                                                 >  2:619310/1‑76 (MQ=255)
      cAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATaa                                                                 <  1:638728/76‑1 (MQ=255)
       ataTGCTCCGCTGACCACAACTCATGAAAACGTGTTTCTNGATGATCCCAGATTTNATCGGCAATACCGGTATAAc                                                                >  2:1693519/3‑76 (MQ=255)
           gCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTg                                                            >  2:1200627/1‑76 (MQ=255)
              ccGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCg                                                         <  2:912239/76‑1 (MQ=255)
              ccGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCg                                                         <  2:818646/76‑1 (MQ=255)
                  tGACCAGAACTCTTCGAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCg                                                     <  1:1415187/76‑1 (MQ=255)
                   gACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCgg                                                    >  2:1280939/1‑76 (MQ=255)
                   gACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCgg                                                    >  2:733394/1‑76 (MQ=255)
                   gACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCgg                                                    >  2:383268/1‑76 (MQ=255)
                       aGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTc                                                <  1:447651/76‑1 (MQ=255)
                       aGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTc                                                <  2:617317/76‑1 (MQ=255)
                       aGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTc                                                <  2:162293/76‑1 (MQ=255)
                        gAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCa                                               <  1:617807/76‑1 (MQ=255)
                           ctctTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATc                                            >  1:196758/1‑76 (MQ=255)
                            tctTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCg                                           >  1:380517/1‑76 (MQ=255)
                             ctTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGc                                          >  2:638284/1‑76 (MQ=255)
                              ttCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCa                                         <  2:933851/76‑1 (MQ=255)
                               tCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCAt                                        >  2:666215/1‑76 (MQ=255)
                                  aaaCGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgt                                     <  1:422511/76‑1 (MQ=255)
                                  aaaCGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgt                                     >  1:286747/1‑76 (MQ=255)
                                  aaaCGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgt                                     <  1:1297006/76‑1 (MQ=255)
                                    aCGTGTTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgtcg                                   <  2:646328/76‑1 (MQ=255)
                                        gtTTCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaa                               <  1:296464/76‑1 (MQ=255)
                                          ttCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaaaa                             <  1:765210/76‑1 (MQ=255)
                                          ttCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaaaa                             >  2:1143450/1‑76 (MQ=255)
                                          ttCTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaaaa                             >  1:1151050/1‑76 (MQ=255)
                                            ntGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACg                           >  1:473979/2‑76 (MQ=255)
                                            cTGGATGATCCCAGATTTGATCGGGAATACCGGTATAACGTTGGCGATCGGCTGTACTGGCGTCGTCGATaaaaag                           >  1:76753/1‑74 (MQ=255)
                                            cTGGATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACg                           >  1:1354511/1‑76 (MQ=255)
                                              ggATGATCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGAt                         <  1:579535/76‑1 (MQ=255)
                                               gatgatCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATa                        <  1:1362448/76‑1 (MQ=255)
                                                atgatCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAg                       >  2:747042/1‑76 (MQ=255)
                                                atgatCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAg                       <  1:1353303/76‑1 (MQ=255)
                                                 tgatCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGa                      <  1:496756/76‑1 (MQ=255)
                                                  gatCCCAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGAt                     >  1:1188898/1‑76 (MQ=255)
                                                     cccAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTc                  <  1:1143452/76‑1 (MQ=255)
                                                       cAGATTTGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCt                <  1:1200629/76‑1 (MQ=255)
                                                          aTTTGATCGGCAATACCGGTATAACGTTGGCGATCAGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCa             <  2:1370849/76‑1 (MQ=255)
                                                           tttGATCGGCAATACCGGCATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCAt            >  1:1003658/1‑76 (MQ=255)
                                                            ttGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATc           <  2:1421959/76‑1 (MQ=255)
                                                             tGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCNTGCATCa          <  2:1375690/76‑1 (MQ=255)
                                                             tGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCa          <  1:1423155/76‑1 (MQ=255)
                                                             tGATCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCa          <  1:765274/76‑1 (MQ=255)
                                                                tCGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGAt       >  1:1205495/1‑76 (MQ=255)
                                                                 cGGCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATa      >  2:787513/1‑76 (MQ=255)
                                                                   gCAATACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATAcc    >  2:1314661/1‑76 (MQ=255)
                                                                     aaTACCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTc  >  2:161582/1‑76 (MQ=255)
                                                                         |                                                                       
CGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTC  >  REL606/1400291‑1400435

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