Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,558,907	T→A	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,558,907	0	T	A	94.2%	144.9 / ‑5.9	53	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base A (25/24);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCC  >  NC_000913/4558663‑4559148                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  cccGTTGCCGTCGGAGCTGAGGGGGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATGCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                                                                                                                                                                                                             <  1:272144/250‑1 (MQ=255) aCCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                                                                                                                                                                                                             >  1:15655/1‑251 (MQ=255)    gTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATCCCCGCCTTAACGGCGCGGGCAATGCCGTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                                                                                                                                                                                                          <  1:212310/251‑1 (MQ=255)        ccGTCGGAGCGGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                                                                                                                                                                                                      <  1:232618/251‑1 (MQ=255)                  tGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                                                                                                                                                                                                                            >  1:70181/1‑251 (MQ=255)                    aGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                                                                                                                                                                                                          <  1:100886/251‑1 (MQ=255)                    aGGGTGACGCGAGCCAGCGGCATACCCCCCTGAACGGCGCGGGCAATCCCTTAGGCAGGGGCGACCGGTTCGTCAATGCCGCTGGTGATCTCGATCGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTCCGTTGCGGTTCACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCCGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGCAacac                                                                                                                                                                                                                          >  2:343915/1‑251 (MQ=255)                          aCGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGAc                                                                                                                                                                                                                    <  2:90443/251‑1 (MQ=255)                              gTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                                                                                                                                                                                                 >  1:52504/1‑250 (MQ=255)                              gTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTAGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                                                                                                                                                                                                 >  1:218816/1‑250 (MQ=255)                                       gAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg                                                                                                                                                                                                       <  1:107214/251‑1 (MQ=255)                                          tACCCGCCTGAAAGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGCACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg                                                                                                                                                                                                    >  2:41297/1‑251 (MQ=255)                                             ccGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                                                                                                                                                                                                 <  1:126612/251‑1 (MQ=255)                                             ccGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                                                                                                                                                                                                 >  1:306318/1‑251 (MQ=255)                                              cGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTGGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGTTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATATGCTGTAACGCCATTTTTCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                                                                                                                                                                                                <  1:227520/251‑1 (MQ=255)                                                       ggCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCCTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGATTACCGCCGAGCAAACCGCCa                                                                                                                                                                                        >  1:196464/1‑250 (MQ=255)                                                        gcgcGGGCAATCCCTCAGGCAGGGGCGACCGGTGCGACAATGCTGTTGGTAATGTCGATGTTTCCGCCTTTGCGCGCGAACTCCAGTGCTTGTCCAAGCAACTGTAGGTTGCGGTTAACGGTGGTCGGCAGCAGCTGGCTAAGCGGCACAGCGCAGTTTTCCAGCAGGTTATAAATAGGCTGTAACTCTTTTTTACTGTCGCGCTGGTTGAACGCGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAc                                                                                                                                                                                      <  1:163221/251‑1 (MQ=255)                                                           cgGGCAATCCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCGGGTTTACCGCCGAGCAAACCGCCAAcgcg                                                                                                                                                                                   <  2:30247/251‑1 (MQ=255)                                                           cgGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgcg                                                                                                                                                                                   >  2:92498/1‑251 (MQ=255)                                                                   ccctccGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACAGCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCt                                                                                                                                                                            <  2:52504/245‑1 (MQ=255)                                                                           cAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCat                                                                                                                                                                    >  2:128095/1‑250 (MQ=255)                                                                           cAGGGGCGACCGGTTCGTCAATGCTGCGGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACACCGGTCCGTTGCGGTTAACGGGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATTTGGAACACGGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCata                                                                                                                                                                   <  1:171011/251‑1 (MQ=255)                                                                                  gACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGACAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGcc                                                                                                                                                             >  1:65382/1‑250 (MQ=255)                                                                                     cGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAACCAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAggg                                                                                                                                                         >  2:207689/1‑250 (MQ=255)                                                                                                        gTGAGATCGATGTGGCCGCCTTTGCCCGCGACCTCCAGCGCCTGCTCAAACAACGGTCCGTTGGGGTTAACGTGGGTCGGCAGCAGCTTGCTGATGGGAACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTCCTGCCGCCCATGTGGAACCCGGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATCCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGc                                                                                                                                      <  2:245109/251‑1 (MQ=255)                                                                                                          gATAGCGATGGTGCCGCCTTTGCCCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTTTCGCCCATGTGGAACACGGTGGCGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGAGTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAg                                                                                                                                    <  1:238873/251‑1 (MQ=255)                                                                                                                 aTGGTGCCGCCTGCGCACGCGCAATCCAGCGCCTGCTCAAACAACGGTCCGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGAGCGGCAGAACGGTg                                                                                                                             >  2:26658/1‑251 (MQ=255)                                                                                                                    gTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTTTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                          <  1:207689/251‑1 (MQ=255)                                                                                                                      gccgcctcagcgagcgaccaccAGAGCCTGCTCAAACAACGGTACCTTGCGGTTAACGTGGGTCGGCAGCAACTTGCTGATCGGCACATCGCAGTTTTCCCGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGGACGGGGATc                                                                                                                          >  2:341343/21‑249 (MQ=255)                                                                                                                          ccTTGGCGCGGGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGGGGGTCGGCAGCAGCTTGCTGAGCGGCACCTCGCAGTTTTCCAGCAGGTCATAAATGGGCTGGACCGCCAGTTTCCTGCCGCCCATGGGGACCCCGTTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGAt                                                                                                                    <  2:151399/251‑1 (MQ=255)                                                                                                                                    cgAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGCTTCTGCCGCCATATTGGCCAGGCGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTc                                                                                                          >  2:223091/1‑251 (MQ=255)                                                                                                                                           cAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCa                                                                                                    >  2:128933/1‑250 (MQ=255)                                                                                                                                             gcgcCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCAAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATc                                                                                                 <  1:258673/251‑1 (MQ=255)                                                                                                                                                   gCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTCCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCacac                                                                                             <  1:147550/249‑1 (MQ=255)                                                                                                                                                       aaaCAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCCCATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatcga                                                                                       <  1:13980/251‑1 (MQ=255)                                                                                                                                                          caacGGTACGTTGCGGTTAACGTGGGTAGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCCGGTGATAAACGTCCGGTGCGGCCGAACGGTGATCCGAGATGGCGCCTTTTACGCCacgcaccagatcgtat                                                                                     >  2:1630/1‑235 (MQ=255)                                                                                                                                                               gTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGcc                                                                               >  1:223328/1‑251 (MQ=255)                                                                                                                                                               gTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATAGACACATCGAAGTTTTCCAGCAGGTCATAAATAGGCTGTAAAGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGACAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCAACCGATCGATAATc                                                                                  >  2:129411/1‑248 (MQ=255)                                                                                                                                                               gTACGTTGCGGTTAACGGGGGTCGGCAGCAGTTTGCTGATCGGCACATCGCAGTTTCCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACGGTCGCCCATGTGGAACACGGTGACGCCAGGTTTGCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATc                                                                                  <  1:129411/248‑1 (MQ=255)                                                                                                                                                                        ggTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCtttt                                                                      >  1:126709/1‑251 (MQ=255)                                                                                                                                                                         gTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCttttt                                                                     >  1:186892/1‑251 (MQ=255)                                                                                                                                                                         gTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGGGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCACTCACAAGATCGCTAATCGCCACGTCttttt                                                                     >  2:15111/1‑251 (MQ=255)                                                                                                                                                                           tatcGTGGGTCGGCAGCGGCTGGCTGATCGGCACATCGCAGTTTTCCAGCAGGCCATAAATAGGCGGTAACGCCATTTTACGGTCGCCCATGTGAAACACGGTGACGCCAGGTTTACGGCCGAGCAAACCGCCAACTCGGGTTTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTc                                                                    <  1:171047/247‑1 (MQ=255)                                                                                                                                                                             aCGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCAc                                                                 >  1:149731/1‑251 (MQ=255)                                                                                                                                                                                  ggTCGTCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTACCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAc                                                            <  1:230170/251‑1 (MQ=255)                                                                                                                                                                                                          cacaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACAt                                    >  1:96878/1‑251 (MQ=255)                                                                                                                                                                                                              gcgcatgtgcccAGAAGGTCAAAAATAGGTGGAAACGCCATTTTACTGTCGCCATTGGGGAACGCGGTGACGCCAGGTTTCCCGCCGACCAACCCGCAACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATa                                <  2:15655/241‑1 (MQ=255)                                                                                                                                                                                                                         cagcagGTCATAAATAGGCTGGAACGCCATTTTACTGTCGCCCATGGGGAACAGGGTGCCGCCATGTTTCCCGCCGAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCTGGAAGGGACATGATAAGCGCCGGTCa                     <  2:223328/251‑1 (MQ=255)                                                                                                                                                                                                                             agGTCATAAATAGGCGGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCAt                 <  1:128095/251‑1 (MQ=255)                                                                                                                                                                                                                                      aaggtgttaaccGCCTTTTGACGGGCGCCCAGGGGGAACGCGGTGCCGCCAGGTGTACCGCGGAGCAAACCCCCAACGCGGGATCCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAAGCACACGATCGATAATCGCCGCGTCTTTTTCCACGGAACCCGTAAGGGTGCGGGAAGGGACATGATAAGCGCCGGGCAGCATCCAGGCGCt        <  2:148367/240‑1 (MQ=255)                                                                                                                                                                                                                                           cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGc   >  2:307306/1‑251 (MQ=255)                                                                                                                                                                                                                                            tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGcc  >  1:264305/1‑251 (MQ=255)                                                                                                                                                                                                                                             gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGcc  <  1:223091/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  ACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCC  >  NC_000913/4558663‑4559148

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