Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	975,083	T→A	15.7%	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	15.7%	179.3 / 24.6	70	D255E (GAT→GAA)	mukF	chromosome condensin MukBEF, kleisin‑like subunit, binds calcium
Reads supporting (aligned to +/- strand):  ref base T (36/23);  new base A (5/6);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCG  >  NC_000913/974984‑975175                                                                                                    |                                                                                             ggAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAt                                                                                              <  1:2257419/100‑1 (MQ=255)  gAAGCGGCAGGCGACAAATTGCAGTCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATc                                                                                             <  2:1024900/100‑1 (MQ=255)  gAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATc                                                                                             <  1:421070/100‑1 (MQ=255)   aaGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCt                                                                                            <  2:460516/100‑1 (MQ=255)       ggCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCag                                                                                        <  2:355612/100‑1 (MQ=255)        gCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCaga                                                                                       >  2:1947512/1‑100 (MQ=255)            gCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCaa                                                                                   >  2:1298138/1‑100 (MQ=255)            gCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCa                                                                                    >  2:1296809/1‑99 (MQ=255)                cAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTc                                                                               >  2:1711646/1‑100 (MQ=255)                cAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTc                                                                               >  1:327399/1‑100 (MQ=255)                  aaTTGCAGGCTAATCTGTTGCGGATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGa                                                                             >  2:1391931/1‑100 (MQ=255)                  aaTTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGa                                                                             >  2:1520484/1‑100 (MQ=255)                  aaTTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGa                                                                             >  1:1431731/1‑100 (MQ=255)                  aaTTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGa                                                                             >  2:1890044/1‑100 (MQ=255)                     tGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCg                                                                          <  2:545852/100‑1 (MQ=255)                       cAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGta                                                                        >  1:518800/1‑100 (MQ=255)                        aGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGtat                                                                       >  2:2094170/1‑100 (MQ=255)                         ggCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGtatt                                                                      >  1:2064434/1‑100 (MQ=255)                         ggCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGtatt                                                                      >  1:465603/1‑100 (MQ=255)                          gCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGtatta                                                                     <  1:242441/100‑1 (MQ=255)                           cTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGtattat                                                                    <  1:889350/100‑1 (MQ=255)                              aTCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAg                                                                 <  2:2049300/100‑1 (MQ=255)                                 tGTTGCGCATTCAGGATGCGACGGTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTg                                                              <  2:1208467/100‑1 (MQ=255)                                 tGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTg                                                              <  1:975393/100‑1 (MQ=255)                                    tGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTgggg                                                           <  2:2348505/100‑1 (MQ=255)                                      cgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGcc                                                         >  2:1700978/1‑100 (MQ=255)                                      cgcATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGc                                                          <  2:2342887/99‑1 (MQ=255)                                        cATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAg                                                       >  1:2148075/1‑100 (MQ=255)                                        cATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAg                                                       >  1:1674680/1‑100 (MQ=255)                                         aTTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGc                                                      <  2:1557110/100‑1 (MQ=255)                                          ttCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCa                                                     <  1:1338467/100‑1 (MQ=255)                                           tCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTACCAGTTGGGGCCAGCaa                                                    >  2:1397366/1‑100 (MQ=255)                                            cAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAAt                                                   >  2:247104/1‑100 (MQ=255)                                            cAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAAt                                                   <  1:1128276/100‑1 (MQ=255)                                            cAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAAt                                                   >  1:1082121/1‑100 (MQ=255)                                             aGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATc                                                  >  1:53535/1‑100 (MQ=255)                                             aGGATGCGACGATGACCCATGACGAGCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATc                                                  <  2:1634956/100‑1 (MQ=255)                                                 tGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCAtn                                              <  2:1104272/100‑2 (MQ=255)                                                  gCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCg                                             >  2:27217/1‑100 (MQ=255)                                                    gacgaTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATGCATCGAc                                           <  1:1255234/100‑1 (MQ=255)                                                    gacgaTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGAc                                           >  1:2250295/1‑100 (MQ=255)                                                     acgaTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACt                                          >  2:2119980/1‑100 (MQ=255)                                                        aTGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTtgt                                       <  2:308119/100‑1 (MQ=255)                                                         tGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTtgtg                                      >  2:1754687/1‑100 (MQ=255)                                                         tGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTtgtg                                      <  1:247107/100‑1 (MQ=255)                                                           aCCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGa                                    <  2:782839/100‑1 (MQ=255)                                                            cccATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGAt                                   >  2:509427/1‑100 (MQ=255)                                                             ccATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGAtt                                  >  2:6444/1‑100 (MQ=255)                                                              cATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCAAGCAATCCATCGACTTGTGGATTg                                 >  2:2028047/1‑100 (MQ=255)                                                              cATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTg                                 <  2:2228307/100‑1 (MQ=255)                                                                tGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGc                               >  2:1306474/1‑100 (MQ=255)                                                                 gACGATCTGCATTTCGTCGATCGTCTGGTGTTCGCACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTAGTGCATTGGCt                              <  2:1259926/100‑1 (MQ=255)                                                                     aTCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGa                          >  1:1928322/1‑100 (MQ=255)                                                                       cTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGAcc                        >  1:1478476/1‑100 (MQ=255)                                                                       cTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGATGGGTCCAGCAATCCATCGACTTGTGGATTGGCTACGAcc                        >  2:723014/1‑100 (MQ=255)                                                                         gCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGc                      >  2:2301461/1‑100 (MQ=255)                                                                         gCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGc                      >  1:1865398/1‑100 (MQ=255)                                                                         gCATTTCGTCGATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGc                      >  1:489829/1‑100 (MQ=255)                                                                          cATTTCGTCGATCGTTTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATAAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGcc                     <  1:642053/100‑1 (MQ=255)                                                                               cgtcgATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTa                >  2:2015214/1‑100 (MQ=255)                                                                               cgtcgATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTa                >  1:1316859/1‑100 (MQ=255)                                                                               cgtcgATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTa                >  1:2389627/1‑100 (MQ=255)                                                                               cgtcgATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTa                >  1:617163/1‑100 (MQ=255)                                                                               cgtcgATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTa                >  2:1032397/1‑100 (MQ=255)                                                                                gtcgATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTac               >  1:1249920/1‑100 (MQ=255)                                                                                gtcgATCGTCTGGTGTTCGAACTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTac               >  1:1004557/1‑100 (MQ=255)                                                                                      cGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAAt         <  1:1290814/100‑1 (MQ=255)                                                                                          tGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTAt     <  1:1391945/100‑1 (MQ=255)                                                                                            gtgtTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTc   <  1:2128708/100‑1 (MQ=255)                                                                                            gtgtTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTc   <  1:940495/100‑1 (MQ=255)                                                                                             tgtTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCg  <  2:503833/100‑1 (MQ=255)                                                                                                    |                                                                                             GGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCG  >  NC_000913/974984‑975175

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