Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,504,889 (A)8→7 intergenic (+633/‑92) insB‑25 → / → ECB_04162 b4563(b4576); IS1 protein insB/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,504,8820A.92.7% 287.1 / 18.9 68intergenic (+626/‑99)insB‑25/ECB_04162b4563(b4576); IS1 protein insB/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base . (31/32);  total (33/35)
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

AGATAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCAAAA  >  REL606/4504807‑4504951
                                                                           |                                                                     
aGATAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTcaca                                                                       <  1:219852/76‑1 (MQ=255)
aGATAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTcaca                                                                       >  2:1616687/1‑76 (MQ=255)
 gATAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaa                                                                      >  2:49908/1‑76 (MQ=255)
 gATAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaa                                                                      >  2:791468/1‑76 (MQ=255)
  aTAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaa                                                                     >  1:502737/1‑76 (MQ=255)
  aTAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaa                                                                     <  2:1124594/76‑1 (MQ=255)
    aaTGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaaaa                                                                   >  1:1443857/1‑76 (MQ=255)
         ttttAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaaaaaacca                                                              >  1:335700/1‑73 (MQ=255)
           ttAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCAtt                                                           <  2:390594/76‑1 (MQ=255)
             aGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAg                                                         <  2:1170057/76‑1 (MQ=255)
              gTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGt                                                        <  1:1284287/76‑1 (MQ=255)
               tAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTc                                                       <  1:1432619/76‑1 (MQ=255)
               tAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTc                                                       <  1:1139869/76‑1 (MQ=255)
                aaaaaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCt                                                      <  1:1247942/76‑1 (MQ=255)
                aaaaaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCt                                                      >  1:1076374/1‑76 (MQ=255)
                 aaaaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTg                                                     <  1:36632/76‑1 (MQ=255)
                 aaaaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTg                                                     <  1:36627/76‑1 (MQ=255)
                   aaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCt                                                   <  1:393905/76‑1 (MQ=255)
                   aaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCt                                                   <  1:612494/76‑1 (MQ=255)
                   aaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCt                                                   <  1:612504/76‑1 (MQ=255)
                   aaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCt                                                   >  2:938303/1‑76 (MQ=255)
                    aaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCtt                                                  <  1:1393420/76‑1 (MQ=255)
                       agTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCt                                               >  1:44217/1‑76 (MQ=255)
                       agTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCt                                               <  2:807397/76‑1 (MQ=255)
                       agTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCt                                               <  2:807412/76‑1 (MQ=255)
                        gTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTg                                              <  2:975772/76‑1 (MQ=255)
                        gTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTg                                              <  2:113731/76‑1 (MQ=255)
                        gTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTg                                              <  2:113738/76‑1 (MQ=255)
                          tGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGtt                                            <  2:258711/76‑1 (MQ=255)
                          tGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGtt                                            <  2:1210787/76‑1 (MQ=255)
                              tGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTa                                        <  2:295067/76‑1 (MQ=255)
                              tGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTa                                        <  1:1609211/76‑1 (MQ=255)
                               gAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACCATTAGTCTGCTTCCTGTTCTTa                                        >  1:757664/1‑76 (MQ=255)
                                aaTATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACATTAGTCTGCTTCCTGTTCTTACa                                      <  1:204522/76‑1 (MQ=255)
                                aaTATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACATTAGTCTGCTTCCTGTTCTTACa                                      <  1:204513/76‑1 (MQ=255)
                                aaTATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACa                                      >  2:196599/1‑76 (MQ=255)
                                aaTATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACa                                      >  1:1182591/1‑76 (MQ=255)
                                aaTATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACa                                      >  1:898588/1‑76 (MQ=255)
                                aaTATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACa                                      >  1:898609/1‑76 (MQ=255)
                                  tataATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaat                                    <  2:1396551/76‑1 (MQ=255)
                                    taataCAGCCTCTTCCTACATANGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaataa                                  >  2:231879/1‑76 (MQ=255)
                                    taataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaata                                   >  1:610557/1‑76 (MQ=255)
                                    taataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaataa                                  <  1:1502036/76‑1 (MQ=255)
                                     aataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaataat                                 >  2:670437/1‑76 (MQ=255)
                                     aataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaataat                                 <  1:1283734/76‑1 (MQ=255)
                                       taCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATc                                <  2:677210/76‑1 (MQ=255)
                                       taCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCa                               >  2:72027/1‑76 (MQ=255)
                                         cAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGt                             <  1:951052/76‑1 (MQ=255)
                                         cAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGt                             <  1:327557/76‑1 (MQ=255)
                                            ccTCTTCCTACATAAGAAGGAATGTGCTCGCAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtatt                          <  2:292070/76‑1 (MQ=255)
                                            ccTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtatt                          <  1:1154616/76‑1 (MQ=255)
                                            ccTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtatt                          <  1:327816/76‑1 (MQ=255)
                                            ccTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtatt                          <  1:273436/76‑1 (MQ=255)
                                             ctctTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtatta                         >  1:824254/1‑76 (MQ=255)
                                                   cTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatt                   >  1:1170391/1‑76 (MQ=255)
                                                     aCATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatttt                 >  2:1179869/1‑76 (MQ=255)
                                                     aCATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatttt                 >  2:70426/1‑76 (MQ=255)
                                                      cATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatttta                >  1:134501/1‑76 (MQ=255)
                                                       aTAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttat               >  2:303163/1‑76 (MQ=255)
                                                       aTAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttat               >  2:304686/1‑76 (MQ=255)
                                                        tAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttatt              >  2:1148839/1‑76 (MQ=255)
                                                         aagaagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttattt             >  2:512742/1‑76 (MQ=255)
                                                         aagaagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttattt             >  2:972426/1‑76 (MQ=255)
                                                         aagaagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttattt             <  2:628550/76‑1 (MQ=255)
                                                         aagaagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttattt             >  2:333879/1‑76 (MQ=255)
                                                         aagaagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttattt             >  2:333866/1‑76 (MQ=255)
                                                         aagaagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttattt             <  2:26620/76‑1 (MQ=255)
                                                            aagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtata          >  2:1177258/1‑76 (MQ=255)
                                                             agGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtatat         >  2:33126/1‑76 (MQ=255)
                                                              ggAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtatata        >  1:144208/1‑76 (MQ=255)
                                                               gAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAg       >  1:1581207/1‑76 (MQ=255)
                                                                aaTGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGc      >  2:1127713/1‑76 (MQ=255)
                                                                aaTGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGc      >  2:1127702/1‑76 (MQ=255)
                                                                aaTGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGc      >  2:932907/1‑76 (MQ=255)
                                                                aaTGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGc      >  2:1385801/1‑76 (MQ=255)
                                                                    tgCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCaaaa  >  1:806076/1‑76 (MQ=255)
                                                                           |                                                                     
AGATAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCAAAA  >  REL606/4504807‑4504951

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