Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	234,095	C→T	E230K (GAA→AAA)	gloB ←	hydroxyacylglutathione hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	234,095	0	C	T	98.0%	174.5 / ‑4.2	51	E230K (GAA→AAA)	gloB	hydroxyacylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/20);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTT  >  NC_000913/233999‑234190                                                                                                 |                                                                                                gAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtatt                                                                                              >  1:548625/1‑100 (MQ=255) gAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtatt                                                                                              <  1:546253/100‑1 (MQ=255) gAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtatt                                                                                              >  2:1485481/1‑100 (MQ=255) gAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTGCTTcatt                                                                                              <  1:827870/100‑1 (MQ=255)  aaCGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTttatta                                                                                             <  1:755820/100‑1 (MQ=255)  aaCGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTttatta                                                                                             <  1:755847/100‑1 (MQ=255)     ggAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtattaatt                                                                                          <  2:1241952/100‑1 (MQ=255)       aGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtattaattac                                                                                        <  2:1176520/100‑1 (MQ=255)        tctgttggcAAGAAAGTTTTATCAGAATCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtattaatt                                                                                          <  1:1414393/89‑1 (MQ=255)        gCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtattaattaca                                                                                       >  2:1304768/1‑100 (MQ=255)        gCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtattaatta                                                                                         <  2:1224321/98‑1 (MQ=255)          cGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtattaattacatt                                                                                     <  2:702800/100‑1 (MQ=255)          cGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtattaattacatt                                                                                     >  1:1437414/1‑100 (MQ=255)           gATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGTTCTTCAGGTTGTTGCAATAATGTTTCTTtattaattacatta                                                                                    <  2:684929/100‑1 (MQ=255)              gACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTtattaattacattaatt                                                                                 <  2:128279/100‑1 (MQ=255)               aCAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCCGGTTGTTGCAATAATGTTTCTTtattacttacatttatta                                                                                >  2:449843/1‑100 (MQ=255)                cAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACAttaattaa                                                                               >  2:725895/1‑100 (MQ=255)                 aaGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTaaa                                                                              >  1:1416904/1‑100 (MQ=255)                  aGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAAt                                                                             >  2:158653/1‑100 (MQ=255)                     aaGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCaa                                                                          >  1:1226848/1‑100 (MQ=255)                      aGTTTTATCAGAAACTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTAAAGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCaa                                                                          <  2:1077778/99‑1 (MQ=255)                          ttATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCt                                                                     >  1:1485596/1‑100 (MQ=255)                            aTCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGTTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCtt                                                                    >  2:587680/1‑99 (MQ=255)                            aTCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTc                                                                   >  1:763261/1‑100 (MQ=255)                             tCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTcc                                                                  >  1:480459/1‑100 (MQ=255)                              cAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTcc                                                                  >  2:1544132/1‑99 (MQ=255)                              cAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCg                                                                 >  2:1063907/1‑100 (MQ=255)                               aGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGt                                                                >  1:60640/1‑100 (MQ=255)                                    cTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTTTTa                                                           >  2:175319/1‑100 (MQ=255)                                      aTCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTaaa                                                         >  1:1528554/1‑100 (MQ=255)                                      aTCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTaaa                                                         >  1:1029689/1‑100 (MQ=255)                                         tttctttGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTaaaaaa                                                      <  1:1165701/100‑1 (MQ=255)                                           tctttGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAAc                                                    >  1:353092/1‑100 (MQ=255)                                           tctttGACCTTAACAATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAAc                                                    <  2:344314/100‑1 (MQ=255)                                            ctttGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACa                                                   >  2:209292/1‑100 (MQ=255)                                            ctttGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACa                                                   <  2:1520892/100‑1 (MQ=255)                                             tttGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACAt                                                  >  1:1243172/1‑100 (MQ=255)                                              ttGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACAtt                                                 >  2:952879/1‑100 (MQ=255)                                                gACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTaa                                               >  2:1380774/1‑100 (MQ=255)                                                  ccTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAAtt                                             <  1:280721/100‑1 (MQ=255)                                                     tAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGc                                          <  1:1507060/100‑1 (MQ=255)                                                     tAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGc                                          >  1:1281066/1‑100 (MQ=255)                                                     tAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGc                                          <  1:273799/100‑1 (MQ=255)                                                      aaCCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGcc                                         >  2:219615/1‑100 (MQ=255)                                                      aaCCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGcc                                         <  1:1304774/100‑1 (MQ=255)                                                       aCCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCg                                        <  1:1544594/100‑1 (MQ=255)                                                           tGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCa                                    >  1:980933/1‑100 (MQ=255)                                                             cAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCAtt                                  >  1:611237/1‑100 (MQ=255)                                                               aaaCGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCAtttt                                <  2:480458/100‑1 (MQ=255)                                                                 aCGCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTc                              >  1:529882/1‑100 (MQ=255)                                                                   gCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAg                            <  2:871582/100‑1 (MQ=255)                                                                   gCTCTTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAg                            <  1:1228861/100‑1 (MQ=255)                                                                      ctTCAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAACAACACATTAATTTGCCGCTCATTTTTCa                             >  2:947773/1‑96 (MQ=255)                                                                         cAGGTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTAc                      >  2:498614/1‑100 (MQ=255)                                                                           ggTTGTTGCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACgg                    >  1:549082/1‑100 (MQ=255)                                                                             ttgttgCAATAATTTTTCTATTTTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAGCATTAATTTGCCGCTCATTTTTCAGAATTACg                     <  2:1360301/97‑1 (MQ=255)                                                                                ttgCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAgt               <  1:722431/100‑1 (MQ=255)                                                                                  gCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAgtgt             <  1:734349/100‑1 (MQ=255)                                                                                  gCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAgtgt             >  1:405138/1‑100 (MQ=255)                                                                                  gCAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAgtgt             <  1:378750/100‑1 (MQ=255)                                                                                   cAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGtt            >  1:1398578/1‑100 (MQ=255)                                                                                   cAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGtt            <  2:605630/100‑1 (MQ=255)                                                                                   cAATAATGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGtt            >  1:1539696/1‑100 (MQ=255)                                                                                    aataatGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTa           <  2:505340/100‑1 (MQ=255)                                                                                       aatGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTAttt        >  1:502270/1‑100 (MQ=255)                                                                                       aatGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTAttt        <  2:827039/100‑1 (MQ=255)                                                                                       aatGTTTCTTTATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTAttt        <  2:1281061/100‑1 (MQ=255)                                                                                           ttgctttATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGAtt    <  2:752003/97‑1 (MQ=255)                                                                                             tctttATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGAtttt  >  1:119449/1‑100 (MQ=255)                                                                                                 |                                                                                                GAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTT  >  NC_000913/233999‑234190

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: ‑