Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,028,245 Δ1 bp coding (217/777 nt) ECB_02828 ← Polysialic acid transport protein kpsM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,028,2430G.93.1% 232.8 / 18.0 61coding (219/777 nt)ECB_02828Polysialic acid transport protein kpsM
Reads supporting (aligned to +/- strand):  ref base G (1/2);  major base . (25/29);  minor base A (4/0);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

CGATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAACCAAA  >  REL606/3028169‑3028314
                                                                          |                                                                       
cGATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAgg                                                                        >  1:1120431/1‑76 (MQ=255)
cGATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAgg                                                                        >  1:112766/1‑76 (MQ=255)
  aTAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAggaa                                                                      >  1:1315365/1‑74 (MQ=255)
     aTTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGGAATc                                                                   <  1:199607/76‑1 (MQ=255)
     aTTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCa                                                                  <  1:1196315/76‑1 (MQ=255)
      ttAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAg                                                                 >  1:237962/1‑76 (MQ=255)
      ttAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAg                                                                 >  1:268822/1‑76 (MQ=255)
      ttAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAg                                                                 >  1:945075/1‑76 (MQ=255)
        aaaCAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGGAATCAgg                                                                >  2:260425/1‑76 (MQ=255)
          acaacCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCa                                                             <  2:1057293/76‑1 (MQ=255)
           caacCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCAt                                                            >  2:1244119/1‑76 (MQ=255)
            aacCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCAtt                                                           >  2:1409684/1‑76 (MQ=255)
             acCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTa                                                          <  1:580464/76‑1 (MQ=255)
              ccccTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTaa                                                         <  2:1184870/76‑1 (MQ=255)
              ccccTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTaa                                                         >  2:127217/1‑76 (MQ=255)
              ccccTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTaa                                                         >  2:1406557/1‑76 (MQ=255)
              ccccTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTaa                                                         >  2:1443498/1‑76 (MQ=255)
              ccccTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTaa                                                         >  2:996889/1‑76 (MQ=255)
              ccccTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTaa                                                         <  1:1296459/76‑1 (MQ=255)
                ccTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGt                                                       <  2:497249/76‑1 (MQ=255)
                 cTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTa                                                      <  2:704044/76‑1 (MQ=255)
                   gATTCGCTTCAATAGCACCTATAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTaaa                                                    >  1:811991/1‑76 (MQ=255)
                     ttCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTaaaaa                                                  <  1:51943/76‑1 (MQ=255)
                        gCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAacac                                               <  1:1086857/76‑1 (MQ=255)
                          ttCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCg                                             >  1:1573921/1‑76 (MQ=255)
                            cAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGGAATCAGGCCATTAAGTAAAAACACCgg                                            <  1:137524/76‑1 (MQ=255)
                            cAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCggg                                           >  1:352426/1‑76 (MQ=255)
                              aTAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGaa                                         <  2:1250476/76‑1 (MQ=255)
                               tAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAAt                                        <  2:333109/76‑1 (MQ=255)
                                aGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATg                                       <  2:743149/76‑1 (MQ=255)
                                aGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATg                                       <  1:553512/76‑1 (MQ=255)
                                 gCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGaa                                         >  1:681911/1‑73 (MQ=255)
                                 gCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGaa                                         <  2:681911/73‑1 (MQ=255)
                                   aCCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATgaga                                    >  1:1512532/1‑76 (MQ=255)
                                     cTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATg                                  <  2:1606243/76‑1 (MQ=255)
                                     cTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATg                                  <  2:959383/76‑1 (MQ=255)
                                       aCAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTc                                >  2:1483515/1‑76 (MQ=255)
                                        cAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCt                               >  2:1148798/1‑76 (MQ=255)
                                         aGAACGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTg                              <  2:785099/76‑1 (MQ=255)
                                           aaCGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGc                            <  1:1267563/76‑1 (MQ=255)
                                           aaCGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGc                            <  1:19557/76‑1 (MQ=255)
                                           aaCGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGc                            <  1:224411/76‑1 (MQ=255)
                                             cGATTGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCAt                          >  2:1237195/1‑76 (MQ=255)
                                                 tGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTg                      <  2:789332/76‑1 (MQ=255)
                                                 tGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTg                      <  2:246124/76‑1 (MQ=255)
                                                 tGCTGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTg                      >  1:1117560/1‑76 (MQ=255)
                                                    tGATACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCgg                   >  1:1095376/1‑76 (MQ=255)
                                                       tACTGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGc                <  1:455642/76‑1 (MQ=255)
                                                          tGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTATAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATa             <  2:1421117/76‑1 (MQ=255)
                                                          tGCTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATa             <  1:1072184/76‑1 (MQ=255)
                                                            cTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGcatcat           >  1:272899/1‑76 (MQ=255)
                                                            cTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAAt           >  1:272914/1‑76 (MQ=255)
                                                            cTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAAt           >  1:1198994/1‑76 (MQ=255)
                                                            cTAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAAt           <  2:613476/76‑1 (MQ=255)
                                                             tAAAGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATg          <  2:55636/76‑1 (MQ=255)
                                                              aaaGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGt         >  2:537796/1‑76 (MQ=255)
                                                              aaaGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGt         <  2:268822/76‑1 (MQ=255)
                                                              aaaGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGt         >  2:1534424/1‑76 (MQ=255)
                                                              aaaGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGt         >  2:1115741/1‑76 (MQ=255)
                                                                aGATAAAAAAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTaa       <  2:127363/76‑1 (MQ=255)
                                                                   ataaaaaaGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAAcc     <  2:462068/74‑1 (MQ=255)
                                                                   ataaaaaaGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAAcc     <  1:1054805/74‑1 (MQ=255)
                                                                   ataaaaaaGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAAcc     <  1:1647008/74‑1 (MQ=255)
                                                                   ataaaaaaGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAAcc     <  1:1054814/74‑1 (MQ=255)
                                                                   ataaaaaaGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAAcc     <  1:524131/74‑1 (MQ=255)
                                                                   ataaaaaaGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAAcc     >  1:825842/3‑76 (MQ=255)
                                                                    taaaaaaGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAACCa    <  1:433079/75‑1 (MQ=255)
                                                                    taaaaaaGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAACCa    >  2:779789/2‑76 (MQ=255)
                                                                     aaaaaaGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAACCaa   >  1:524043/1‑76 (MQ=255)
                                                                      aaaaAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAACCaaa  >  2:689339/1‑76 (MQ=255)
                                                                      aaaaAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAACCaaa  >  2:640516/1‑76 (MQ=255)
                                                                      aaaaAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAACCaaa  >  2:760760/1‑76 (MQ=255)
                                                                      aaaaAGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAACCaaa  >  2:552595/1‑76 (MQ=255)
                                                                          |                                                                       
CGATAATTAAACAACCCCTGATTCGCTTCAATAGCACCTACAGAACGATTGCTGATACTGCTAAAGATAAAAAAGGGAATCAGGCCATTAAGTAAAAACACCGGGAATGAGATGTCTGGCATCGTGCGGTGCATAATGTAACCAAA  >  REL606/3028169‑3028314

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