Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,202,158 A→C I134M (ATT→ATG iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,202,1580AC97.6% 168.1 / NA 42I134M (ATT→ATGiclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (18/23);  minor base T (0/1);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAAATT  >  REL606/4202085‑4202229
                                                                         |                                                                       
aaGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGatcat                                                                       >  1:513047/1‑76 (MQ=255)
aaGCGGGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGatcat                                                                       <  1:195020/76‑1 (MQ=255)
   ggggCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATcataat                                                                    <  1:59335/73‑1 (MQ=255)
   cGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATcataat                                                                    >  2:1043790/1‑76 (MQ=255)
     tGCATCGGCAATTTCCCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCg                                                                  <  1:87489/76‑1 (MQ=255)
       cATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCt                                                                >  1:566822/1‑76 (MQ=255)
       cATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCt                                                                >  2:446110/1‑76 (MQ=255)
          cGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCg                                                             <  2:1727763/76‑1 (MQ=255)
               ttttACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATc                                                        <  1:97834/75‑1 (MQ=255)
                tttACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCg                                                       >  2:1107766/1‑76 (MQ=255)
                tttACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCg                                                       >  2:13935/1‑76 (MQ=255)
                    ccgccgATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCttt                                                   <  1:1578507/76‑1 (MQ=255)
                     cgccgATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTg                                                  <  2:307607/76‑1 (MQ=255)
                     cgccgATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTg                                                  <  1:137205/76‑1 (MQ=255)
                       ccgATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGAt                                                <  2:1338436/76‑1 (MQ=255)
                       ccgATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGAt                                                <  2:1322749/76‑1 (MQ=255)
                        cgATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATc                                               >  2:371234/1‑76 (MQ=255)
                         gATAGGGGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCa                                              <  1:1424889/76‑1 (MQ=255)
                         gATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCa                                              <  1:1049518/76‑1 (MQ=255)
                         gATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCa                                              <  2:996305/76‑1 (MQ=255)
                          aTAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCaa                                             <  1:90492/76‑1 (MQ=255)
                             ggCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCa                                          <  2:442047/76‑1 (MQ=255)
                             ggCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCa                                          >  1:989856/1‑76 (MQ=255)
                              gcgcGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCAc                                         <  2:1435255/76‑1 (MQ=255)
                               cgcgGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCAcc                                        <  2:806808/76‑1 (MQ=255)
                                 cgGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGc                                      <  2:1345783/76‑1 (MQ=255)
                                  ggACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGcc                                     <  2:595540/76‑1 (MQ=255)
                                   gACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCa                                    >  2:454145/1‑76 (MQ=255)
                                   gACATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCa                                    <  1:758001/76‑1 (MQ=255)
                                    aCATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCat                                   <  1:426552/76‑1 (MQ=255)
                                     cATTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCata                                  >  1:1352591/1‑76 (MQ=255)
                                      aTTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCatat                                 <  1:1692992/76‑1 (MQ=255)
                                      aTTCGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCatat                                 >  1:655420/1‑76 (MQ=255)
                                         cGCATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGa                              <  2:133057/76‑1 (MQ=255)
                                           cATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGAcc                            >  1:1555680/1‑76 (MQ=255)
                                           cATCAGATGCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGAcc                            >  2:1342474/1‑76 (MQ=255)
                                                   gCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGcc                    <  1:989368/76‑1 (MQ=255)
                                                   gCGTACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGcc                    <  1:687698/76‑1 (MQ=255)
                                                     gTACGCTGTCCCGGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAg                  <  1:776471/76‑1 (MQ=255)
                                                     gTACACTGTCCCGGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAg                  <  1:731149/76‑1 (MQ=255)
                                                      tACACTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGa                 >  1:1307585/1‑76 (MQ=255)
                                                       acacTGTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGAc                >  1:432548/1‑76 (MQ=255)
                                                        cacTTTCCCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGAct               <  2:1748000/76‑1 (MQ=255)
                                                        cacTTTACTTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGAct               <  1:1160722/76‑1 (MQ=255)
                                                          cTGTCCCTGGTCGATCATATTCGCTTCGTGATCGTTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGActct             <  1:1700267/76‑1 (MQ=255)
                                                            gTACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTc           >  1:589403/1‑76 (MQ=255)
                                                             tACCTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTcc          >  2:1229478/1‑76 (MQ=255)
                                                              cccTGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCCCCAGACTCTTCCa         <  1:1606553/75‑1 (MQ=255)
                                                               ccTGGTCGATTATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCAt        <  2:510029/76‑1 (MQ=255)
                                                                 tGGTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTa      <  2:432548/76‑1 (MQ=255)
                                                                 tGGTCGATCATAACCGCTTCGTAATCGCTTTGATAAAGCCCCGCCATATTGACCGTTTCGCCAGACTCTTCCATTa      <  1:416933/76‑1 (MQ=255)
                                                                  ggTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTaa     >  2:1123314/1‑76 (MQ=255)
                                                                  ggTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTaa     <  2:1273582/76‑1 (MQ=255)
                                                                   gTCGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTaaa    >  1:925537/1‑76 (MQ=255)
                                                                     cGATCATAATCGCTTCGTTATCCCTTTGAACAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAAAtt  <  2:1765064/76‑1 (MQ=255)
                                                                     cGATCATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAAAtt  >  2:1410642/1‑76 (MQ=255)
                                                                         |                                                                       
AAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAAATT  >  REL606/4202085‑4202229

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