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	98.4%	219.2 / ‑4.4	63	A177T (GCA→ACA)	pldB	lysophospholipase L2
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/37);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

ATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGACGCCACGCCAACGCTTCTCTTGCAGGCTGAAGAGGAACGCGTGGTGGA  >  NC_000913/4009405‑4009990                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                     aTTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGANGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGc                                                                                                                                                                                                                                                                                                <  1:497999/300‑1 (MQ=255)       ggTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATcc                                                                                                                                                                                                                                                                                          <  1:453868/300‑1 (MQ=255)            gATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAAc                                                                                                                                                                                                                                                                                     <  1:52798/300‑1 (MQ=255)             aTGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACt                                                                                                                                                                                                                                                                                    <  1:412289/300‑1 (MQ=255)                 cTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGCTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGc                                                                                                                                                                                                                                                                                >  1:620009/1‑300 (MQ=255)                              cATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCa                                                                                                                                                                                                                                                                   <  2:721546/300‑1 (MQ=255)                                    gggCAGGGACGTTCCTGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTGAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATNTACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACACCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTc                                                                                                                                                                                                                                                             <  2:168140/300‑1 (MQ=255)                                     ggCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCNTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCAATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGNCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTcc                                                                                                                                                                                                                                                            >  2:683867/1‑300 (MQ=255)                                                  cGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGc                                                                                                                                                                                                                                               <  1:243516/300‑1 (MQ=255)                                                     tCGCCTGGTAGCCGATCCATATCTCGTGCATGTTAATCGCTTTAATGACTATGTTGATGATCTGGCGGCATTCTGGCAGCAGTAGGTTCAGTCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGAGNGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAAt                                                                                                                                                                                                                                            <  2:173867/300‑1 (MQ=255)                                                       gCCTGTNAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGNTAATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAg                                                                                                                                                                                                                                          <  2:590128/300‑1 (MQ=255)                                                        ccTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAgg                                                                                                                                                                                                                                         >  2:710973/1‑300 (MQ=255)                                                        ccTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCAGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAgg                                                                                                                                                                                                                                         >  2:275423/1‑300 (MQ=255)                                                           gTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCTGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTACTTCTGCAACGCCATCCAGGTGTATTTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCTTGATGGCTATGCAATAGGCAc                                                                                                                                                                                                                                      <  2:276632/300‑1 (MQ=255)                                                           gTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCAc                                                                                                                                                                                                                                      <  1:721003/300‑1 (MQ=255)                                                                ccGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGgc                                                                                                                                                                                                                                 <  1:609137/300‑1 (MQ=255)                                                                     ccGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTgg                                                                                                                                                                                                                            <  1:145045/300‑1 (MQ=255)                                                                            tcGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCCTTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCAGTGATGGCTATGCAATAGGCACCGAGCGCTGGCGCGCGt                                                                                                                                                                                                                     >  2:270162/1‑300 (MQ=255)                                                                                  aTGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGt                                                                                                                                                                                                               >  1:352620/1‑300 (MQ=255)                                                                                            cTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCaa                                                                                                                                                                                                     >  2:298982/1‑300 (MQ=255)                                                                                                  tGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACt                                                                                                                                                                                               <  1:44091/300‑1 (MQ=255)                                                                                                  tGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCATCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACt                                                                                                                                                                                               <  1:43414/300‑1 (MQ=255)                                                                                                         gggggggggggggcggcATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAGCGCGATATACTGGCATAGTCGGTGGGCGGTGCGATCTCCACGTTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCTTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAGCGTACTGACCcac                                                                                                                                                                                        <  2:388428/289‑1 (MQ=255)                                                                                                         gTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCcac                                                                                                                                                                                        <  2:388558/300‑1 (MQ=255)                                                                                                         gTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCcac                                                                                                                                                                                        >  1:630100/1‑300 (MQ=255)                                                                                                            gatgatTTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTTGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCAcagc                                                                                                                                                                                     <  2:245424/300‑1 (MQ=255)                                                                                                                 tCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACa                                                                                                                                                                                >  1:203364/1‑300 (MQ=255)                                                                                                                   tGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGGAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGc                                                                                                                                                                              <  2:562808/300‑1 (MQ=255)                                                                                                                     gcggcATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGa                                                                                                                                                                            <  1:699477/300‑1 (MQ=255)                                                                                                                     acggcATTCTGGCAGCAGGAGGTTCAGCCCGNTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCANTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCACCGTACTGACCCACAGCAGACCGCGa                                                                                                                                                                            >  2:521806/2‑300 (MQ=255)                                                                                                                          aTTCTTGCATCAGGATGTTCGGCCCGGTTCGTGGCGTAAATGCTCTATACTGGCACATTCGATGGGCGGTGCGAGCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAATTGCGCCAATGGTGGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCTTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTCACCCACAACAGACACCGATATca                                                                                                                                                                       <  2:467592/300‑2 (MQ=255)                                                                                                                                          gTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTc                                                                                                                                                       >  1:363210/1‑300 (MQ=255)                                                                                                                                               gCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGc                                                                                                                                                  <  1:576789/300‑1 (MQ=255)                                                                                                                                                    gTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtgatg                                                                                                                                             <  1:406033/300‑1 (MQ=255)                                                                                                                                                      ccGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGTGATCTCCACATTATGTCTGCAACGCCAGCCAGTTGTATGTGACGCCATTGCGCTAACGGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGAc                                                                                                                                           <  2:398714/300‑1 (MQ=255)                                                                                                                                                           gCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAAc                                                                                                                                      >  1:596271/1‑300 (MQ=255)                                                                                                                                                                 aCGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTcg                                                                                                                                >  2:456211/1‑300 (MQ=255)                                                                                                                                                                 aCGCTATATACTGACACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTCTGCTGATGACCCAACGATTcg                                                                                                                                >  2:723005/1‑300 (MQ=255)                                                                                                                                                                      atataCTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTCTGACACGGCAGATCCTCAACTGGGACGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCCCCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTAATGGCCCACAGCGGACAGCGATATCGAAGTAACTTACGCTTCTCTGCTGATGActacacgctttgcgtcg                                                                                                                           >  2:10712/1‑284 (MQ=255)                                                                                                                                                                                   ttCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTa                                                                                                              >  1:131047/1‑300 (MQ=255)                                                                                                                                                                                       aTGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCAt                                                                                                          <  1:709876/300‑1 (MQ=255)                                                                                                                                                                                       aTGGGCGGTGCGATCTCCACATTATTTCGGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATTACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCAt                                                                                                          <  2:711430/300‑1 (MQ=255)                                                                                                                                                                                           gtggTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGTGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTggg                                                                                                      <  2:302749/298‑1 (MQ=255)                                                                                                                                                                                              gTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGCTGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTAc                                                                                                   >  1:276268/1‑300 (MQ=255)                                                                                                                                                                                              gTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTAc                                                                                                   >  1:276812/1‑300 (MQ=255)                                                                                                                                                                                                     ctcCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAg                                                                                            >  1:738015/1‑300 (MQ=255)                                                                                                                                                                                                         acaTTATTTCTGCAACGCCATCCAGGTGTATNTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCNTCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGAACTACCATTGGGTACGCGAAAGTAtt                                                                                        >  2:514099/1‑300 (MQ=255)                                                                                                                                                                                                          caTTATTTCTGCAACGCCATCCAGGTGTATTTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTc                                                                                       <  2:29797/300‑1 (MQ=255)                                                                                                                                                                                                                        cGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAg                                                                         <  1:325231/300‑1 (MQ=255)                                                                                                                                                                                                                            aTCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGgtgt                                                                     >  1:417387/1‑300 (MQ=255)                                                                                                                                                                                                                                aGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCTAAGCGCATCCACGTTTCCTTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGc                                                                 <  1:365437/300‑1 (MQ=255)                                                                                                                                                                                                                                 ggTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGcc                                                                <  1:662587/300‑1 (MQ=255)                                                                                                                                                                                                                                      aTGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGc                                                           >  1:289150/1‑300 (MQ=255)                                                                                                                                                                                                                                      aTGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTATATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGc                                                           >  1:286334/1‑300 (MQ=255)                                                                                                                                                                                                                                                aTTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCTTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGac                                                 <  2:411898/300‑1 (MQ=255)                                                                                                                                                                                                                                                aTTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGac                                                 <  1:362924/300‑1 (MQ=255)                                                                                                                                                                                                                                                             cgcCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGACGCCACGCCAACGc                                    <  1:384244/300‑1 (MQ=255)                                                                                                                                                                                                                                                                    gTTTGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGACGCCACGCCAACGCTTCTCtt                             >  1:7245/1‑300 (MQ=255)                                                                                                                                                                                                                                                                       tGGGATCGTGATTCGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGACGCCACGCCAACGCTTCTCTTGCa                          <  1:159266/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                    cGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGACGCCACGCCAACGCTTCTCTTGCAGGCTGAAGAGGaa             <  1:396649/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                    cGTATGCCGTCATTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGACGCCACGCCAACGCTTCTCTTGCAGGCTGAAGAGGaa             <  2:716231/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                               aTTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGACGCCACGCCAACGCTTCTCTTGCAGGCTGAAGAGGAACGCGTGGTGGa  >  1:93411/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                               aTTTATGACACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGACGCAAAGCCAACGCTTCTCTTGCAGGCTGAAGAGGAACGCGTGGTGGa  >  2:230182/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                     ATTTGGGGTTTGATGTCTTAATCATCGACCATCGCGGGCAGGGACGTTCCGGTCGCCTGTTAGCCGATCCGCATCTCGGGCATGTTAATCGCTTTAATGATTATGTTGATGATCTGGCGGCATTCTGGCAGCAGGAGGTTCAGCCCGGTCCGTGGCGTAAACGCTATATACTGGCACATTCGATGGGCGGTGCGATCTCCACATTATTTCTGCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAGCCGGTGCGGGTGATGACGCCACGCCAACGCTTCTCTTGCAGGCTGAAGAGGAACGCGTGGTGGA  >  NC_000913/4009405‑4009990

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