Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 975,083 T→A D255E (GAT→GAA)  mukF → chromosome condensin MukBEF, kleisin‑like subunit, binds calcium

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913975,0830TA90.3% 229.3 / 14.3 72D255E (GAT→GAA) mukFchromosome condensin MukBEF, kleisin‑like subunit, binds calcium
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base A (35/30);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGT  >  NC_000913/974988‑975176
                                                                                               |                                                                                             
gCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTgn                                                                                           <  2:735170/100‑2 (MQ=255)
  ggCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCag                                                                                         <  2:1064779/100‑1 (MQ=255)
    cAGGCGACAAATTGCAGGCTACTCTGTTGCGCCTTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACtcc                                                                                           >  1:143042/1‑94 (MQ=255)
     aGGCTACAAAGGTCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAAATGAAGAGc                                                                                      >  1:357084/1‑100 (MQ=255)
      ggCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCa                                                                                     <  1:515490/100‑1 (MQ=255)
       gCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGc                                                                                      >  1:1877978/1‑98 (MQ=255)
       gCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCaa                                                                                    >  2:1208035/1‑100 (MQ=255)
         gACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAAc                                                                                  <  2:598789/100‑1 (MQ=255)
         gACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAAc                                                                                  <  2:316011/100‑1 (MQ=255)
          aCAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACt                                                                                 >  1:65482/1‑100 (MQ=255)
             aaTTGCAGGCTAATCTGTTGCGCATTCAGGGTGCGACTATGACCCATGACGATCTGCATTTCGTCGATCGGCTGGTGTTCGAACTGCAGAGCAAACTCg                                                                               <  1:1158796/99‑1 (MQ=255)
             aaTTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGa                                                                              >  2:1195572/1‑100 (MQ=255)
              aTTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGAt                                                                             <  2:1526071/100‑1 (MQ=255)
              aTTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGAt                                                                             <  2:580965/100‑1 (MQ=255)
                tGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCg                                                                           <  1:1643552/100‑1 (MQ=255)
                 gCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGt                                                                          >  2:704866/1‑100 (MQ=255)
                    ggCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGtatt                                                                       >  1:1211115/1‑100 (MQ=255)
                    ggCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGtatt                                                                       >  1:204632/1‑100 (MQ=255)
                    ggCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGtatt                                                                       <  1:689404/100‑1 (MQ=255)
                       tAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATc                                                                    <  1:1491459/100‑1 (MQ=255)
                        aaTCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCa                                                                   >  1:1592862/1‑100 (MQ=255)
                        aaTCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCa                                                                   >  2:1849279/1‑100 (MQ=255)
                        aaTCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCa                                                                   >  2:1803877/1‑100 (MQ=255)
                        aaTCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCa                                                                   <  2:1592873/100‑1 (MQ=255)
                             gTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTgg                                                              >  2:162894/1‑100 (MQ=255)
                               tGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTgggg                                                            >  1:1440396/1‑100 (MQ=255)
                               tGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTgggg                                                            >  1:196766/1‑100 (MQ=255)
                                gcgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGc                                                           >  2:289877/1‑100 (MQ=255)
                                gcgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGc                                                           >  1:1658041/1‑100 (MQ=255)
                                gcgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGc                                                           >  1:1003087/1‑100 (MQ=255)
                                 cgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGcc                                                          <  1:1540103/100‑1 (MQ=255)
                                 cgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGcc                                                          >  2:589287/1‑100 (MQ=255)
                                 cgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGcc                                                          <  2:1119721/100‑1 (MQ=255)
                                   cATTCAGGATGCGACGATGATCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAg                                                        >  1:1188016/1‑100 (MQ=255)
                                   cATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAg                                                        <  2:777575/100‑1 (MQ=255)
                                    aTTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGc                                                       <  1:1010407/100‑1 (MQ=255)
                                    aTTCAGGATGCGACGATGACCCATGACCATCTTCCTTTCGTCGATGGTCAGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTGGGGGCCAg                                                        >  1:1879954/1‑99 (MQ=255)
                                     ttCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCa                                                      >  1:1072316/1‑100 (MQ=255)
                                        aGGATGCGACGATGATCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATc                                                   <  2:1914746/100‑1 (MQ=255)
                                        aGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATc                                                   <  2:642646/100‑1 (MQ=255)
                                         ggATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATcc                                                  <  1:715806/100‑1 (MQ=255)
                                         ggATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATc                                                   <  2:79271/99‑1 (MQ=255)
                                            tGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATc                                               <  2:301642/100‑1 (MQ=255)
                                            tGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATc                                               <  2:193436/100‑1 (MQ=255)
                                               gacgaTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGAc                                            >  1:596839/1‑100 (MQ=255)
                                                acgaTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACt                                           <  2:1930715/100‑1 (MQ=255)
                                                  gaTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTtg                                         <  2:596841/100‑1 (MQ=255)
                                                      aCCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGa                                     >  2:102594/1‑100 (MQ=255)
                                                      aCCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGa                                     <  2:695643/100‑1 (MQ=255)
                                                      aCCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGa                                     <  1:1803866/100‑1 (MQ=255)
                                                       cccATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGAt                                    >  2:1538372/1‑100 (MQ=255)
                                                         cATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCAATCGACTTGTGGATTg                                  >  2:1591648/1‑100 (MQ=255)
                                                            gACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCt                               <  1:1260399/100‑1 (MQ=255)
                                                              cGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTAc                             <  1:1867062/100‑1 (MQ=255)
                                                              cGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTAc                             <  1:1634176/100‑1 (MQ=255)
                                                                aTCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGa                           >  1:1530837/1‑100 (MQ=255)
                                                                 tCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGAc                          <  2:374375/100‑1 (MQ=255)
                                                                 tCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGAc                          >  2:1334509/1‑100 (MQ=255)
                                                                    gCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGTCCAGCTATCCATCGACTTGTGGATTGGCTACGACCGc                       >  1:1308668/1‑100 (MQ=255)
                                                                          cgtcgATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTa                 >  1:1265217/1‑100 (MQ=255)
                                                                           gtcgATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGCATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTac                >  2:285186/1‑100 (MQ=255)
                                                                            tcgATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTaca               <  2:1567163/100‑1 (MQ=255)
                                                                            tcgATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTaca               >  2:702117/1‑100 (MQ=255)
                                                                             cgATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTacac              <  2:143043/100‑1 (MQ=255)
                                                                                tCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACaa            >  1:1945474/1‑99 (MQ=255)
                                                                                  gTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAAtt         <  2:1031434/100‑1 (MQ=255)
                                                                                  gTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAAtt         >  1:1541525/1‑100 (MQ=255)
                                                                                     tGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTAt      <  1:1849267/100‑1 (MQ=255)
                                                                                     tGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTAt      >  2:1564076/1‑100 (MQ=255)
                                                                                      ggTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTAtt     <  1:188890/100‑1 (MQ=255)
                                                                                      ggTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTAtt     <  2:982116/100‑1 (MQ=255)
                                                                                         gtTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGt  >  2:1510137/1‑100 (MQ=255)
                                                                                               |                                                                                             
GCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGT  >  NC_000913/974988‑975176

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