Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,029,249 Δ1 bp intergenic (‑788/+273) ECB_02828 ← / ← yghD Polysialic acid transport protein kpsM/predicted secretion pathway M‑type protein, membrane anchored

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,029,2460T.93.8% 345.9 / 18.1 81intergenic (‑785/+276)ECB_02828/yghDPolysialic acid transport protein kpsM/predicted secretion pathway M‑type protein, membrane anchored
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base . (25/51);  total (25/56)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGATTATCTCCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACTTAAG  >  REL606/3029174‑3029317
                                                                        |                                                                       
tGATTATCTCCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGtttg                                                                      >  2:555741/1‑75 (MQ=255)
 gATTATCTCCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGtttgg                                                                     >  2:1465990/1‑74 (MQ=255)
  aTTATCTCCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGtttggt                                                                    >  1:1492910/1‑73 (MQ=255)
  aTTATCTCCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGtttggt                                                                    >  1:1492902/1‑73 (MQ=255)
    tATCTCCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCt                                                                 <  1:25958/76‑1 (MQ=255)
      tctcCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTAt                                                               <  1:607675/76‑1 (MQ=255)
        tcCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCa                                                             >  2:66804/1‑76 (MQ=255)
        tcCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCa                                                             >  2:254749/1‑76 (MQ=255)
        tcCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCa                                                             >  1:85955/1‑76 (MQ=255)
            ggAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGATTTATCCTATATGTTTGGTCTATCAATTg                                                         <  1:627816/76‑1 (MQ=255)
            ggAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTg                                                         <  1:627814/76‑1 (MQ=255)
              aaTTTTACTAATTTCCAAGCTACACAAAATACAGATCTCACAAGAATTAGCCTGNATGTTTGGTCTATCAATTGtt                                                       >  2:278310/1‑76 (MQ=255)
              aaTTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGtt                                                       >  1:779676/1‑76 (MQ=255)
              aaTTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGtt                                                       >  2:1494718/1‑76 (MQ=255)
                 tttACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACt                                                    >  1:998055/1‑76 (MQ=255)
                  ttACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTc                                                   <  2:1475173/76‑1 (MQ=255)
                  ttACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTc                                                   >  2:1327293/1‑76 (MQ=255)
                   tACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCa                                                  <  1:108351/76‑1 (MQ=255)
                   tACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCa                                                  <  2:488626/76‑1 (MQ=255)
                   tACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCa                                                  <  2:348731/76‑1 (MQ=255)
                   tACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCa                                                  <  2:348728/76‑1 (MQ=255)
                   tACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCa                                                  <  2:348699/76‑1 (MQ=255)
                   tACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCa                                                  <  1:849314/76‑1 (MQ=255)
                   tACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCa                                                  <  2:609951/76‑1 (MQ=255)
                   tACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCa                                                  <  2:760836/76‑1 (MQ=255)
                    aCTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCaa                                                 <  1:246251/76‑1 (MQ=255)
                    aCTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCaa                                                 <  1:246238/76‑1 (MQ=255)
                    aCTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCaa                                                 >  2:352777/1‑76 (MQ=255)
                       aaTTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGa                                              >  2:823024/1‑76 (MQ=255)
                         tttCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGAtt                                            >  2:806550/1‑76 (MQ=255)
                         tttCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGAtt                                            >  1:1339539/1‑76 (MQ=255)
                          ttCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGAttt                                           <  2:1642933/76‑1 (MQ=255)
                          ttCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGAttt                                           <  1:802303/76‑1 (MQ=255)
                           tCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATACGTTTTGGTCTATCAATTGTTACTCAATGAttt                                           <  2:598517/76‑1 (MQ=255)
                               aGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGc                                      <  2:670985/76‑1 (MQ=255)
                                   acacAAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACt                                  <  1:1292077/76‑1 (MQ=255)
                                    cacaAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTa                                 <  1:426077/76‑1 (MQ=255)
                                     acaAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTAt                                <  2:741698/76‑1 (MQ=255)
                                     acaAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTAt                                <  1:1327586/76‑1 (MQ=255)
                                     acaAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTAt                                <  2:50470/76‑1 (MQ=255)
                                     acaAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTAt                                <  2:50450/76‑1 (MQ=255)
                                      gaaaaTACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTAtt                               <  1:1085851/75‑1 (MQ=255)
                                      caAAATACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTAtt                               >  2:333168/1‑76 (MQ=255)
                                       aaaaTACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATta                              <  2:1529967/76‑1 (MQ=255)
                                        aaaTACAGATCACACAAGAATTATCCTATATGTTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATta                              <  1:980641/76‑1 (MQ=255)
                                        aaaTACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATtaa                             <  2:627582/76‑1 (MQ=255)
                                         aaTACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATtaat                            >  2:83302/1‑76 (MQ=255)
                                          aTACAGCTCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATtaata                           >  2:993922/1‑76 (MQ=255)
                                          aTACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATtaata                           >  2:94480/1‑76 (MQ=255)
                                          aTACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATtaata                           >  2:1086160/1‑76 (MQ=255)
                                          aTACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATtaata                           <  1:1251108/76‑1 (MQ=255)
                                           tACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATtaataa                          <  1:1606791/76‑1 (MQ=255)
                                           tACAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATtaataa                          <  2:769204/76‑1 (MQ=255)
                                            aCAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTaataat                         <  2:97320/76‑1 (MQ=255)
                                            aCAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTaataat                         >  2:331879/1‑76 (MQ=255)
                                             cAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATc                        <  1:528644/76‑1 (MQ=255)
                                             cAGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATc                        <  1:528649/76‑1 (MQ=255)
                                              aGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCa                       <  1:147124/76‑1 (MQ=255)
                                              aGATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCa                       <  1:147109/76‑1 (MQ=255)
                                               gATCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCAt                      <  1:249834/76‑1 (MQ=255)
                                                 tCACACAAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATaa                    <  2:731580/76‑1 (MQ=255)
                                                 ncacacaAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATaa                    >  1:145533/2‑76 (MQ=255)
                                                  cacacaAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAAt                   >  1:1671797/1‑76 (MQ=255)
                                                  cacacaAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAAt                   <  2:1689291/76‑1 (MQ=255)
                                                   acacaAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAAtg                  >  1:814535/1‑76 (MQ=255)
                                                   acacaAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAAtg                  <  1:150980/76‑1 (MQ=255)
                                                     acaAGAATTATCCTATATGTTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAAtgt                 <  2:352058/76‑1 (MQ=255)
                                                     acaAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAAtgtg                <  1:875044/76‑1 (MQ=255)
                                                     acaAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAAtgtg                <  1:875061/76‑1 (MQ=255)
                                                      caAGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATgtgt               <  2:1139460/76‑1 (MQ=255)
                                                       aaGAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTc              >  1:683086/1‑76 (MQ=255)
                                                         gAATTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCaa            <  1:632585/76‑1 (MQ=255)
                                                          aaTTATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCaaa           <  1:1450380/76‑1 (MQ=255)
                                                             tATCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAaca        >  2:1318977/1‑76 (MQ=255)
                                                             tATCCTATATATTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAac         <  1:710973/76‑1 (MQ=255)
                                                               tCCTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACt      <  1:563227/76‑1 (MQ=255)
                                                                ccTATATGTTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACt      <  2:833969/76‑1 (MQ=255)
                                                                ccTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACtt     >  1:307878/1‑76 (MQ=255)
                                                                ccTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACtt     <  1:705304/76‑1 (MQ=255)
                                                                ccTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACtt     <  2:773416/76‑1 (MQ=255)
                                                                ccTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACtt     <  2:1234082/76‑1 (MQ=255)
                                                                 cTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATANTCATAATGTGTCAAAACACTTa    <  2:230505/76‑1 (MQ=255)
                                                                 cTATATGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACTTa    <  1:1042839/76‑1 (MQ=255)
                                                                   atatGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACTTAAg  <  1:1299661/76‑1 (MQ=255)
                                                                   atatGTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACTTAAg  >  2:1394366/1‑76 (MQ=255)
                                                                        |                                                                       
TGATTATCTCCAGGAATTTTACTAATTTCCAAGCTACACAAAATACAGATCACACAAGAATTATCCTATATGTTTTGGTCTATCAATTGTTACTCAATGATTTATTGCAACTATTAATAATCATAATGTGTCAAAACACTTAAG  >  REL606/3029174‑3029317

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: