Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	400,214	T→C	D237G (GAC→GGC)	ddlA ←	D‑alanine‑D‑alanine ligase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	400,214	0	T	C	93.8%	255.5 / 4.7	81	D237G (GAC→GGC)	ddlA	D‑alanine‑D‑alanine ligase A
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base C (32/44);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

TACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCG  >  NC_000913/400117‑400312                                                                                                  |                                                                                                   tACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCg                                                                                                  <  2:316696/100‑1 (MQ=255) tACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCg                                                                                                  <  2:1215518/100‑1 (MQ=255)  aCCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGt                                                                                                 >  2:408267/1‑100 (MQ=255)  aCCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGt                                                                                                 >  2:615906/1‑100 (MQ=255)  aCCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGt                                                                                                 >  1:900267/1‑100 (MQ=255)   ccACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGt                                                                                                 <  1:1227115/99‑1 (MQ=255)    cACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTg                                                                                               <  1:1089684/100‑1 (MQ=255)         tCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCGCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAg                                                                                          <  1:1047341/100‑1 (MQ=255)          cgcgCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAg                                                                                          <  2:139735/99‑1 (MQ=255)           gcgcCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGaa                                                                                        >  1:955155/1‑100 (MQ=255)            cgcCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAAc                                                                                       <  2:301594/100‑1 (MQ=255)              ccATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGGAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTg                                                                                     <  1:1449902/100‑1 (MQ=255)              ccATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTg                                                                                     <  1:1449787/100‑1 (MQ=255)                aTCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTgcg                                                                                   <  2:137416/100‑1 (MQ=255)                 tCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTgcgc                                                                                  <  1:458872/100‑1 (MQ=255)                    tCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCAtt                                                                               >  1:1117936/1‑100 (MQ=255)                     cGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTc                                                                              >  1:62134/1‑100 (MQ=255)                      gTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTc                                                                              >  2:588164/1‑99 (MQ=255)                       tCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGa                                                                            <  1:579347/100‑1 (MQ=255)                        cAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGAt                                                                           <  2:866776/100‑1 (MQ=255)                        cAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGAt                                                                           >  1:9184/1‑100 (MQ=255)                             tACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCAc                                                                      >  1:1796673/1‑100 (MQ=255)                             tACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCAc                                                                      >  1:1796691/1‑100 (MQ=255)                              aCTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACg                                                                     >  2:1201650/1‑100 (MQ=255)                                     gtCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTg                                                              <  1:1404440/100‑1 (MQ=255)                                     gtCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTg                                                              >  2:890420/1‑100 (MQ=255)                                        gTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATc                                                           >  1:1771632/1‑100 (MQ=255)                                          aGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGCACTGCGCATTCGATCTCACGACCTTTGATccc                                                         >  1:398320/1‑100 (MQ=255)                                           ggCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCt                                                        >  2:512236/1‑100 (MQ=255)                                           ggCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCt                                                        >  1:1834273/1‑100 (MQ=255)                                            gCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCtt                                                       >  2:105838/1‑100 (MQ=255)                                             cATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTg                                                      <  2:900268/100‑1 (MQ=255)                                              aTAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGc                                                     <  2:469641/100‑1 (MQ=255)                                                 gAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCa                                                  >  1:1248630/1‑100 (MQ=255)                                                    aTCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACg                                               <  2:77229/100‑1 (MQ=255)                                                     tCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGa                                              <  2:1228065/100‑1 (MQ=255)                                                      cGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGan                                             <  2:615236/100‑2 (MQ=255)                                                      cGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGAt                                             >  2:106219/1‑100 (MQ=255)                                                       gatcttGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATc                                            <  2:1595153/95‑1 (MQ=255)                                                         tGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCAc                                          <  2:934843/100‑1 (MQ=255)                                                          ggTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACt                                         >  2:907491/1‑100 (MQ=255)                                                           gTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACtt                                        >  2:583128/1‑100 (MQ=255)                                                            tGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCagatc                                       >  1:1595142/1‑96 (MQ=255)                                                            tGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACttt                                       <  1:980687/100‑1 (MQ=255)                                                             gAGTACGATCTCTCCACAGGTGCTGTCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTa                                      <  1:320163/100‑1 (MQ=255)                                                              aGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTAt                                     >  1:1535715/1‑100 (MQ=255)                                                               gTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATg                                    <  1:1600976/100‑1 (MQ=255)                                                                tACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGa                                   <  1:623447/100‑1 (MQ=255)                                                                 aCGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAt                                  <  1:408265/100‑1 (MQ=255)                                                                 aCGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAt                                  <  1:615905/100‑1 (MQ=255)                                                                  cGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtc                                 <  1:945689/100‑1 (MQ=255)                                                                  cGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtc                                 >  1:1910369/1‑100 (MQ=255)                                                                  cGATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAt                                  >  1:1378443/1‑99 (MQ=255)                                                                   gATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGTACTGCGCATTCGATCTCACGACATTTGATCCCTTGCTCAACGATCACTTTATGAtcg                                <  1:1362980/100‑1 (MQ=255)                                                                   gATCTCTCCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGGACTGCGCATTCGATCTCGCGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcg                                <  1:1363456/100‑1 (MQ=255)                                                                      ctctcCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaac                             >  2:1069763/1‑100 (MQ=255)                                                                         tcCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcg                          <  1:1266053/100‑1 (MQ=255)                                                                         tcCACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcg                          >  2:473074/1‑100 (MQ=255)                                                                          ccACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcga                         <  1:687844/100‑1 (MQ=255)                                                                          ccACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcga                         >  1:1646082/1‑100 (MQ=255)                                                                          ccACAGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTCATCCCTTGCTCAACGATCACTTTATGAtcgaactcga                         <  1:1041159/100‑1 (MQ=255)                                                                           cacaGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcg                          <  2:1010420/98‑1 (MQ=255)                                                                            acaGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaactcgaac                       <  2:76883/100‑1 (MQ=255)                                                                             caGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACg                      <  1:972872/100‑1 (MQ=255)                                                                             caGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACg                      <  2:1196884/100‑1 (MQ=255)                                                                             caGGTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACg                      >  2:1615543/1‑100 (MQ=255)                                                                               ggTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGcc                    <  2:1917378/100‑1 (MQ=255)                                                                                gTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCa                   <  2:1796684/100‑1 (MQ=255)                                                                                gTGCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCa                   <  2:1796702/100‑1 (MQ=255)                                                                                  gCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCGTGCTCACCGATCACTTTATGATCGAACTCGAACGCCAGa                 >  1:1366730/1‑100 (MQ=255)                                                                                  gCTGGCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGa                 <  2:1725310/100‑1 (MQ=255)                                                                                   cTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGAt                >  2:1812813/1‑100 (MQ=255)                                                                                      gCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGa             <  2:1528832/100‑1 (MQ=255)                                                                                      gCTTGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGa             >  2:934637/1‑100 (MQ=255)                                                                                         tGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACgg          <  2:1458626/100‑1 (MQ=255)                                                                                         tGCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACgg          <  1:192734/100‑1 (MQ=255)                                                                                          gCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGc         >  1:1473301/1‑100 (MQ=255)                                                                                          gCGGATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGc         <  2:1117937/100‑1 (MQ=255)                                                                                             gATTGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAAt      >  2:804229/1‑100 (MQ=255)                                                                                               ttGCCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTg    >  2:1802576/1‑100 (MQ=255)                                                                                                 gTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCg  >  2:737216/1‑100 (MQ=255)                                                                                                  |                                                                                                   TACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCG  >  NC_000913/400117‑400312

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