Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,386	C→T	S12S (TCC→TCT)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,386	0	C	T	96.4%	230.2 / ‑0.2	84	S12S (TCC→TCT)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (36/44);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

CAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTA‑TTGGCTGCT‑TGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAA  >  NC_000913/71139‑71625                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTAcc                                                                                                                                                                                                                                                   >  2:817870/1‑251 (MQ=255)  aaTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATCATTGGAAACGAATCAGACAATTGACGGCTTGACGGAGACGCCTAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACCTTCTGCACGCATTGCTGGAACACTTTATTACACAATCTACCg                                                                                                                                                                                                                                                  >  2:917478/1‑251 (MQ=255)  aaTGTTGGTTTTGACGATCAACGTTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTTGAACACTTTATTACCCAATCTACCg                                                                                                                                                                                                                                                  <  2:712257/251‑1 (MQ=255)   aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGACTGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATACCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTCTTACCCAATCTACCgt                                                                                                                                                                                                                                                 >  2:84598/1‑251 (MQ=255)               aCGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGAt                                                                                                                                                                                                                                     <  1:618891/251‑1 (MQ=255)               aCGATCAACTCTAGTTCGCGCTGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGATTTCCGTGCCGGTTGTGAATAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTTACGTAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCAtttat                                                                                                                                                                                                                                     <  1:618886/251‑4 (MQ=255)                cGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCACTCTACCGTGTATTCATTGATg                                                                                                                                                                                                                                    >  2:507518/1‑251 (MQ=255)                    cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCgg                                                                                                                                                                                                                                <  1:814199/251‑1 (MQ=255)                     aaCTCTATTTCTCGCGGGTATTTAAAACATGCACCGGTGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACTAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCggt                                                                                                                                                                                                                               <  2:663476/251‑1 (MQ=255)                       ctctATTTCTCGCGGGTATTTAAAACATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCggtgg                                                                                                                                                                                                                             >  1:693328/1‑251 (MQ=255)                           aTTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCggtggt                                                                                                                                                                                                                            <  1:491342/248‑1 (MQ=255)                           aTTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAAGCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCggtggt                                                                                                                                                                                                                            <  2:289789/248‑1 (MQ=255)                                   cgggTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCt                                                                                                                                                                                                                 >  1:588197/4‑251 (MQ=255)                                   cgggTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGAATTGCTGGACCACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCt                                                                                                                                                                                                                 >  2:390593/4‑251 (MQ=255)                                    gggTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCtt                                                                                                                                                                                                                >  1:24023/3‑251 (MQ=255)                                     ggTATTTAAAAAAGGCACCGGGGCCAGCCCGAGCGAGGTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATTCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCttt                                                                                                                                                                                                               <  2:6988/250‑1 (MQ=255)                                      gTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCtttt                                                                                                                                                                                                              >  1:870156/1‑250 (MQ=255)                                            aaaaaaTGCACCGGGGCCAGCCCGGGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTGGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGt                                                                                                                                                                                                        <  1:737033/251‑1 (MQ=255)                                             aaaaaTGCACCGGGGCCAGCCCGAGCGCGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGGGTc                                                                                                                                                                                                       >  1:361317/1‑251 (MQ=255)                                               aaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTcgc                                                                                                                                                                                                     >  1:811334/1‑251 (MQ=255)                                                 aTGCGCCGGGGCTAGCCCGAGCGAGTTCCGTGCCGGTTGTGAATAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTTTTGGTGGCCTTTCTGGATTcgctg                                                                                                                                                                                                   <  2:770741/251‑1 (MQ=255)                                                 aTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGGGTTGGGGGCCTTTCTGGAGTcgctg                                                                                                                                                                                                   >  1:838870/1‑251 (MQ=255)                                                       cGGGGCCAGCCCGTGCGAGTTCCGTGCCGGGTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTcgctggcgctg                                                                                                                                                                                             <  2:147989/251‑1 (MQ=255)                                                         gggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTcgctggcgctggg                                                                                                                                                                                           >  1:747747/1‑250 (MQ=255)                                                              aGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGAGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTTGGCGGtt                                                                                                                                                                                      >  1:266521/1‑251 (MQ=255)                                                                cccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTg                                                                                                                                                                                    <  2:811334/251‑1 (MQ=255)                                                                  cgagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGAt                                                                                                                                                                                  >  1:935048/1‑251 (MQ=255)                                                                   gagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGAtt                                                                                                                                                                                 >  1:75641/1‑251 (MQ=255)                                                                     gggagtgccGTGCCGGTTGTGAAGAAAAAGTTAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGCCAGGCACATTATGCAAGCATTGCTGGAACTCTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTTTTGGTGGCCTTTCTGGAGTTGCTGGCGCTGGTCGGTTTGATTCt                                                                                                                                                                               <  2:609489/243‑1 (MQ=255)                                                                     gcgaggtccgggCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCt                                                                                                                                                                               <  2:90757/240‑1 (MQ=255)                                                                     gcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGGTTGATTCt                                                                                                                                                                               >  1:47003/1‑251 (MQ=255)                                                                      cgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTa                                                                                                                                                                              <  1:618818/251‑1 (MQ=255)                                                                       gagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGATGGCGGTGGTAGGTTTgtatctac                                                                                                                                                                             >  1:252988/1‑244 (MQ=255)                                                                                gCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACgg                                                                                                                                                                    >  1:914917/1‑251 (MQ=255)                                                                                gCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGA‑GGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTTGGGGCCTTTCTGGAGTCGCTGGCGCTGTGCGGGTTGATTCTACCCGGTAcagg                                                                                                                                                                   >  2:271930/1‑248 (MQ=255)                                                                                 ccGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACAATTTATTACCCACTCTACCGTGTATTCATTGATGGCGGTGGTGTTGGGGGCATTTCTGGAGTCGCTGGCGCTGGTCGGGTTGATTCTACCCGGTACGGt                                                                                                                                                                   >  1:585359/1‑251 (MQ=255)                                                                                    gTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGGGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGGACGGGGCt                                                                                                                                                                >  2:410903/1‑251 (MQ=255)                                                                                         gaagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATgg                                                                                                                                                           >  1:491225/1‑251 (MQ=255)                                                                                         gaagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGGAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATg                                                                                                                                                            <  2:809574/250‑1 (MQ=255)                                                                                          aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCgt                                                                                                                                                                                                                                                 >  1:410144/1‑164 (MQ=255)                                                                                          aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTAGTACCCAACCCACCgt                                                                                                                                                                                                                                                 <  2:410144/164‑1 (MQ=255)                                                                                           agaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCg                                                                                                                                                         <  1:910336/251‑1 (MQ=255)                                                                                             aaaaaGTGAATGATTTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTTACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTTACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCggg                                                                                                                                                       <  1:557859/251‑1 (MQ=255)                                                                                                aattGAATGATGGAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGAACATTGACGGAGTGACGGAGGAGTATTGGGTTTGCAGAATCCCTG‑TTCGGTCCATTTGACAGGCACAGTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTTATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTTCTGATGGCGGGGCt                                                                                                                                                    <  2:228556/248‑1 (MQ=255)                                                                                                   tGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TTCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGGTTGATTCTACCCGGTACGGGGCTGATGGCGGGGGTggg                                                                                                                                                 >  1:924766/1‑251 (MQ=255)                                                                                                      atgttgTAGCCTTCAAGTTTTCATAATGGGTAACGAATCAGACAATGGACGGCTTGACGGAGTAGCTTAGGTTTTGCAGAATCCCTGCGTCG‑TCCATGTGACAGGCACATTATGCAAGCATTGCTGGAACTCTGTATTACCCAATCTTCCGTGTAGTCATTGATGGCGGTGGTGTTGGTGGCCTTTCGGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAgt                                                                                                                                              <  2:636770/247‑2 (MQ=255)                                                                                                        gatgTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGGACGGTGCTGATGGGGGGGCggggggggc                                                                                                                                            >  1:109175/1‑246 (MQ=255)                                                                                                        gatgGAGCCGTCAAGTTGTCATAATTGGTAACGAATGAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCATGCATTGCTGGAGCACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTTGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAgcg                                                                                                                                             <  1:908978/250‑1 (MQ=255)                                                                                                             aGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGAtt                                                                                                                                       <  1:555233/251‑1 (MQ=255)                                                                                                                   cAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGGGCGCTGATTgggagg                                                                                                                                 >  1:225136/1‑247 (MQ=255)                                                                                                                   cAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACTGTGTATTCATTGATGGCGGTGGGGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGGGTTTACCCGGTACGGTGCTGATGGCggggggcggggggcggtattggagg                                                                                                                                 >  2:141123/1‑230 (MQ=255)                                                                                                                      gTTGTAATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTTGAGTCGATGGCGCTGGTCGGTTTTATTCTACCCGGTACGGGGCTGATGGCGGGGGTgggggagtgggtttggcggggg                                                                                                                              >  2:903995/1‑232 (MQ=255)                                                                                                                          tCATAATTGGTAACGAATCAGACAATTGACGGCTTGCCGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGGTGGGGGGGCTGATTgggccgggggagg                                                                                                                          >  1:748280/1‑240 (MQ=255)                                                                                                                                   gTAACGAATCAGACAATTTACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCt                                                                                                               <  2:75641/251‑1 (MQ=255)                                                                                                                                    tAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGGTGGGAGCGCTGATTTGCAGCGGGGAG‑TTTAAGttttt                                                                                                               >  1:883500/1‑249 (MQ=255)                                                                                                                                    tAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGGGTTGGTGGCCTTTCTGGAGTCGCTGGCGATGGTCGGGTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGGGCGCTGGTTGGGCGCGGGGGG‑TT‑AAGATTCTg                                                                                                              >  1:915667/1‑251 (MQ=255)                                                                                                                                     aaCGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTg                                                                                                              <  1:817870/250‑1 (MQ=255)                                                                                                                                       cGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTAACCGGTACGGTGCTGATGGCGGGGGTGGGAGCGGTGGTTggaaacgggggagtaagatttttgg                                                                                                              >  1:972182/1‑227 (MQ=255)                                                                                                                                            cAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGGGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGGCGGTTTGATTCTACCCGGGACGGTGCTGATGGCGGGGCTGGGGGGGCTGATTTGCAGCGGCGAG‑TT‑AAgattttgggccgccc                                                                                                      >  1:314580/1‑237 (MQ=255)                                                                                                                                                  aTTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCctggctgg                                                                                                <  1:314521/251‑1 (MQ=255)                                                                                                                                                        ggCTTGACGGAGTAGCATATGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACTCTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGa                                                                                          <  1:794977/251‑1 (MQ=255)                                                                                                                                                            tGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCCCATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGGTGGTGGCCTTTCTGGAGTCGCTGGCGCGGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCGGGCTGGCAGGGatt                                                                                        <  2:870156/249‑1 (MQ=255)                                                                                                                                                               cGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGGGCGCTCATTGGCCGCGGCGGG‑GT‑TAGTTTTTGGGACGCCCGGCGGGCAggggttgtggg                                                                                   >  1:657920/1‑243 (MQ=255)                                                                                                                                                                ggAGTAGCATTGGGGTTGCAGACTCCCTGCGTCG‑TCCATTTGACAGGTCCATTAGGCAAGCATTGCTGGACAACGTTAGTTTCCAATCTACCGTGTATTTATTGATGGCGGTGGTGTTGGTGGTCGTTCTGGAGTCGCTGGCGCTGGTAGGTTTGATGCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGTTGGCAGGGattgtcgg                                                                                   <  2:28749/250‑6 (MQ=255)                                                                                                                                                                          aGGGTTTGCAGAATCCCTTCTTCG‑TCCATTTTACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGattgtcggttgtctgctg                                                                         <  2:747747/250‑16 (MQ=255)                                                                                                                                                                             gTTTGCAGAATCCCTGCTTCG‑TCCATTTGCCAGGCACATTAGGCAAGCATTGCGTGAACTCTTTATTACCCAAGCTACCGTGTAGTCATTGATGGCGGTGGTGGTGGTGGCCTTTCTGGAGTCGCTGGCGCGGGTCGTTTTGAGTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGTGCGGATTGGCAGCGGCGAG‑TT‑AGGTTTCTGGCACGCCTGGGTGGCAGGGattgtcggttgtcggctgggcg                                                                     <  2:885068/251‑20 (MQ=255)                                                                                                                                                                                                    ccATTTGATAGGCACATTATGCATGCATTGCAGTAACTATTTATTACCCAATCTACCGTGTATTCAGTGATGGCGGTGGTGTTGGTGGCCGTTCTGGAGGCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGATGGTAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGTCTGGATCTCCTTctggatgg                                              <  2:517905/250‑1 (MQ=255)                                                                                                                                                                                                    ccATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGACTGGATCTCCTTctggctgg                                              <  1:507518/250‑1 (MQ=255)                                                                                                                                                                                                         tgtcAGGCACATTATGCAAGCATGTCTGGAACACTTTATTACCCAACCTACCGTGTATTCATGGATGGCGGTGGTGTTGGTGGACTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATGGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGACTGGATCTCCTTCTGgctgggctgg                                         <  2:266521/247‑1 (MQ=255)                                                                                                                                                                                                         tGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTg                                                                                                                                                                                    <  1:191935/115‑1 (MQ=255)                                                                                                                                                                                                         tGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGGGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTg                                                                                                                                                                                    >  2:191935/1‑115 (MQ=255)                                                                                                                                                                                                            cAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGACTGGATCGCCTTCTGGCTGGGCTGGCG‑tt                                    <  1:279793/251‑1 (MQ=255)                                                                                                                                                                                                                  caTTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTTCCGTGTATTCATTGATGGCGGTGGTGTTTGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTaaaa                              <  1:649164/251‑1 (MQ=255)                                                                                                                                                                                                                         cAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGtt                       <  1:545384/251‑1 (MQ=255)                                                                                                                                                                                                                          aaGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGGGCGCTGATTg                                                                                                                                      >  1:785524/1‑144 (MQ=255)                                                                                                                                                                                                                          aaGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTg                                                                                                                                      <  2:785524/144‑1 (MQ=255)                                                                                                                                                                                                                          aaGCATTGCTGGAACACTGTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGGGCGCTGATTGGCAGCGGCGGGGTT‑AAGTTT‑TGGGCCGCCTGGCTGGGAGGGGTTGGCGGGTTTCT‑T‑CTGGGGGACTGGATTTCCTTCTGGCTGGGGCGGCGGTTTTAAAAA‑CCGTTg                      >  1:542225/1‑251 (MQ=255)                                                                                                                                                                                                                           aGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTcgctggcgctg                                                                                                                                                                                             <  1:162003/88‑1 (MQ=255)                                                                                                                                                                                                                               ttgcgtgAACACTTTTTTACCCAAGCTAGCTTTTATTCATTGATGTCTGTGGTTGTGGTGGCCTTTCTGGTGTCGCTGGCGCTGGTCGGTTTTTGTCTTCCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGTCAGCGGCGAG‑TT‑AAGTTTCTGGCACACCTGGCTGTCAGGGATTG‑TCGGTTGTC‑TGCTGGTCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCAt                   <  2:787465/243‑1 (MQ=255)                                                                                                                                                                                                                                 gCTGGAACATTTTATTACCCAATCTACCGTGTATGCATTGATGGCGTGGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGTTTTGATTCTACCCGGTGCGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGG‑TGGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCT‑TTTTAAAAAGCCGTTGCATCGCt               <  2:24023/251‑1 (MQ=255)                                                                                                                                                                                                                                       atacttgttaaCCCAATCGTCCGTGTATTCATTGATGGCGGTGGGGTTGGGGGCCTGTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AATTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCTCTGGTCAt         <  2:383256/241‑1 (MQ=255)                                                                                                                                                                                                                                        cacTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCAtt        <  1:665083/251‑1 (MQ=255)                                                                                                                                                                                                                                            ttattaCCCAATCTACGGTGTAGTCGTTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCt     <  2:361317/250‑1 (MQ=255)                                                                                                                                                                                                                                              attaCCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTG‑TCGGTTGTC‑TGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTgaa  <  1:559130/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     CAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCG‑TCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAG‑TT‑AAGTTTCTGGCACGCCTGGCTGGCAGGGATTA‑TTGGCTGCT‑TGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAA  >  NC_000913/71139‑71625

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