Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	15,174	G→A	R336K (AGG→AAG)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,174	0	G	A	98.1%	168.1 / NA	53	R336K (AGG→AAG)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (22/30);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAAC  >  NC_000913/14942‑15402                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      gAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCCCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAAGCGTCAAGTCTGTCCGCGGGGGCGCACAGGGCGATTTGCTGTGCCGCGTTGTCGTCGAAACCCCGGTCGGCCTGTACGCGAAGCAGAAACAGCcgctgc                                                                                                                                                                                                                      >  1:188790/1‑250 (MQ=255)   aGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaa                                                                                                                                                                                                                    >  1:435859/1‑250 (MQ=255)      caacaaCCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAGGCAGAAACAGCTgctgcaagagc                                                                                                                                                                                                                >  1:331800/1‑251 (MQ=255)         caacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgc                                                                                                                                                                                                             >  1:128810/1‑251 (MQ=255)              tctATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCCCAGGGGGATTTGCTGTGCCGCGTTGTCGTCGCAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgcaaga                                                                                                                                                                                                         >  2:504419/3‑250 (MQ=255)                           cccGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGc                                                                                                                                                                                           <  1:393782/250‑1 (MQ=255)                                          tATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGATTTGCTGTGCCGCGTTGTCGTCGGAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCCACCggtgggcac                                                                                                                                                                            >  1:106283/1‑243 (MQ=255)                                          tATGGCGGCGCTGGGGGGCGCAATCGAAGTACCCACCCTTGTTAGGCGCGTCCCACTGTAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGAGCGGTAAAGGCGGCAAGGCTGTCCGCGGGGGCGCACAGGGGGGTTTGCAGTTCCGCGGTGTCGGCGACACCCCGGGAGGCCTGAACGAGGAGCAGAAACAACTGCTGCAAGAACTGCAAGAAAGATTCGAGTGCCCCACCCggggacacc                                                                                                                                                                            >  2:263975/1‑244 (MQ=255)                                             ggcggcGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTTAAAGGGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCacaaca                                                                                                                                                                        <  1:138933/251‑1 (MQ=255)                                               cggcgCGGGGTGGCGAAATCGAAGTACCGACCCGGGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGc                                                                                                                                                                      <  2:284245/251‑1 (MQ=255)                                                 gcgctggtgggggcattCGAAGTACCGACCCTTGATGGTCGCTTCCATCTGAAAGTTCCGGGCGAACCCCATACCGGTAATCTTTTCCGGATGCGCGGTAAAGGCGTCAAGTGTGTCCGCGGTGGCGCACAGTGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCATAAACATCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGTCCCAACCGGTTAGCACCACAGccc                                                                                                                                                                    <  2:262905/235‑1 (MQ=255)                                                 gcgcTGGGTGGCGAACTCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGccc                                                                                                                                                                    <  1:261103/251‑1 (MQ=255)                                                       ggTGGCGAATTCGAAGTACCGACCCTTCATGGTCGCGTCAAACCGAAAGGGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTTGCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTc                                                                                                                                                              <  2:344743/251‑1 (MQ=255)                                                        gTGGCGAAAGCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCa                                                                                                                                                             <  2:422052/251‑1 (MQ=255)                                                            cGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAGGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAaga                                                                                                                                                         <  1:38990/251‑1 (MQ=255)                                                                   gAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGcttctt                                                                                                                                                  >  2:124476/1‑251 (MQ=255)                                                                      gTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGa                                                                                                                                               <  2:188178/251‑1 (MQ=255)                                                                         ccGACCCTTGAGGTTCGCGTCAAACTGAAAGTGCCTGGAGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGGAAATGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGGAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATgg                                                                                                                                            <  2:476917/251‑1 (MQ=255)                                                                            gcccTTGATGGCCGCCGAAAACGGAAATTGCTTGGAGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGgtgt                                                                                                                                         <  2:137942/250‑1 (MQ=255)                                                                            aCCCTTGGTGGTCGCGTCAAACTGACAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAGACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGgtgg                                                                                                                                         >  2:247032/1‑250 (MQ=255)                                                                            aCCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGgtgg                                                                                                                                         >  2:31475/1‑250 (MQ=255)                                                                             cccgtGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACACAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGtgtg                                                                                                                                        <  1:401246/247‑1 (MQ=255)                                                                             cccTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGtgtg                                                                                                                                        <  2:331800/251‑1 (MQ=255)                                                                                 tGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGGCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTgaa                                                                                                                                      <  1:265518/249‑1 (MQ=255)                                                                                 cgATGGACGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCTCAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTgaag                                                                                                                                     <  1:24310/249‑1 (MQ=255)                                                                                   aTGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAACGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTgaaga                                                                                                                                    >  2:72484/1‑249 (MQ=255)                                                                                    tGGTCGCGTCAAACTGAAGGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGt                                                                                                                                 >  1:123011/1‑251 (MQ=255)                                                                                     ggTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTTTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGGTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGt                                                                                                                                 <  1:59544/250‑1 (MQ=255)                                                                                                 cTGAAAGTGCCTGGAGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCAAACCGGTGAGCACAACAGACAGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCt                                                                                                                    >  2:466806/1‑251 (MQ=255)                                                                                                               cGAAACCCAGACCGGTAAGCTTTTCCGTATGAGCGGTAAAGGCGTCAAGTCTGGCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTTATGGTGTGAAGAAGTTTTTAGACGACCTGACCCGCTAACCTc                                                                                                      <  2:84457/251‑1 (MQ=255)                                                                                                                 aaaccaatgccGGTAAGCTTTGCAGTATGCGCGGTCAAGGCGTCAAGGCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGGGCCGCGTGGTCGGCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGATGCAAGAAAGCGTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGGTTTCTGACGACCTGACCCGTTAACCTccc                                                                                                    <  2:72753/242‑1 (MQ=255)                                                                                                                  aaCCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTcccc                                                                                                   >  1:48852/1‑251 (MQ=255)                                                                                                                    cccAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCa                                                                                                  >  1:403648/1‑250 (MQ=255)                                                                                                                        gACCGGTAAGCTGTGCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCTCGCGCTCAAAGAGCTTCTTTGATGGTTTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGc                                                                                             <  1:297420/251‑1 (MQ=255)                                                                                                                                aGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAACCTGCCCgg                                                                                     >  2:265387/1‑250 (MQ=255)                                                                                                                                    tttCCGTATGCGCGGTAAAGGCGGCGAGTCGGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGAGGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTgg                                                                                   <  2:338085/248‑1 (MQ=255)                                                                                                                                       ccggaTGCGCGGGAAAGGCGTCAAGTCTGTCCGCGGGGGCGCACAGGGTGAGTTGCTGTGCCGCGTTGGCGTCGAAACACCGGTAGGCCAGAACGAGAAGCAGAAACAGCTGCTTCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTGTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAg                                                                               <  2:312258/246‑1 (MQ=255)                                                                                                                                                      aGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAACGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTaaaatt                                                                 >  1:48785/1‑247 (MQ=255)                                                                                                                                                         cGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAgggg                                                            >  1:412761/1‑250 (MQ=255)                                                                                                                                                            cAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGc                                                          <  1:72484/250‑1 (MQ=255)                                                                                                                                                              aggcTGTCAGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGATCTTCTTTTATGGTGTGAAGAAGTTTTTTGACGTCCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGTCTTGGGTAAAATTAGGGTGCGtt                                                       <  2:430518/248‑1 (MQ=255)                                                                                                                                                                 cTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCACGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGa                                                     >  1:430310/1‑250 (MQ=255)                                                                                                                                                                          ggggCGCACAGGGTGATTTGCGGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGa                                           <  2:22421/249‑1 (MQ=255)                                                                                                                                                                                      gTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGCTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTgaaa                                 >  1:52118/1‑246 (MQ=255)                                                                                                                                                                                      gTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTgaaa                                 >  2:175414/1‑246 (MQ=255)                                                                                                                                                                                                gtgCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGAt                     <  2:106283/251‑1 (MQ=255)                                                                                                                                                                                                  gCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCa                   <  1:175414/251‑1 (MQ=255)                                                                                                                                                                                                         tGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCCAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACGccccc            >  1:58386/1‑249 (MQ=255)                                                                                                                                                                                                              gcgAAACCCCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCTCTCCCATAAgcg       <  2:159008/250‑1 (MQ=255)                                                                                                                                                                                                                 aaaCACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTa    <  1:436699/251‑1 (MQ=255)                                                                                                                                                                                                                   acacCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAAc  <  2:58386/251‑1 (MQ=255)                                                                                                                                                                                                                      ccGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAg                                                                               >  2:501051/1‑171 (MQ=255)                                                                                                                                                                                                                      ccGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAg                                                                               <  1:501051/171‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      GAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCG‑GTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAAC  >  NC_000913/14942‑15402

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