Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,335,418	A→G	S522G (AGC→GGC)	acnA →	aconitate hydrase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,335,418	0	A	G	98.9%	228.5 / NA	88	S522G (AGC→GGC)	acnA	aconitate hydrase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (35/52);  minor base T (0/1);  total (35/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAA  >  NC_000913/1335348‑1335489                                                                       |                                                                        gtACGGTTGTACCTCCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAACCGGCAATCAAAAAAGGCGAt                                                                    <  2:4025625/75‑1 (MQ=255) aTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAt                                                                    <  2:2520003/76‑1 (MQ=255) aTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAt                                                                    <  2:2374224/76‑1 (MQ=255) aTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAt                                                                    <  1:2723579/76‑1 (MQ=255)    cGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTa                                                                 >  2:3906238/1‑76 (MQ=255)     ggTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTaa                                                                <  2:3347789/76‑1 (MQ=255)      gTTGTACCGCCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAAc                                                               <  2:792926/76‑1 (MQ=255)      gTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAAc                                                               <  2:3705417/76‑1 (MQ=255)       ttGTGCCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAAcc                                                              <  1:2515218/76‑1 (MQ=255)       ttGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAAcc                                                              <  2:4676413/76‑1 (MQ=255)       ttGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAAcc                                                              <  1:1081677/76‑1 (MQ=255)        tGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGCTTTAAccc                                                             >  2:3616561/1‑75 (MQ=255)        tGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCg                                                             <  1:2754223/76‑1 (MQ=255)        tGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCg                                                             <  1:1877814/76‑1 (MQ=255)        tGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCg                                                             <  2:4570899/76‑1 (MQ=255)        tGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCg                                                             <  1:2171522/76‑1 (MQ=255)        tGTACCACCTGTATTGGTAACTCTGGCCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCg                                                             >  2:2305788/1‑76 (MQ=255)          tACCGCCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                           <  1:3675/76‑1 (MQ=255)          tACCGCCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                           <  2:2379795/76‑1 (MQ=255)          tACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                           <  2:4579426/76‑1 (MQ=255)          tACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                           >  1:3958721/1‑76 (MQ=255)          tACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                           <  1:3403747/76‑1 (MQ=255)            ccaccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcgg                                                         >  1:1734055/1‑76 (MQ=255)            ccaccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcgg                                                         >  2:1234811/1‑76 (MQ=255)             caccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggt                                                        <  1:152012/76‑1 (MQ=255)              cccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtg                                                       <  1:3334386/75‑1 (MQ=255)               ccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgc                                                      <  1:3747251/76‑1 (MQ=255)               ccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgc                                                      <  1:2230661/76‑1 (MQ=255)               ccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgc                                                      <  1:1851759/76‑1 (MQ=255)                cTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcg                                                     <  1:4474857/76‑1 (MQ=255)                cTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcg                                                     <  1:87352/76‑1 (MQ=255)                   tATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcggtg                                                  >  1:2395744/1‑76 (MQ=255)                    aTTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCggtgcggtgc                                                 >  2:1623004/1‑76 (MQ=255)                    aTTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGGCggtgcggtgc                                                 <  1:3940536/76‑1 (MQ=255)                      tGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTg                                               <  2:2100356/76‑1 (MQ=255)                       ggTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTgc                                              >  2:3361987/1‑75 (MQ=255)                       ggTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGt                                              <  1:1674517/76‑1 (MQ=255)                        gTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTc                                             >  2:879825/1‑76 (MQ=255)                        gTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTc                                             >  1:4492710/1‑76 (MQ=255)                          aaCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCg                                           >  2:3716423/1‑76 (MQ=255)                           aCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCgg                                          <  2:1616980/76‑1 (MQ=255)                             tctGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCa                                        >  1:3711986/1‑76 (MQ=255)                             tcggggCCGCTGCCCGATCCTATAGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGAGCGGCGCTGTCCAGCa                                        >  1:2064349/4‑76 (MQ=255)                              ctGGCCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCaa                                       <  2:4361919/76‑1 (MQ=255)                               tGGGCCGCTGCGCGATCCTATCGAAACGGCAATGAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAc                                      <  1:3551840/76‑1 (MQ=255)                               tGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAc                                      >  1:4652405/1‑76 (MQ=255)                               tGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAc                                      <  1:3052867/76‑1 (MQ=255)                                gggCGGCTGCCCGATCCTATGGAGCGGGCAATAAAAAACTGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAcc                                     <  2:2064348/76‑1 (MQ=255)                                 ggCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCg                                    >  1:3453178/1‑76 (MQ=255)                                  gCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGt                                   <  1:1341301/76‑1 (MQ=255)                                   ccGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTa                                  >  1:2605573/1‑76 (MQ=255)                                    cGCTGCGCGCTCCTTTCGGAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGGGGTGCTGTCCGGCAACCGTaa                                 <  2:1573449/76‑1 (MQ=255)                                    cGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTaa                                 >  1:4473220/1‑76 (MQ=255)                                      cggCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACt                               <  2:1185507/74‑1 (MQ=255)                                      cTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACt                               >  1:86489/1‑76 (MQ=255)                                      cTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACt                               >  2:1087400/1‑76 (MQ=255)                                      cTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGGAACt                               >  2:1645018/1‑76 (MQ=255)                                       tGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGCCGGTGCGGGGCTGTCCGGCAACCGCAACtt                              >  2:4630919/1‑76 (MQ=255)                                        gCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACttt                             <  2:1055847/76‑1 (MQ=255)                                         cccGATCCTAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTg                            <  1:3090975/76‑1 (MQ=255)                                          ccGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGa                           <  2:2601824/76‑1 (MQ=255)                                           cGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGGTGTCCGGCAACCGTAACTTTGaa                          <  2:1669024/76‑1 (MQ=255)                                           cGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGaa                          >  1:2352470/1‑76 (MQ=255)                                              tacTACCGAAACGGCAATCAAAAAAGCCGATTTAACCGACGCTGCGGTGCCGCCCCACAAGCGTAACTTGCAAGGc                       >  2:2643528/3‑76 (MQ=255)                                              tCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGc                       <  2:3281400/76‑1 (MQ=255)                                               ccTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGcc                      >  2:3675018/1‑76 (MQ=255)                                               ccTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGcc                      <  1:4399368/76‑1 (MQ=255)                                               catATCGAAACGGCAATAAAAAAAGGCGATTTAACCGTAGGTGCGGTGCTGTCCGGCAACCGTAAATCTGAAGGcc                      >  1:3195393/3‑76 (MQ=255)                                                cTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCg                     >  1:1213710/1‑76 (MQ=255)                                                cTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCg                     <  1:717845/76‑1 (MQ=255)                                                 tATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGt                    <  2:4652405/76‑1 (MQ=255)                                                   tCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTat                  >  1:4611185/1‑76 (MQ=255)                                                   tCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGGAACTTTGAAGGCCGTat                  >  1:2013082/1‑76 (MQ=255)                                                   tCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGCCAACCCTCACTTTGAAGGCCgcac                  >  2:279823/1‑73 (MQ=255)                                                    cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                 <  1:2094956/76‑1 (MQ=255)                                                     gAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGCatcc                >  2:874692/1‑76 (MQ=255)                                                      taaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatcca               >  2:3327074/2‑76 (MQ=255)                                                       aaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccat              >  2:2392344/1‑76 (MQ=255)                                                         cGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccatcc            <  2:835939/76‑1 (MQ=255)                                                         cGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccatcc            <  2:1030117/76‑1 (MQ=255)                                                          ggCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCg           <  1:2606577/76‑1 (MQ=255)                                                          ggCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCg           <  1:1400503/76‑1 (MQ=255)                                                             aaTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTg        >  1:3304797/1‑76 (MQ=255)                                                             aaTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTg        >  2:4687966/1‑76 (MQ=255)                                                              aTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTgg       <  2:1352072/76‑1 (MQ=255)                                                              aTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTgg       <  1:4568767/76‑1 (MQ=255)                                                              aTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTgg       <  2:3498880/76‑1 (MQ=255)                                                                cgaaaaagGCGACTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGtt     >  1:2261780/3‑76 (MQ=38)                                                                cAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGtt     <  1:323705/76‑1 (MQ=255)                                                                 aaaaaagGCGATTTACCCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTa    <  1:2156048/76‑1 (MQ=255)                                                                   aaaagGCGATTTAACCGTCGGTGCGGTGCGGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTaaa  <  2:4041121/76‑1 (MQ=255)                                                                       |                                                                        ATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAA  >  NC_000913/1335348‑1335489

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