Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,105,116 (A)8→7 intergenic (+77/+135) ECB_01993 → / ← ECB_01994 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,105,1090A.95.3% 188.8 / 5.3 44intergenic (+70/+142)ECB_01993/ECB_01994hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base . (22/19);  minor base C (0/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

ATATAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACGGGGAGGAC  >  REL606/2105034‑2105180
                                                                            |                                                                       
atatAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCa                                                                         >  2:453628/1‑76 (MQ=255)
atatAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCa                                                                         >  2:1870054/1‑76 (MQ=255)
atatAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCa                                                                         <  2:242685/76‑1 (MQ=255)
atatAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCa                                                                         <  2:1201227/76‑1 (MQ=255)
 tataAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaa                                                                        <  2:282432/76‑1 (MQ=255)
 tataAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaa                                                                        <  2:40800/76‑1 (MQ=255)
  ataAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaa                                                                       <  2:108417/76‑1 (MQ=255)
  aaaaaacCACTAAACTACAGATCCTCCTAACTTGGTACCTGCCCCCCCCCCCCAAAAATTAGGGGTGAGGGGCaaa                                                                       <  2:1485762/74‑1 (MQ=255)
  aaaaaaaCACTTAACTACAGATCCCCCTATTTTTGTTTCGGCCCCCCCCCCCCAAAAATAAGGGGGGGGGGGCaaa                                                                       <  2:1084856/74‑1 (MQ=255)
  aaaaaaaCACTTAACTACAGATACTCCTAACTTTGTATCTGACCCCCCTCACCCTAAATTAAGGGTGAGGGGCaaa                                                                       <  2:446534/74‑1 (MQ=255)
   taAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaaa                                                                      <  2:196653/76‑1 (MQ=255)
   aaaaaaCACTTAACTACAAATACCCCTAACTTTGTATCTGCCGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaaa                                                                      <  2:204602/75‑1 (MQ=255)
   aaaaaaCACTAAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCCCCCAAAATAAAGGGTGAGGGGCaaaa                                                                      <  2:1234072/75‑1 (MQ=255)
    aaaaaCACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaaaa                                                                     <  2:1496705/76‑1 (MQ=255)
      aaaCACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaaaaaa                                                                   <  2:1442527/76‑1 (MQ=255)
       aaCACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaaaaaat                                                                  <  2:12899/76‑2 (MQ=255)
        ccacTTAACAACAAATCCCCCAAACTTTGTATCTGACCCCCCCCCCCCAAAATAAGGGGTGAGGGGCaaaaaaata                                                                 <  2:556519/75‑3 (MQ=255)
        acacTTAACTACAGATCCTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaaaaaata                                                                 <  2:1849730/76‑3 (MQ=255)
         cacTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCCCCCTAAATTAAGGGTGAGGGGCaaaaaaataa                                                                <  2:1666215/76‑4 (MQ=255)
         cacTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCCCACCCCAAATTAAGGGGGAGGGGCaaaaaaaaaa                                                                >  2:475366/1‑74 (MQ=255)
           cTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCCCACCCCCAAATAAAGGGGAGGGGGAAAAAAaaaaat                                                              >  1:1862531/1‑76 (MQ=255)
              aaCTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAattatt                                                          >  2:88928/1‑76 (MQ=255)
               aCTACAGATCCCCCTATTTTTGTATCTGCCCCCCCCCCCCAAAAATAAGGGGGGGGGGGCAAAAAAATAATTAttt                                                         <  2:908612/76‑1 (MQ=255)
                cTACAGATACTCCTAACTTTGTATCTGACGCCCCCCCCCCCTAAATTAAGGGGAAGGGGCAAAAAAAAAATAAttt                                                         >  1:1467878/1‑76 (MQ=255)
                cTACAGATACTCCTAACTTTGTATCTGACGCCCCCCACCCTAAATTAAGGGGGAGGGGCAAAAAAATAATTAtttt                                                        >  2:1570150/1‑75 (MQ=255)
                  aCAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCCAAATTAAGGGGGAGGGGCAAAAAAATAATTATTTGAt                                                      >  2:400032/1‑76 (MQ=255)
                    aGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGt                                                    >  2:649612/1‑76 (MQ=255)
                      aTACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGtta                                                  >  2:251467/1‑76 (MQ=255)
                       tACTCCTAACTTTGTATCTGACGCCCCTCACCCTCAATTAAAGGGGAGGGGGAAAAAAATAATTATTTGATGttat                                                 >  2:855870/1‑76 (MQ=255)
                        aCTCCTAACTTTGTATCTGACGCCCCTCACCCTCAATTTAGGGGGAGGGGCAAAAAAAAAATTATTTGATGttttt                                                >  2:484214/1‑73 (MQ=37)
                          tCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatatt                                              <  2:12512/76‑1 (MQ=255)
                           ccTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatattt                                             <  2:1827683/76‑1 (MQ=255)
                            cTAACTTTGTATCTGACGCCCCCCCCCCCTAAATTAAGGTGAGGGGGCAAAAAAAAAATTTTTTTTGttatttttt                                            >  2:1070660/1‑75 (MQ=255)
                            cTAACTTTGTATCTGACGCCCCCCACCCTAAATTAAAGGGGAGGGGCAAAAAAAAAATTATTTGATGttattttta                                            >  2:596067/1‑76 (MQ=255)
                             tAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatatttat                                           <  2:240964/76‑1 (MQ=255)
                              aaCTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGGGAGGGGCAAAAAAATAATTATTTGATGttatatttata                                          >  2:1389913/1‑76 (MQ=39)
                               aCTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatatttatat                                         <  2:1481367/76‑1 (MQ=255)
                                ttttGTATCTGCCCCCCCTCCCCAAAAATTAGGGGGGGGGGGCAAAAAAATAATTATTTGATGTTATATTTATATg                                        <  2:701784/75‑1 (MQ=255)
                                cTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATg                                        <  2:659355/76‑1 (MQ=255)
                                cTTTGTATCTGACGCCCCCCCCCCCAAATTAAAGGGGAGGGGCAAAAAAAAAATTATTTTATGttatatttatatt                                        >  2:929575/1‑75 (MQ=255)
                                 tttGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAAATAATTATTTGATGTTATATTTATATg                                        <  2:155988/76‑1 (MQ=255)
                                 tttGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGc                                       <  2:279811/76‑1 (MQ=255)
                                 tttGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGc                                       <  2:1845252/76‑1 (MQ=255)
                                 tttGTATCTGACGCCCCTCACCCTAAATTAAGGGGGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGc                                       >  2:219311/1‑76 (MQ=39)
                                 tttGTATCTGACGCCCCTCACCCTAAATTAAAGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGc                                       >  2:1275587/1‑76 (MQ=39)
                                  ttGTTTCGGCCGCCCCTCCCCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGcc                                      <  2:204780/76‑1 (MQ=255)
                                  ttGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGcc                                      <  2:204770/76‑1 (MQ=255)
                                   tGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCa                                     >  2:55943/1‑76 (MQ=255)
                                   tGTATCTGACGCCCCTCACCCTAAATTAAGGGGGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCa                                     >  2:1755513/1‑76 (MQ=255)
                                    gTATCTGACGCCCCCCCCCCCCAATTAAGGGTGAGGGGGAAAAAAAAAAATATTTGATTTTATATTTATATGCCAt                                    >  2:318053/1‑76 (MQ=255)
                                     tATCTGACGCCCCTCACCCTAAATTAAGGGGGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATc                                   >  2:146871/1‑76 (MQ=255)
                                       tCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATcac                                 <  2:240891/76‑1 (MQ=255)
                                          gACGCCCCTCACCCCAAATTAAAGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATCCCATCACATc                              >  2:358396/1‑76 (MQ=255)
                                          gACGCCCCCCACCCTAAATTAAAGGTGAGGGGCAAAAAAATAATTATTTGATGGTATATTTATATTCCATCACATc                              >  2:1127229/1‑76 (MQ=255)
                                            cGCCCCTCACCCAAAATAAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCtt                            <  2:480380/76‑1 (MQ=255)
                                                 ccTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTaa                        <  2:149510/76‑1 (MQ=255)
                                                 ccTCACCCTAAATTAAGGGGGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTaa                        >  2:1072523/1‑76 (MQ=255)
                                                    cACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTc                     >  2:544182/1‑76 (MQ=255)
                                                      cccTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCa                   <  1:1870054/76‑1 (MQ=255)
                                                        caaaaaTAAGGGGGGGGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAg                  <  1:1414145/74‑1 (MQ=255)
                                                        cTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAgt                 >  2:1595904/1‑76 (MQ=255)
                                                         tAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAgtt                >  1:1252871/1‑76 (MQ=255)
                                                         tAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAgtt                <  2:398157/76‑1 (MQ=255)
                                                           aaTTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAgttgt              >  2:1472832/1‑76 (MQ=255)
                                                           aaTTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAgttgt              >  1:1240983/1‑76 (MQ=255)
                                                             ttAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTa            <  2:571907/76‑1 (MQ=255)
                                                             ttAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTa            >  2:98340/1‑76 (MQ=255)
                                                              tAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTAc           <  2:180471/76‑1 (MQ=255)
                                                              tAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTAc           >  1:868606/1‑76 (MQ=255)
                                                              tAAGGGTGAGGGGC‑‑AAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACg          >  1:514835/1‑76 (MQ=255)
                                                               aaGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACg          >  1:892956/1‑76 (MQ=255)
                                                                aGGGTGAGGGGCAAAAAAATAATTATTTGAGGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACgg         <  2:29677/76‑1 (MQ=255)
                                                                 gggTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACggg        >  1:859467/1‑76 (MQ=255)
                                                                       ggggCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACGGGGAGGAc  <  1:1669729/76‑1 (MQ=255)
                                                                            |                                                                       
ATATAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACGGGGAGGAC  >  REL606/2105034‑2105180

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