Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,602,152 2 bp→AT intergenic (+175/+52) pinQ → / ← ECB_01505 predicted site‑specific recombinase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,602,1520TA98.8% 331.6 / ‑3.5 83intergenic (+175/+53)pinQ/ECB_01505predicted site‑specific recombinase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (35/47);  total (35/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
*REL6061,602,1530GT98.8% 327.0 / ‑3.7 83intergenic (+176/+52)pinQ/ECB_01505predicted site‑specific recombinase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (35/47);  total (35/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TGCAGTTGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCTGGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCAA  >  REL606/1602105‑1602222
                                               ||                                                                     
tgCAGTTGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTtcctcc                                            <  2:6196064/76‑1 (MQ=255)
tgCAGTTGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTtcctcc                                            <  2:770563/76‑1 (MQ=255)
tgCAGTTGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTtcctcc                                            <  2:770564/76‑1 (MQ=255)
  cAGTTGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGt                                          >  1:1063823/1‑76 (MQ=255)
   aGTTGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGtt                                         <  2:2315904/76‑1 (MQ=255)
   aGTTGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGtt                                         >  1:4763606/1‑76 (MQ=255)
     ttGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGc                                       <  2:6300231/76‑1 (MQ=255)
      tGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCa                                      >  2:3411361/1‑76 (MQ=255)
      tGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCa                                      <  1:282968/76‑1 (MQ=255)
        aCGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCata                                    <  2:5817806/76‑1 (MQ=255)
         cGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCatat                                   <  2:2375910/76‑1 (MQ=255)
         cGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCatat                                   <  1:2240931/76‑1 (MQ=255)
         cGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCatat                                   >  1:2522390/1‑76 (MQ=255)
         cGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCatat                                   <  1:5892438/76‑1 (MQ=255)
           tATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATAc                                 <  1:1601603/76‑1 (MQ=255)
           tATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATAc                                 >  2:1838676/1‑76 (MQ=255)
           tATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATAc                                 <  1:2728526/76‑1 (MQ=255)
           tATCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATAc                                 >  2:255326/1‑76 (MQ=255)
            aTCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACa                                <  1:518610/76‑1 (MQ=255)
            aTCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACa                                >  2:4074005/1‑76 (MQ=255)
            aTCAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACa                                <  1:4298858/76‑1 (MQ=255)
              cAACACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGc                              >  1:2390782/1‑76 (MQ=255)
               aaCACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCt                             <  1:2599392/76‑1 (MQ=255)
               aaCACGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCt                             >  1:1567369/1‑76 (MQ=255)
                acacGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTc                            <  1:3370861/76‑1 (MQ=255)
                acacGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTc                            <  2:3882051/76‑1 (MQ=255)
                 cacGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCa                           <  2:4404085/76‑1 (MQ=255)
                 cacGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCa                           >  2:876187/1‑76 (MQ=255)
                 cacGGTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCa                           <  2:4332160/76‑1 (MQ=255)
                     gTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCtt                       >  2:3326213/1‑76 (MQ=255)
                     gTTCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCtt                       <  2:3484692/76‑1 (MQ=255)
                      ttCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCttt                      >  1:1732151/1‑76 (MQ=255)
                      ttCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCttt                      <  2:4417051/76‑1 (MQ=255)
                       tCAGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCtttt                     <  2:5668884/76‑1 (MQ=255)
                         aGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCtttttt                   <  1:5345227/76‑1 (MQ=255)
                         aGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCtttttt                   >  2:5723004/1‑76 (MQ=255)
                         aGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCtttttt                   <  2:6044238/76‑1 (MQ=255)
                         aGCAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCtttttt                   >  1:3327984/1‑76 (MQ=255)
                          gCAACACCCGATACTTTTTCCTGGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTc                  <  1:2231130/76‑1 (MQ=255)
                           cAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCa                 <  2:3888384/76‑1 (MQ=255)
                           cAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCa                 >  2:4189275/1‑76 (MQ=255)
                           cAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCa                 <  1:5629392/76‑1 (MQ=255)
                           cAACACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCa                 >  2:1181244/1‑76 (MQ=255)
                            aaCACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCAc                >  2:622940/1‑76 (MQ=255)
                            aaCACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCAc                >  1:1139251/1‑76 (MQ=255)
                            aaCACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCAc                <  1:5360790/76‑1 (MQ=255)
                            aaCACCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCAc                >  1:5239519/1‑76 (MQ=255)
                             acacCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCAcc               >  2:3835753/1‑76 (MQ=255)
                              cacCCGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCAccc              <  1:4374304/76‑1 (MQ=255)
                              cacCCGATACTTTTTCCATGCTTAAAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCAccc              >  2:4194965/1‑76 (MQ=11)
                              cacCCGATACTTGTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCAccc              <  2:667387/76‑1 (MQ=11)
                               ccccGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCa             <  2:984581/75‑1 (MQ=255)
                                cccGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCAc            <  1:5106451/76‑1 (MQ=255)
                                 ccGATACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACg           <  1:3375390/76‑1 (MQ=255)
                                    aTACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGAtt        <  2:532670/76‑1 (MQ=255)
                                    aTACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGAtt        >  2:2960614/1‑76 (MQ=255)
                                    aTACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGAtt        >  2:4356431/1‑76 (MQ=255)
                                     tACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTa       <  1:5081699/76‑1 (MQ=255)
                                     tACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTa       <  2:3809602/76‑1 (MQ=255)
                                     tACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTa       >  2:2701736/1‑76 (MQ=255)
                                     tACTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTa       <  1:6214923/76‑1 (MQ=255)
                                      aCTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTaa      <  2:1893008/76‑1 (MQ=255)
                                      aCTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTaa      <  1:241400/76‑1 (MQ=255)
                                      aCTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTaa      >  1:1545774/1‑76 (MQ=255)
                                      aCTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTaa      <  1:729286/76‑1 (MQ=255)
                                      aCTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTaa      >  2:924768/1‑76 (MQ=255)
                                      aCTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTaa      <  2:1992867/76‑1 (MQ=255)
                                       cTTTTTCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAAc     >  2:249290/1‑76 (MQ=255)
                                        tttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAAcc    >  2:2981666/1‑76 (MQ=255)
                                        tttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAAcc    <  2:451087/76‑1 (MQ=255)
                                        tttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAAcc    <  1:5267713/76‑1 (MQ=255)
                                        tttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAAcc    <  1:3475471/76‑1 (MQ=255)
                                        tttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAAcc    >  1:1846850/1‑76 (MQ=255)
                                        tttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAAcc    >  1:1717509/1‑76 (MQ=255)
                                        tttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAAcc    >  1:1606863/1‑76 (MQ=255)
                                         ttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCa   >  1:5806226/1‑76 (MQ=255)
                                         ttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCa   <  2:4086093/76‑1 (MQ=255)
                                         ttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCa   <  1:5109913/76‑1 (MQ=255)
                                         ttttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCa   >  1:2581430/1‑76 (MQ=255)
                                          tttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCaa  >  1:3955730/1‑76 (MQ=255)
                                          tttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCaa  >  1:3356258/1‑76 (MQ=255)
                                          tttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCaa  <  1:4003975/76‑1 (MQ=255)
                                          tttCCATGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCaa  <  2:3878946/76‑1 (MQ=255)
                                               ||                                                                     
TGCAGTTGACGTATCAACACGGTTCAGCAACACCCGATACTTTTTCCTGGCTTCCAGCAACGGGGTTTCTTCCTCCGTTGCATATACAGCTCACCTTTTTTCACCCACGATTAACCAA  >  REL606/1602105‑1602222

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