Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,015,124 G→A G182D (GGC→GAC)  ECB_02816 → KpsD protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,015,1240GA99.0% 1200.5 / 1.2 290G182D (GGC→GAC) ECB_02816KpsD protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (131/156);  total (132/158)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

Only 100 of 291 total aligned reads displayed.

TACGCCTCCTTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGTTA  >  REL606/3015051‑3015196
                                                                         |                                                                        
tACGCCTCCTTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACg                                                                        <  1:2806644/76‑1 (MQ=255)
tACGCCTCCTTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACg                                                                        >  1:4306794/1‑76 (MQ=255)
 aCGCCTCCTTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACgg                                                                       >  1:1036537/1‑76 (MQ=255)
   gCCTCCTTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGgtg                                                                     >  2:3723898/1‑76 (MQ=255)
    cctcctTATTCCGGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGgtgt                                                                    <  1:3411755/76‑1 (MQ=255)
    cctcctTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGgtgt                                                                    >  1:6729964/1‑76 (MQ=255)
     ctcctTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGtgtg                                                                   <  1:3554069/76‑1 (MQ=255)
      tcctTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGa                                                                  <  1:7771800/76‑1 (MQ=255)
       cctTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGAc                                                                 >  2:7672115/1‑76 (MQ=255)
         ttattaCAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGt                                                               >  2:2181138/1‑76 (MQ=255)
          tattaCAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTc                                                              >  2:6529247/1‑76 (MQ=255)
          tattaCAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTc                                                              >  2:5015754/1‑76 (MQ=255)
          tattaCAGGCCCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTc                                                              <  1:7628796/76‑1 (MQ=255)
           attaCAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCt                                                             >  2:7584081/1‑76 (MQ=255)
           attaCAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCt                                                             >  2:201158/1‑76 (MQ=255)
            ttaCGGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTg                                                            <  2:6355120/76‑1 (MQ=255)
            ttaCATGCGCAGCCAGGAAAAGTGTACGAGACCGGATTTGTGCGGCATCCTGGTCTGTATGACGGAGTGACGTCTg                                                            >  2:3217483/1‑76 (MQ=255)
            ttaCAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTg                                                            <  2:4749286/76‑1 (MQ=255)
               cAGCCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGAtt                                                         <  1:5320342/76‑1 (MQ=255)
                aGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTc                                                        <  1:2008224/76‑1 (MQ=255)
                    gcAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTa                                                    >  1:1629546/1‑76 (MQ=255)
                    ccAGCCAGTAAAAGTGTACGTGTCCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTa                                                    <  1:3123917/75‑1 (MQ=255)
                     cAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTAc                                                   >  1:4470158/1‑76 (MQ=255)
                       gCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTc                                                 >  1:5946338/1‑76 (MQ=255)
                        ccAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCa                                                >  1:6863611/1‑76 (MQ=255)
                        ccAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCa                                                <  1:7558209/76‑1 (MQ=255)
                        ccAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCa                                                <  1:1141918/76‑1 (MQ=255)
                         cAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCaa                                               <  1:1564822/76‑1 (MQ=255)
                          aGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAAt                                              >  2:6239482/1‑76 (MQ=255)
                           gTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAAtt                                             <  2:5215605/76‑1 (MQ=255)
                            tAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTa                                            >  2:3462080/1‑76 (MQ=255)
                             aaaaGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTAt                                           >  1:6342260/1‑76 (MQ=255)
                             aaaaGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTAt                                           >  1:4584962/1‑76 (MQ=255)
                             aaaaGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTAt                                           >  2:2774627/1‑76 (MQ=255)
                             aaaaGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTAt                                           <  2:4954913/76‑1 (MQ=255)
                                  tgtACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGAt                                      <  2:2631896/76‑1 (MQ=255)
                                  tgtACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGAt                                      <  2:4037615/76‑1 (MQ=255)
                                  tgtACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGAt                                      >  2:956802/1‑76 (MQ=255)
                                   gtACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATc                                     >  2:2439346/1‑76 (MQ=255)
                                   gtACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATc                                     >  1:2594442/1‑76 (MQ=255)
                                   gtACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATc                                     >  2:7465476/1‑76 (MQ=255)
                                   gtACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATc                                     <  2:7008335/76‑1 (MQ=255)
                                   gtACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATc                                     >  2:6888150/1‑76 (MQ=255)
                                    tACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCa                                    <  1:6358134/76‑1 (MQ=255)
                                    tACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCa                                    >  1:2216842/1‑76 (MQ=255)
                                    tACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCa                                    <  1:574883/76‑1 (MQ=255)
                                     aCGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCATTACTCAATTATCTGATCaa                                   <  1:5011216/76‑1 (MQ=255)
                                      cGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAg                                  <  1:7181447/76‑1 (MQ=255)
                                      cGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAg                                  <  2:835667/76‑1 (MQ=255)
                                        tGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGc                                >  1:2972611/1‑76 (MQ=255)
                                        tGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGc                                >  2:3844604/1‑76 (MQ=255)
                                        tGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGc                                >  2:3126139/1‑76 (MQ=255)
                                        tGACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGc                                <  1:4055618/76‑1 (MQ=255)
                                         gACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCt                               >  1:5672150/1‑76 (MQ=255)
                                         gACCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCt                               <  2:1284331/76‑1 (MQ=255)
                                          aCCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTg                              <  2:7879716/76‑1 (MQ=255)
                                          aCCGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTg                              >  2:4366367/1‑76 (MQ=255)
                                           ccGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTgg                             >  1:1420573/1‑76 (MQ=255)
                                            cGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTggg                            >  2:795077/1‑75 (MQ=255)
                                            cGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTggg                            >  2:4624546/1‑75 (MQ=255)
                                            cGGATTTGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTggc                            <  1:232742/76‑1 (MQ=255)
                                                tttGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTggcgggg                        >  1:5788071/1‑74 (MQ=255)
                                                tttGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGgcggcg                        >  1:196322/1‑76 (MQ=255)
                                                 ttGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTgggggggg                       >  2:5497927/1‑73 (MQ=255)
                                                 ttGTGCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGt                       >  2:7502137/1‑76 (MQ=255)
                                                   gtgCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTg                     <  2:1385590/76‑1 (MQ=255)
                                                    tgCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGa                    <  2:6041029/76‑1 (MQ=255)
                                                     gCGTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGAt                   >  1:5592201/1‑76 (MQ=255)
                                                       gTAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATcc                 >  1:3944908/1‑76 (MQ=255)
                                                        tAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGGGGGGTTGATTCa                >  1:6627696/1‑76 (MQ=255)
                                                        tAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCa                <  1:5189043/76‑1 (MQ=255)
                                                        tAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCa                <  1:370983/76‑1 (MQ=255)
                                                        tAATCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGGTGATTCa                >  2:3062773/1‑76 (MQ=255)
                                                         aaTCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATcccg               >  1:6153058/1‑74 (MQ=255)
                                                         aaTCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCag               >  1:4834256/1‑76 (MQ=255)
                                                         aaTCCTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCag               <  2:1162735/76‑1 (MQ=255)
                                                             cTGGTCTGTATGACGTTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGGGGTGGTCCagaggg           >  1:2079765/1‑74 (MQ=255)
                                                             cTGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGGGGGGGTGATCTAGAgcg           >  1:3254152/1‑76 (MQ=255)
                                                              tGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGGGGCGTTGATCCCGAgggc          >  2:6625924/1‑76 (MQ=255)
                                                              tGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTTATCCCGAgcgg          >  1:2780235/1‑75 (MQ=255)
                                                              tGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAgcgc          <  2:39206/76‑1 (MQ=255)
                                                              tGGTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAgcgc          <  1:3818527/76‑1 (MQ=255)
                                                                gTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCgg        <  2:3685807/76‑1 (MQ=255)
                                                                gTCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCgg        <  2:3608814/76‑1 (MQ=255)
                                                                 tCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGa       <  2:5058733/76‑1 (MQ=255)
                                                                 tCTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGa       >  2:7708095/1‑76 (MQ=255)
                                                                 tCTGTATGACGGTGTGAAGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCCGa       >  2:2910800/1‑76 (MQ=255)
                                                                  cTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGGTCCAGAGCGCggca      >  2:5514590/1‑74 (MQ=255)
                                                                  cTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGGGCggga      >  1:894746/1‑74 (MQ=255)
                                                                  cTGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGaa      <  2:2496762/76‑1 (MQ=255)
                                                                   tGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAg     <  1:3786498/76‑1 (MQ=255)
                                                                   tGTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAg     <  2:3386396/76‑1 (MQ=255)
                                                                    gTATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGt    <  2:7241898/76‑1 (MQ=255)
                                                                     tATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGGTCCCGAGCGCGGAAGtt   >  1:6229159/1‑76 (MQ=255)
                                                                     tATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGtt   <  1:7269395/76‑1 (MQ=255)
                                                                     tATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGtt   <  1:3014762/76‑1 (MQ=255)
                                                                     tATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGtt   <  2:2354553/76‑1 (MQ=255)
                                                                     tATGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGtt   <  1:6418956/76‑1 (MQ=255)
                                                                      aTGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGTTa  <  2:5918981/76‑1 (MQ=255)
                                                                      aTGACGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGTTa  <  2:2843655/76‑1 (MQ=255)
                                                                         |                                                                        
TACGCCTCCTTATTACAGGCGCAGCCAGTAAAAGTGTACGTGACCGGATTTGTGCGTAATCCTGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGTTA  >  REL606/3015051‑3015196

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