Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,009,698	G→A	A177T (GCA→ACA)	pldB →	lysophospholipase L2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,009,698	0	G	A	96.8%	95.2 / NA	31	A177T (GCA→ACA)	pldB	lysophospholipase L2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/14);  minor base C (0/1);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCG  >  NC_000913/4009455‑4009934                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              ggTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCa                                                                                                                                                                                                                                        <  1:146956/251‑1 (MQ=255)                gATCCCCAGCTCGTGCATGTTAATCGTTGTAATGATTATGTTAATGATCTGGCGGCATTCTGGCAGCAGTAGGTTCAGCCCGGTCCGTGGCTTTAACGCTATATACTGGCACTGTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACTCCATCCAGGTGTATGTGACGCCATTGCGCTATCTTCGCCAATGTTTGGGAGCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGc                                                                                                                                                                                                                         <  2:121578/251‑1 (MQ=255)                                aTGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGt                                                                                                                                                                                                         >  1:93724/1‑251 (MQ=255)                                aTGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGAGCGATCTCCACATTATTTCTGCAACGCCATCCCGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGt                                                                                                                                                                                                         >  2:106259/1‑251 (MQ=255)                                        cGCTTTAATGATGAGGTTGAGGATCTGGCGGTCTTCTGGCAGCAGGACGATCAGCCCGGTCGGTGGCGTAAACGCTATATACTGGCTCATTCGAGGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTTCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTACCGTg                                                                                                                                                                                                  <  2:137090/250‑1 (MQ=255)                                                                cTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGACACCGGgcg                                                                                                                                                                        >  1:247336/1‑251 (MQ=255)                                                                  ggcgtcATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCgg                                                                                                                                                                           <  1:238741/246‑1 (MQ=255)                                                                     ggcATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGgc                                                                                                                                                                   <  1:219777/251‑1 (MQ=255)                                                                                        gTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTAt                                                                                                                                                >  1:122659/1‑251 (MQ=255)                                                                                          tCAGCCCGGTCCGTGTCGTAAACGCTATATTCTGGCACATTCGATGGGCGGTGCGATCTCCACATTAGTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATc                                                                                                                                               <  1:151196/250‑1 (MQ=255)                                                                                                     cGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGAcc                                                                                                                                   <  2:103315/251‑1 (MQ=255)                                                                                                      gTGGCGTAAACGCTATATACGGGCACATGCGATGGGCGGTGCGATCTCCCCATTATTTCTGCAACGCCATCCAGGTTTATGTGACGCCATTGCGCTAACTGCGCCTATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTGTGCTATCAACGTACTGAccc                                                                                                                                  <  1:299491/251‑1 (MQ=255)                                                                                                                        aCTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATc                                                                                                                <  1:186759/251‑1 (MQ=255)                                                                                                                             cacaTTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGt                                                                                                           >  1:32789/1‑251 (MQ=255)                                                                                                                                aTTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAAc                                                                                                        >  1:270569/1‑251 (MQ=255)                                                                                                                                               cGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCt                                                                                         >  1:213007/1‑251 (MQ=255)                                                                                                                                                   ctcCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtgatg                                                                                     >  1:259963/1‑251 (MQ=255)                                                                                                                                                    cccACATTATTTCTGCAACGCCATCCAGGTGTATGTGTCGCCATTGCGCTAACTGCGCCAATGGTTGGGATCGTGATTCGTATGCCGGCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtgatga                                                                                    <  2:270569/250‑1 (MQ=255)                                                                                                                                                     ccACAGTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGAACGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtgatg                                                                                     <  2:85424/249‑1 (MQ=255)                                                                                                                                                            attgcTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACg                                                                             <  1:319801/246‑1 (MQ=255)                                                                                                                                                                          aTCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTggg                                                              >  1:156909/1‑251 (MQ=255)                                                                                                                                                                                 tgtATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCt                                                       >  2:243783/1‑251 (MQ=255)                                                                                                                                                                                   tATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTAc                                                     >  1:221609/1‑251 (MQ=255)                                                                                                                                                                                       tgACGCCATTGCGCTAACTGCGCCAATGTTTGGGCTCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGAACAGACCCACAGCAGGAAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCACCGATTCGCGTCGGTGGGCCGACCTACCAtt                                                 >  2:323753/1‑251 (MQ=255)                                                                                                                                                                                                         tGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTAATTTATGACACGGCAGCTCCTCAACTGGGCCGAAGCACACCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGATCTTGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCAAAGAAGACAGCGATATCGACGCAACTTACGCTTCTCTGCTGATGACCCAACGATTCGCATCGGTGGGCCGACCTACCATTGGGTACGCGAACGTAttt                               >  2:281453/1‑250 (MQ=255)                                                                                                                                                                                                                            gTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCCACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCCCAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGACAGTATTCTGGCTGGCGAACAGGTGtt            >  1:88770/1‑251 (MQ=255)                                                                                                                                                                                                                             ggatgcgtttGCCGTCATTTAGTACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACTTTGCCGTGTG‑GGTTAGGCAATAGGCACCGGTCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGGCGGGGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTa           <  2:213007/242‑1 (MQ=255)                                                                                                                                                                                                                                 tCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGGGTTAGCCg       >  1:2691/1‑251 (MQ=255)                                                                                                                                                                                                                                       tgcctaaTTTATGCCACGGCAGAGCCTCAACTGGGCCGAAGCGCATCCACGTGTCCG‑GTTAGGCTTTTCAAGTGGCGCCGGGCGCTGGCGCGCGTTGCCGTTTGCTACGAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCGGATGACCCAACGATTCGCGTCGGTGGGCCGACCTATCATTGGGTACGCTAAAGTATTCTGGCTGGCTAACAGGTGTTAGCCGGTGCg  <  2:27357/244‑1 (MQ=255)                                                                                                                                                                                                                                               ttATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACt                                                                                                       <  1:90105/141‑1 (MQ=255)                                                                                                                                                                                                                                               ttATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACt                                                                                                       >  2:90105/1‑141 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              GGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGAT‑GGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCG  >  NC_000913/4009455‑4009934

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