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.7% 201.1 / 3.4 47intergenic (+70/+142)ECB_01993/ECB_01994hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (19/25);  total (20/27)
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

ATATAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACGGGGAGGAC  >  REL606/2105034‑2105180
                                                                             |                                                                       
atatAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCa                                                                         <  2:1127454/76‑1 (MQ=255)
atatAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCa                                                                         <  2:458876/76‑1 (MQ=255)
 tataAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGCCGCCCCTCACCCTAAATTAAGGGTGAGGGGCaa                                                                        <  2:915404/76‑1 (MQ=255)
 tataAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaa                                                                        <  2:291681/76‑1 (MQ=255)
 tataAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaa                                                                        <  2:706777/76‑1 (MQ=255)
 tataAAAACACTTAACGACAGATACTCCTAACTTGGTATCTGACGCCCCTCACCCTAAATTAGGGGTGAGGGGCaa                                                                        <  2:1052940/76‑1 (MQ=255)
  ataAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaa                                                                       <  2:712520/76‑1 (MQ=255)
  ataAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaa                                                                       <  2:1121289/76‑1 (MQ=255)
  ataAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaa                                                                       >  2:303459/1‑76 (MQ=255)
   aaaaaaCCCTTAACTCAACATCCTCCTAACTTTATTTCGCACCCCCCCCCCCCAAAATTAAGGGTGAGGGGCaaaa                                                                      <  2:674124/75‑1 (MQ=255)
    aaaaaCACTTAATAGCATAGACTTNTAACTTTGTATCTGACGCCCCTCACCCTAAATTNAGGGTGAGGGGCaaaaa                                                                     <  2:957666/76‑1 (MQ=255)
    aaaaaCACTTAACTCCAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaaaa                                                                     >  2:464799/1‑76 (MQ=255)
    aaaaaCACTTAACTACAGATAATCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaaaa                                                                     >  2:333874/1‑76 (MQ=255)
      aaaCACTAAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCaaaaaaa                                                                   <  2:862642/76‑1 (MQ=255)
           cTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAatt                                                             <  2:942006/76‑1 (MQ=255)
            ttAACTACAGATCCTCCTACCTTTGTATCTGCCGCCCCTCACCCAAAATTAAGGGTGAGGGGCAAAAAAATAatta                                                            <  2:850175/76‑1 (MQ=255)
             tAACTACAGATCCTCCTAACTTTGTATCTGACGCCCCTCCCCCAAAATAAAGGGTGAGGGGCAAAAAAATAattat                                                           <  2:412691/76‑1 (MQ=255)
             tAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAattat                                                           <  2:317503/76‑1 (MQ=255)
             tAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGAGGAGGGGCAAAAAAATAattat                                                           >  1:833496/1‑76 (MQ=255)
             tAACTACAGATACCCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAattat                                                           >  2:479979/1‑76 (MQ=255)
              atcTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTNGGGGGCAAAAAAATAattatt                                                          >  2:971267/3‑76 (MQ=255)
              aaCTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACACTAAATTAAGGGTGAGGGGCAAAAAAATAattatt                                                          >  2:271070/1‑76 (MQ=255)
                 tACAGATACTCCTAACTTTGTATCTGACGCCCCCCCCCCCCAAATTAAAGGGAGGGGGGAAAAAAAAAAAatttta                                                          >  1:733868/1‑73 (MQ=255)
                 tACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGa                                                       >  2:879259/1‑76 (MQ=255)
                 tACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGa                                                       >  2:985466/1‑76 (MQ=255)
                 tACAGATACTCCTAACTTTGTATCTGACGCCCCCCCCCCCCCAATAATAGGGAGGGGGGAAAAAAAAAataattaa                                                        >  1:523458/1‑74 (MQ=255)
                  aCAGATACTCCTAACTTTGTATCTGACGCCCCCCACCCCAAATTAAAGGGGAGGGGGAAAAAAAAAAATAAtttat                                                       >  2:570403/1‑74 (MQ=255)
                  aCAGATACTCCTAACTTTATATCTGACGCCCCTCACCCTAAATTAAGCGTGAGGGGCAAAAAAATAATTATTTGAt                                                      >  2:1051520/1‑76 (MQ=38)
                      aTACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGtta                                                  >  2:1081966/1‑76 (MQ=255)
                       tACTCCTAACTTTGTATCTGTCGCCCCTCCCCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttat                                                 <  2:962522/76‑1 (MQ=38)
                        aCTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttata                                                <  2:833496/76‑1 (MQ=255)
                        aCTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttata                                                <  1:333874/76‑1 (MQ=255)
                        aCTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttata                                                <  1:464799/76‑1 (MQ=255)
                         cTCCTAACTTTGTATCTGCCGCCCCTCACCCTAAATTAACGGGGAGGGGATAACAAATAATTATTTGATGttatat                                               >  2:1063675/1‑76 (MQ=255)
                         cTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatat                                               <  2:703648/76‑1 (MQ=255)
                          tCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatatt                                              <  2:78920/76‑1 (MQ=255)
                          tCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatatt                                              <  2:789578/76‑1 (MQ=255)
                          tCCTAACTTTGTATCTGACGCCCATCACCCTAAATTAAGGGTGAGGGGCAAAATAAATAATTATTTGATGttattt                                               <  1:995062/76‑3 (MQ=39)
                           ccTAANTTTGTATCTGACGNCCCTCACCCNAAATTAAGGGTGCGGGGCAAAAAAATAAGTATNTGATNttatattt                                             <  2:190519/76‑1 (MQ=255)
                           ccTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatattt                                             >  2:693864/1‑76 (MQ=255)
                           ccTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatattt                                             >  2:565693/1‑76 (MQ=255)
                           ccTAACTTTGTATCTGACGCCCCTCACCCCAAAATAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatattt                                             >  2:921232/1‑76 (MQ=255)
                             tAACTTTGTATCTGCCCCCCCCCCCCAAAAATAAGGGGTGAGGGGCAAAAAAATAATTATTTGATGttatatttat                                           <  2:582152/76‑1 (MQ=255)
                             tAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGttatatttat                                           >  1:923128/1‑76 (MQ=255)
                              aaCTTTGTATCTGACGCCCCCCACCCTAAATTAAAGGTGAGGGGCAAAAAAATAATTATTTTATGttatatttata                                          >  1:945965/1‑76 (MQ=255)
                                 tttGTATCTGACGCCCTTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGc                                       <  2:82263/76‑1 (MQ=255)
                                 tttGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCNAAAAAATAATTATTTGATGTTATATTTATATGc                                       >  2:1054573/1‑76 (MQ=255)
                                  ttGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGcc                                      >  2:1067174/1‑76 (MQ=255)
                                  ttGTATCTGACGCCCCTCACCCTAAATTAAGGGGGAGGGGCAAAAAAAAAATTATTTGATGTTATATTTATATGcc                                      >  2:107541/1‑76 (MQ=255)
                                   tGTATCTGACGCCCCCCACCCTAAATTAAAGGGGAGGGGCAAAAAAATAATTATTTTATGTTATATTTATATGCCa                                     >  2:787157/1‑76 (MQ=255)
                                    gTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCAt                                    >  2:447325/1‑76 (MQ=255)
                                     tATCTGACGCCCCCCACCCCCAAATAAAGGTGAGGGGGAAAAAAAAAATTATTTGATGTTATATTTATATTTCATc                                   >  2:452414/1‑76 (MQ=255)
                                       tCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATcac                                 <  2:130698/76‑1 (MQ=255)
                                       tCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATcac                                 >  2:106838/1‑76 (MQ=255)
                                          gACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATc                              >  2:502757/1‑76 (MQ=255)
                                          gACGCCCCTCACACTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATc                              >  2:187225/1‑76 (MQ=39)
                                           aCGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATGTGATGTTATATTTATATGCCATCACATCt                             >  1:1044724/1‑76 (MQ=255)
                                            cGCCCCTCCCCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCtt                            <  1:219399/76‑1 (MQ=39)
                                            cGCCCCTCCCCCAAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCtt                            <  2:379028/76‑1 (MQ=255)
                                            cGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCtt                            >  2:275035/1‑76 (MQ=255)
                                               ccccTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCtttt                          <  2:159376/76‑1 (MQ=255)
                                                cccTCACCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTa                         <  2:758321/76‑1 (MQ=255)
                                                     aCCCTAAATTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTcc                    >  2:404053/1‑76 (MQ=255)
                                                          aaaTTAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAgttg               <  2:631086/76‑1 (MQ=255)
                                                          aaaTAAAGGGTGAGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAgttg               <  1:348866/76‑1 (MQ=255)
                                                             ttAAGGGTGAGGGGCAAAAAAATAATTATNTGATGTTGTATTGATNTGCCATCACAGCTTTTAACTCCAGTNGNTa            <  2:975177/76‑1 (MQ=255)
                                                                gaGGGGGGGGGGGCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTAc           <  1:570403/75‑1 (MQ=255)
                                                                   ggTGAGGGGCAAAAAAAGAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACgggg       <  2:386660/76‑1 (MQ=255)
                                                                      ggggggCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACGGggagg    <  1:452414/74‑1 (MQ=255)
                                                                      gAGGGGCAAAAAAATAATTATTTGATGTTATATCTATATGCCATCACATCTTTTAACTCCAGTTGTTGCGGggagg    <  1:914407/76‑1 (MQ=255)
                                                                        ggggCAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACGGGGAGGAc  <  2:805639/76‑1 (MQ=255)
                                                                             |                                                                       
ATATAAAAACACTTAACTACAGATACTCCTAACTTTGTATCTGACGCCCCTCACCCTAAATTAAGGGTGAGGGGCAAAAAAAATAATTATTTGATGTTATATTTATATGCCATCACATCTTTTAACTCCAGTTGTTACGGGGAGGAC  >  REL606/2105034‑2105180

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