Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 670,789 A→G A452A (GCA→GCG ECB_00624 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606670,7890AG98.3% 706.9 / ‑1.4 177A452A (GCA→GCGECB_00624conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (48/125);  total (50/127)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 179 total aligned reads displayed.

ATGTGGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCACTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAATACCATC  >  REL606/670716‑670861
                                                                         |                                                                        
aTGTGGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCt                                                                        <  1:5482913/76‑1 (MQ=255)
aTGTGGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCt                                                                        <  1:4014906/76‑1 (MQ=255)
aTGTGGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCt                                                                        <  2:5127481/76‑1 (MQ=255)
aTGTGGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCt                                                                        <  1:4207677/76‑1 (MQ=255)
aTGTGGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCt                                                                        <  1:3051016/76‑1 (MQ=255)
aTGTGGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCt                                                                        <  1:2588079/76‑1 (MQ=255)
 tgtgGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCtt                                                                       <  2:3471004/76‑1 (MQ=255)
 tgtgGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCtt                                                                       <  2:2884358/76‑1 (MQ=255)
  gtgGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCttt                                                                      <  1:1470981/76‑1 (MQ=255)
    ggCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTcc                                                                    >  1:3358545/1‑76 (MQ=255)
    ggCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTcc                                                                    <  2:3862692/76‑1 (MQ=255)
    ggCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTcc                                                                    <  1:3153536/76‑1 (MQ=255)
     gCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTccc                                                                   <  1:4235684/76‑1 (MQ=255)
       ggATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCg                                                                 <  1:3096396/76‑1 (MQ=255)
        gATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCgg                                                                <  2:4296564/76‑1 (MQ=255)
          tatCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGAt                                                              <  2:605589/76‑1 (MQ=255)
           atCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATg                                                             <  2:1215639/76‑1 (MQ=255)
           atCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATg                                                             <  2:4162412/76‑1 (MQ=255)
           atCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATg                                                             <  1:4966279/76‑1 (MQ=255)
           atCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATg                                                             <  1:5522557/76‑1 (MQ=255)
             cccGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATAGCGCTTTCCCCGGATGaa                                                           >  1:3323441/1‑76 (MQ=255)
                gCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGaaaaa                                                        >  1:2754116/1‑76 (MQ=255)
                gCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGaaaaa                                                        <  2:6236398/76‑1 (MQ=255)
                 cTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGaaaaaa                                                       <  1:5771462/76‑1 (MQ=255)
                 cTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGaaaaaa                                                       <  2:5387487/76‑1 (MQ=255)
                 cTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGaaaaaa                                                       <  2:631583/76‑1 (MQ=255)
                 cTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGaaaaaa                                                       >  2:5086843/1‑76 (MQ=255)
                  tGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAat                                                      >  2:2234344/1‑76 (MQ=255)
                  tGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAat                                                      <  1:5754877/76‑1 (MQ=255)
                  tGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAat                                                      <  1:3664450/76‑1 (MQ=255)
                  tGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAat                                                      <  2:2536811/76‑1 (MQ=255)
                  tGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAat                                                      <  2:2668826/76‑1 (MQ=255)
                    gagaAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAatat                                                    >  1:2279809/1‑76 (MQ=255)
                     agaAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATAtt                                                   >  1:3627602/1‑76 (MQ=255)
                       aaaTACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATAttgt                                                 <  1:3221140/76‑1 (MQ=255)
                       aaaTACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATAttgt                                                 <  1:2002494/76‑1 (MQ=255)
                       aaaTACGACGGGACGGGGACGAGAAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGCTGAAAAAATAttgt                                                 >  1:2586277/1‑76 (MQ=255)
                         aTACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATAttgttg                                               <  1:1766599/76‑1 (MQ=255)
                          tACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGa                                              >  2:4711591/1‑76 (MQ=255)
                          tACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGa                                              <  1:5366746/76‑1 (MQ=255)
                          tACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGa                                              <  2:4472090/76‑1 (MQ=255)
                           acgacgGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGAt                                             <  1:544392/76‑1 (MQ=255)
                             gacgggacggGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGTTTTCCCCGGATGGAAAAATATTGTTGATaa                                           >  1:1145799/1‑75 (MQ=255)
                             gacgggacggGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATAc                                           >  1:4851380/1‑76 (MQ=255)
                              acgggacgggGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATAcc                                          <  1:3899494/76‑1 (MQ=255)
                              acgggacgggGACGAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATAcc                                          <  1:103541/76‑1 (MQ=255)
                                gggacggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCg                                        >  1:2310657/1‑76 (MQ=255)
                                 ggacggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCACTTTCCCCGGATGAAAAAATATTGTTGATACCCgg                                       >  1:6165397/1‑76 (MQ=255)
                                   acggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTa                                     >  2:2450393/1‑76 (MQ=255)
                                   acggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTa                                     <  2:6101953/76‑1 (MQ=255)
                                    cggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAGTATTGTTGATACCCGGTAt                                    <  2:4899789/76‑1 (MQ=255)
                                    cggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTAt                                    <  2:2741452/76‑1 (MQ=255)
                                    cggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTAt                                    <  1:6247029/76‑1 (MQ=255)
                                    cggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTAt                                    <  1:2070978/76‑1 (MQ=255)
                                     ggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATg                                   <  1:111279/76‑1 (MQ=255)
                                     ggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATg                                   <  2:1640164/76‑1 (MQ=255)
                                     ggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATg                                   <  2:5505273/76‑1 (MQ=255)
                                     ggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATg                                   <  1:5107530/76‑1 (MQ=255)
                                     ggggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAACATTGTTGATACCCGGTATg                                   >  2:2754659/1‑76 (MQ=255)
                                       ggacgAGTAATCAACCTGCTGCCTGGGTTATGGCGCGGCCCCCGCATGAAAAAATATTGGTGATACCCGGTATGGt                                 <  2:2894495/76‑1 (MQ=255)
                                       ggacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGt                                 >  2:3645280/1‑76 (MQ=255)
                                        gacgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTc                                <  1:1678976/76‑1 (MQ=255)
                                         acgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCa                               <  1:5761865/76‑1 (MQ=255)
                                         acgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCa                               >  1:1609842/1‑76 (MQ=255)
                                         acgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCa                               <  1:3810527/76‑1 (MQ=255)
                                         acgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCa                               <  1:2673006/76‑1 (MQ=255)
                                          cgAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCag                              <  2:1401050/76‑1 (MQ=255)
                                           gAGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCaga                             <  2:1782909/76‑1 (MQ=255)
                                            aGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCagag                            >  2:2189521/1‑76 (MQ=255)
                                            aGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCagag                            <  2:4265117/76‑1 (MQ=255)
                                            aGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCagag                            >  1:4411557/1‑76 (MQ=255)
                                            aGTAATAAACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCagag                            >  1:1969817/1‑76 (MQ=255)
                                               aataaACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTg                         <  2:4097606/76‑1 (MQ=255)
                                               aataaACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTg                         >  2:5307775/1‑76 (MQ=255)
                                               aataaACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTg                         <  1:1744637/76‑1 (MQ=255)
                                               aataaACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTg                         >  2:1988292/1‑76 (MQ=255)
                                                ataaACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGa                        <  2:5178193/76‑1 (MQ=255)
                                                ataaACCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGa                        <  1:2220986/76‑1 (MQ=255)
                                                  aaaCCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACa                      <  1:5090646/76‑1 (MQ=255)
                                                    aCCTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACAtt                    <  2:1949805/76‑1 (MQ=255)
                                                     ccTGGTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTg                   >  1:4878433/1‑76 (MQ=255)
                                                        ggTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTAc                <  1:2636216/76‑1 (MQ=255)
                                                        ggTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTAc                <  1:5081452/76‑1 (MQ=255)
                                                         gTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACg               <  1:3727982/76‑1 (MQ=255)
                                                         gTGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACg               >  2:5261198/1‑76 (MQ=255)
                                                          tGCCTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGc              <  1:3455262/76‑1 (MQ=255)
                                                            ccTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCAt            <  2:881736/76‑1 (MQ=255)
                                                            ccTGGGTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCAt            <  2:4026675/76‑1 (MQ=255)
                                                               gggTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATGCCCGGTATGGTCAGAGGTGACATTGTACGCATCaa         <  2:3119570/76‑1 (MQ=255)
                                                               gggTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCaa         <  2:4746491/76‑1 (MQ=255)
                                                                ggTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAAt        >  2:3321012/1‑76 (MQ=255)
                                                                ggTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAAt        <  2:1812298/76‑1 (MQ=255)
                                                                ggTTATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAAt        <  1:3768245/76‑1 (MQ=255)
                                                                   tATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAATAcc     <  1:5160376/76‑1 (MQ=255)
                                                                   tATGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAATAcc     <  2:3694970/76‑1 (MQ=255)
                                                                    aTGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAATACCa    >  2:2847848/1‑76 (MQ=255)
                                                                     tGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAATACCAt   <  1:6072124/76‑1 (MQ=255)
                                                                     tGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAATACCAt   <  1:4709314/76‑1 (MQ=255)
                                                                     tGGCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATCGTCAGAGGTGACATTGTACGCATCAATACCAt   <  2:1583683/76‑1 (MQ=255)
                                                                      ggCGCTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAATACCATc  <  1:5030261/76‑1 (MQ=255)
                                                                         |                                                                        
ATGTGGCGGATATCCCGCTGGAGAAATACGACGGGACGGGGACGAGTAATAAACCTGGTGCCTGGGTTATGGCACTTTCCCCGGATGAAAAAATATTGTTGATACCCGGTATGGTCAGAGGTGACATTGTACGCATCAATACCATC  >  REL606/670716‑670861

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