Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,027,164 (A)6→5 coding (507/660 nt) ECB_02827 ← Polysialic acid transport ATP‑binding KpsT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,027,1590A.95.9% 694.2 / 21.5 146coding (512/660 nt)ECB_02827Polysialic acid transport ATP‑binding KpsT
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base . (81/59);  total (83/63)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 168 total aligned reads displayed.

ACCGCGTTATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCACCGACTG  >  REL606/3027084‑3027231
                                                                           |                                                                        
aCCGCGTTATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTa                                                                          <  2:3081674/76‑1 (MQ=255)
aCCGCGTTATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTa                                                                          <  1:852213/76‑1 (MQ=255)
 ccGCGTTATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTaa                                                                         <  1:2103063/76‑1 (MQ=255)
    cgTTATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTaaaaa                                                                      <  1:5256383/76‑1 (MQ=255)
      ttATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTaaaaact                                                                    >  1:501685/1‑74 (MQ=255)
      ttATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTaaaaact                                                                    <  1:5871711/76‑3 (MQ=255)
      ttATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTaaaaact                                                                    >  1:6972551/1‑74 (MQ=255)
      ttATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTaaaaact                                                                    >  2:225087/1‑74 (MQ=255)
       tATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTaaaaacta                                                                   >  2:1835078/1‑73 (MQ=255)
       tATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATCCAAACTATGTGAAACCATTaaaaacta                                                                   >  2:5219596/1‑73 (MQ=255)
           aTTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGAtt                                                              >  1:1752199/1‑76 (MQ=255)
           aTTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGAtt                                                              <  1:661210/76‑1 (MQ=255)
           aTTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGAtt                                                              <  1:3529445/76‑1 (MQ=255)
           aTTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGAtt                                                              <  1:5943175/76‑1 (MQ=255)
             ttttAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCt                                                            <  1:3038305/76‑1 (MQ=255)
              tttAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCtt                                                           <  2:7233639/76‑1 (MQ=255)
               ttAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCttt                                                          >  1:6335743/1‑76 (MQ=255)
               ttAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCttt                                                          >  1:6042181/1‑76 (MQ=255)
               ttAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCttt                                                          >  1:4912005/1‑76 (MQ=255)
                tAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTa                                                         >  1:6832602/1‑76 (MQ=255)
                 aaaTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTAt                                                        <  2:7089334/76‑1 (MQ=255)
                 aaaTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTAt                                                        >  2:5435000/1‑76 (MQ=255)
                 aaaTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTAt                                                        <  1:1134192/76‑1 (MQ=255)
                 aaaTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTAt                                                        >  1:1402641/1‑76 (MQ=255)
                 aaaTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTAt                                                        >  1:2179664/1‑76 (MQ=255)
                  aaTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATg                                                       >  2:4241686/1‑76 (MQ=255)
                  aaTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATg                                                       >  2:4397120/1‑76 (MQ=255)
                   aTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGc                                                      >  2:2700937/1‑76 (MQ=255)
                   aTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGc                                                      >  2:6447183/1‑76 (MQ=255)
                   aTACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGc                                                      >  1:7506529/1‑76 (MQ=255)
                     aCAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCct                                                    <  1:6922577/76‑1 (MQ=255)
                     aCAAAGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCct                                                    >  1:5079954/1‑76 (MQ=255)
                      cAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctt                                                   <  2:6987992/76‑1 (MQ=255)
                      cAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctt                                                   >  1:1586668/1‑76 (MQ=255)
                      cAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctt                                                   >  1:3475078/1‑76 (MQ=255)
                      cAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctt                                                   >  1:329117/1‑76 (MQ=255)
                          ggCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctttctt                                               >  1:2266732/1‑76 (MQ=255)
                          ggCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctttctt                                               >  1:4483795/1‑76 (MQ=255)
                           gCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctttcttt                                              <  1:3116442/76‑1 (MQ=255)
                           gCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctttcttt                                              <  2:744120/76‑1 (MQ=255)
                           gCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctttcttt                                              <  2:4849979/76‑1 (MQ=255)
                           gCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCctttcttt                                              <  2:2333286/76‑1 (MQ=255)
                             cacaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaa                                            >  2:5913769/1‑76 (MQ=255)
                             cacaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaa                                            <  1:2832092/76‑1 (MQ=255)
                             cacaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaa                                            >  2:1241057/1‑76 (MQ=255)
                             cacaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaa                                            <  2:4793364/76‑1 (MQ=255)
                              acaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaaa                                           >  1:6694491/1‑76 (MQ=255)
                              acaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaaa                                           <  2:3213671/76‑1 (MQ=255)
                              acaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaaa                                           <  2:2665324/76‑1 (MQ=255)
                              acaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaaa                                           <  2:4406166/76‑1 (MQ=255)
                              acaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaaa                                           <  1:5932742/76‑1 (MQ=255)
                              acaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaaa                                           <  2:4529008/76‑1 (MQ=255)
                              acaTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTaaa                                           >  1:7395617/1‑76 (MQ=255)
                               caTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAAc                                          <  2:5083541/76‑1 (MQ=255)
                               caTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAAc                                          >  2:6627761/1‑76 (MQ=255)
                                aTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACa                                         >  2:2090003/1‑76 (MQ=255)
                                aTCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACa                                         <  2:5111184/76‑1 (MQ=255)
                                  cGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAAt                                       >  2:5350775/1‑76 (MQ=255)
                                  cGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAAt                                       >  1:3225419/1‑76 (MQ=255)
                                   gCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAAtt                                      <  1:2062286/76‑1 (MQ=255)
                                   gCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAAtt                                      >  1:97422/1‑76 (MQ=255)
                                    cAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTg                                     <  2:50214/76‑1 (MQ=255)
                                      aaaCTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAg                                   <  1:4020977/76‑1 (MQ=255)
                                      aaaCTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAg                                   >  2:1293403/1‑76 (MQ=255)
                                        aCTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAgcg                                 >  1:5013188/1‑76 (MQ=255)
                                          tctTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCa                               >  2:6594402/1‑76 (MQ=255)
                                          tctTTCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCa                               >  1:95118/1‑76 (MQ=255)
                                            tttCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCAtt                             >  2:3527318/1‑76 (MQ=255)
                                              tCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCAtttt                           <  2:6257056/76‑1 (MQ=255)
                                              tCAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCAtttt                           <  1:6033117/76‑1 (MQ=255)
                                               cAATGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCAttttt                          <  2:4640446/76‑1 (MQ=255)
                                                aaTGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTc                         >  2:6237492/1‑76 (MQ=255)
                                                aaTGAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTc                         >  2:1659653/1‑76 (MQ=255)
                                                   gAATTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCttt                      >  2:855305/1‑76 (MQ=255)
                                                    aaTTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTa                     <  1:630819/76‑1 (MQ=255)
                                                    aaTTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTa                     <  2:1501532/76‑1 (MQ=255)
                                                     aTTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTaa                    >  1:4151539/1‑76 (MQ=255)
                                                     aTTCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTaa                    >  1:5743788/1‑76 (MQ=255)
                                                      ttCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTaaa                   >  1:5524204/1‑76 (MQ=255)
                                                       tCAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAAc                  >  2:3041722/1‑76 (MQ=255)
                                                        cAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAAcc                 >  1:67905/1‑76 (MQ=255)
                                                        cAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAAcc                 <  2:5789676/76‑1 (MQ=255)
                                                        cAAACTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAAcc                 >  2:5298015/1‑76 (MQ=255)
                                                         aaaCTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCt                >  2:3857681/1‑76 (MQ=255)
                                                         aaaCTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCt                >  2:2554798/1‑76 (MQ=255)
                                                          aaCTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTg               >  1:6203290/1‑76 (MQ=255)
                                                          aaCTATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTg               <  2:2847412/76‑1 (MQ=255)
                                                             tATGTGAAACCATTAAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGc             >  2:5478355/1‑76 (MQ=255)
                                                             tATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCa            >  2:2790273/1‑76 (MQ=255)
                                                             tATGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCa            >  1:6450079/1‑76 (MQ=255)
                                                             tATGTGAAACCATTAAAAACTAGATTCTTTAGGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCa            >  2:5805140/1‑76 (MQ=255)
                                                              aTGTGAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCAt           <  1:1250876/76‑1 (MQ=255)
                                                                gtgAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCa         >  2:4017279/1‑76 (MQ=255)
                                                                  gAAACCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCAcc       >  2:5508415/1‑76 (MQ=255)
                                                                   aaaCCATTAAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCAcc       <  2:1716487/76‑1 (MQ=255)
                                                                    aaCCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCACCGa     <  2:1499975/76‑1 (MQ=255)
                                                                     aCCATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCACCGAc    <  1:2843386/76‑1 (MQ=255)
                                                                      ccATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCACCGACt   <  2:5566628/76‑1 (MQ=255)
                                                                      ccATTAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCACCGACt   <  1:5347097/76‑1 (MQ=255)
                                                                        aTTAAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCACCGACTg  >  2:3437921/1‑76 (MQ=255)
                                                                           |                                                                        
ACCGCGTTATCATTTTTAAATACAATGGCCACATCGCAAAACTCTTTCAATGAATTCAAACTATGTGAAACCATTAAAAAACTAGATTCTTTATGCCTTTCTTTAAACAATTGAGCGCATTTTTCTTTAAACCTGGCATCACCGACTG  >  REL606/3027084‑3027231

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