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_000913	975,083	0	T	A	98.6%	236.8 / ‑3.7	69	D255E (GAT→GAA)	mukF	chromosome condensin MukBEF, kleisin‑like subunit, binds calcium
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (29/39);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

GCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGTACC  >  NC_000913/974988‑975179                                                                                                |                                                                                                 gCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGc                                                                                              <  1:762868/100‑1 (MQ=255) gCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGc                                                                                              <  1:5662/100‑1 (MQ=255)   ggCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCag                                                                                            <  2:919438/100‑1 (MQ=255)   ggCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCag                                                                                            <  1:1238116/100‑1 (MQ=255)    gCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCa                                                                                             <  2:365440/98‑1 (MQ=255)      aGGCGGCGAATAGCAGGCTAATCTGTGGCGCATTCGGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCagag                                                                                          <  2:245119/99‑1 (MQ=255)         cGACAAATTGTAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCaaa                                                                                      <  1:493119/100‑1 (MQ=255)          gACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGCACTGCAGAGCaa                                                                                       >  2:75221/1‑98 (MQ=255)                ttGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATc                                                                               <  2:785461/100‑1 (MQ=255)                 tGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCg                                                                              <  1:564149/100‑1 (MQ=255)                   cAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGta                                                                            <  2:771697/100‑1 (MQ=255)                     ggCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGtatt                                                                          >  1:596233/1‑100 (MQ=255)                     ggCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGtatt                                                                          <  2:395503/100‑1 (MQ=255)                        tAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATc                                                                       <  1:727522/100‑1 (MQ=255)                              gTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTgg                                                                 >  2:1282127/1‑100 (MQ=255)                                tGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTgggg                                                               <  2:1193053/100‑1 (MQ=255)                                  cgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGcc                                                             <  2:95034/100‑1 (MQ=255)                                  cgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGcc                                                             <  2:313404/100‑1 (MQ=255)                                   gcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCa                                                            >  2:890972/1‑100 (MQ=255)                                    cATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAg                                                           <  1:278446/100‑1 (MQ=255)                                    cATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAg                                                           <  1:253380/100‑1 (MQ=255)                                    cATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAg                                                           <  2:1422163/100‑1 (MQ=255)                                     aTTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGc                                                          >  1:483273/1‑100 (MQ=255)                                     aTTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGc                                                          >  2:853232/1‑100 (MQ=255)                                     aTTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGc                                                          >  2:901425/1‑100 (MQ=255)                                     aTTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGc                                                          <  1:1274394/100‑1 (MQ=255)                                        cAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAAt                                                       <  2:852603/100‑1 (MQ=255)                                        cAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAAt                                                       >  1:616804/1‑100 (MQ=255)                                        cAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAAt                                                       >  1:1289475/1‑100 (MQ=255)                                         aGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATc                                                      <  1:734889/100‑1 (MQ=255)                                          ggATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATcc                                                     >  1:362083/1‑100 (MQ=255)                                          ggATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATcc                                                     >  1:1390593/1‑100 (MQ=255)                                             tGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATc                                                  <  1:759804/100‑1 (MQ=255)                                               cgacgaTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGa                                                >  1:1155643/1‑100 (MQ=255)                                               cgacgaTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGa                                                >  1:913931/1‑100 (MQ=255)                                               cgacgaTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGa                                                >  1:878183/1‑100 (MQ=255)                                                     tGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTtgtg                                          <  2:726373/100‑1 (MQ=255)                                                       aCCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGa                                        <  2:155582/100‑1 (MQ=255)                                                        cccATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGAt                                       >  1:959641/1‑100 (MQ=255)                                                        cccATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGAt                                       >  1:424257/1‑100 (MQ=255)                                                        cccATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGAt                                       >  1:51346/1‑100 (MQ=255)                                                         ccATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGAtt                                      <  2:491810/100‑1 (MQ=255)                                                           aTGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTgg                                    <  1:133979/100‑1 (MQ=255)                                                           aTGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTgg                                    >  1:1406497/1‑100 (MQ=255)                                                           aTGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTgg                                    <  2:838944/100‑1 (MQ=255)                                                            tGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGc                                   <  2:1034940/100‑1 (MQ=255)                                                              aCGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCtn                                 <  2:788471/100‑2 (MQ=255)                                                                gATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACg                               <  2:200679/100‑1 (MQ=255)                                                                  tCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGAc                             >  2:888192/1‑100 (MQ=255)                                                                   cTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGAcc                            >  1:698557/1‑100 (MQ=255)                                                                    tGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGAccn                           <  1:1141078/100‑2 (MQ=255)                                                                     gCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGc                          >  2:1101946/1‑100 (MQ=255)                                                                       aTTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCa                        <  1:597891/100‑1 (MQ=255)                                                                       aTTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCa                        >  1:299462/1‑100 (MQ=255)                                                                        tttCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCAc                       >  1:663646/1‑100 (MQ=255)                                                                           cgtcgATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTa                    >  1:494209/1‑100 (MQ=255)                                                                            gtcgATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTac                   >  1:1581501/1‑100 (MQ=255)                                                                            gtcgATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGACCGTATTATCAGTTGGGGCCAGCAATCCATCGCCTTGTGGATTAGCTACGACCGCCACGTac                   >  2:652919/1‑100 (MQ=255)                                                                                aTCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACaa               <  2:635309/100‑1 (MQ=255)                                                                                  cGTCTGGTGTTCGAACTGCAGAACAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAAt             >  2:918950/1‑100 (MQ=255)                                                                                   gTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAAtt            <  2:1406489/100‑1 (MQ=255)                                                                                       ggTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTAtt        <  2:6118/100‑1 (MQ=255)                                                                                        gtgtTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTAtt        <  2:786012/99‑1 (MQ=255)                                                                                        gtgtTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTAtt        <  2:1155638/99‑1 (MQ=255)                                                                                         tgtTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCg      <  1:1196198/100‑1 (MQ=255)                                                                                          gtTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGt     >  1:367204/1‑100 (MQ=255)                                                                                          gtTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGt     >  1:209398/1‑100 (MQ=255)                                                                                            tCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGTAc   <  2:1418069/100‑1 (MQ=255)                                                                                             cGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGTAcc  >  2:440624/1‑100 (MQ=255)                                                                                                |                                                                                                 GCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGTACC  >  NC_000913/974988‑975179

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