Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,009,349	C→T	C60C (TGC→TGT)	pldB →	lysophospholipase L2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,009,349	0	C	T	97.5%	148.6 / ‑6.5	40	C60C (TGC→TGT)	pldB	lysophospholipase L2
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (15/24);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCGTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGCCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAAT  >  NC_000913/4009205‑4009493                                                                                                                                                 |                                                                                                                                                 aaaCGCGTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGCCCgggg                                                                                                                                            >  2:618799/1‑151 (MQ=255)   aCGCGTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCg                                                                                                                                          <  2:360690/151‑1 (MQ=255)      cgTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTAt                                                                                                                                       <  2:90584/151‑1 (MQ=255)       gTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTAtt                                                                                                                                      <  2:99634/151‑1 (MQ=255)       gTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTAtt                                                                                                                                      >  2:373868/1‑151 (MQ=255)       gTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTAtt                                                                                                                                      >  1:86057/1‑151 (MQ=255)            ctgctTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTgag                                                                                                                                   >  2:262644/1‑149 (MQ=255)             tgctTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGagag                                                                                                                                 >  2:115877/1‑150 (MQ=255)                  ttACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACg                                                                                                                            <  2:744639/150‑1 (MQ=255)                             ccGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCgg                                                                                                                >  1:605591/1‑151 (MQ=255)                              cGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGa                                                                                                               <  1:486601/151‑1 (MQ=255)                                  gactgaTTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTg                                                                                                           >  2:439275/1‑151 (MQ=255)                                  gactgaTTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTg                                                                                                           >  2:320582/1‑151 (MQ=255)                                      gaTTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCt                                                                                                       <  2:287456/151‑1 (MQ=255)                                       aTTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGcc                                                                                                        <  1:97554/149‑1 (MQ=255)                                        tttCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTAt                                                                                                     >  1:226071/1‑151 (MQ=255)                                                  cAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTcc                                                                                           <  1:657919/151‑1 (MQ=255)                                                     cGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCAtt                                                                                        >  1:689525/1‑151 (MQ=255)                                                      gTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCAttt                                                                                       <  1:358284/151‑1 (MQ=255)                                                         atgaAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGagag                                                                                                                                 >  2:457013/1‑106 (MQ=255)                                                         atgaAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGagag                                                                                                                                 <  1:456962/106‑1 (MQ=255)                                                               cAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGa                                                                              <  2:163291/151‑1 (MQ=255)                                                                   gTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTc                                                                          <  1:52158/151‑1 (MQ=255)                                                                       aCTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTa                                                                       >  2:715139/1‑150 (MQ=255)                                                                             gtgGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAatcatc                                                                 <  2:624313/150‑1 (MQ=255)                                                                                   gACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATc                                                          <  2:493199/151‑1 (MQ=255)                                                                                        tCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCgg                                                      <  2:721106/150‑1 (MQ=255)                                                                                          cGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCa                                                   <  2:248364/151‑1 (MQ=255)                                                                                            gTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAgg                                                 <  1:731270/151‑1 (MQ=255)                                                                                              gcgcTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAggg                                                >  1:339912/1‑150 (MQ=255)                                                                                                gcTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACg                                             <  1:645063/151‑1 (MQ=255)                                                                                                             tCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCt                                 <  1:4438/150‑1 (MQ=255)                                                                                                                 gcACAGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAg                            <  1:654668/151‑1 (MQ=255)                                                                                                                   acaGCACCATGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGcc                          <  2:168135/151‑1 (MQ=255)                                                                                                                           aTGACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGc                   >  2:758330/1‑150 (MQ=255)                                                                                                                             gACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCAtc                >  1:523008/1‑151 (MQ=255)                                                                                                                             gACCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCAt                 <  2:761757/150‑1 (MQ=255)                                                                                                                              aCCGGGTGGTAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCAtc                <  1:11110/150‑1 (MQ=255)                                                                                                                                   gtggtAGTCATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGc          <  1:94150/151‑1 (MQ=255)                                                                                                                                           cATCTGTCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAAt  >  1:666717/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                 AAACGCGTTTGCTGCTTTTACCATGGGACCGCTGACTGATTTCTGGCGTCAGCGTGATGAAGCAGAGTTTACTGGTGTGGATGACATTCCGGTGCGCTTTGTCCGTTTTCGCGCACAGCACCATGACCGGGTGGTAGTCATCTGCCCGGGGCGTATTGAGAGCTACGTAAAATATGCGGAACTGGCCTATGACCTGTTCCATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAAT  >  NC_000913/4009205‑4009493

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