Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,316	T→C	intergenic (+284/‑318)	cra → / → mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,316	0	T	C	100.0%	253.2 / NA	76	intergenic (+284/‑318)	cra/mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/36);  total (40/36)

GAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGTAACAAG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAA  >  NC_000913/89079‑89559                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     gagTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGtt                                                                                                                                                                                                                                        <  1:986919/251‑1 (MQ=255)         gaatggaaataaatagccTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTAAAACCCGCCCGTGAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAAGTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCAAC‑ACAGACGCAACATG‑GACTGCTAACCgggt                                                                                                                                                                                                                                <  2:151573/234‑2 (MQ=255)            tGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGa                                                                                                                                                                                                                               >  1:719610/1‑249 (MQ=255)                 aTCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATCTCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGAAGCAACAGG‑GACTGTTAACCGGGGAAGATATg                                                                                                                                                                                                                        >  2:317366/1‑251 (MQ=255)                 aTCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATg                                                                                                                                                                                                                        >  1:441213/1‑251 (MQ=255)                     atatCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATCGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCt                                                                                                                                                                                                                    >  1:286547/1‑251 (MQ=255)                       atccTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTACTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAAATATGTCCTa                                                                                                                                                                                                                   >  2:884342/4‑250 (MQ=255)                             aaaGTAAATAAGTCTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGcc                                                                                                                                                                                                            <  2:952396/251‑1 (MQ=255)                                  aaaTAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTc                                                                                                                                                                                                        <  1:945380/250‑1 (MQ=255)                                      aaGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTc                                                                                                                                                                                                   <  1:3094/251‑1 (MQ=255)                                           tATTCCTGCCTATTATTGATACGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTTAAATTCACAATAATTAAGTGATATCGACAGCGCGATTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGCTTTAACAAGCACTTTTCTTTGCTTTTGCGAAAACCAGCTGGCATCAACCGCCAC‑ACTGACGCAACAGG‑GACTGTCAACCGGGGAAGCTATGTCCTATATGGCCGCTCGCGTCGCaa                                                                                                                                                                                               >  2:733099/1‑250 (MQ=255)                                            aTTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCaaa                                                                                                                                                                                              <  1:345035/250‑1 (MQ=255)                                              tCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACt                                                                                                                                                                                            <  1:974290/250‑1 (MQ=255)                                                      attattGATAAGTTCTTTTAACCCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGTTATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACtt                                                                                                                                                                                     <  1:583100/249‑1 (MQ=255)                                                      attattGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACttt                                                                                                                                                                                    <  1:758860/250‑1 (MQ=255)                                                      attattGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACtt                                                                                                                                                                                     <  2:654191/249‑1 (MQ=255)                                                           tGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATAtt                                                                                                                                                                               <  1:993345/250‑1 (MQ=255)                                                             ataagttcTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCCCTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGc                                                                                                                                                                            >  1:750632/8‑251 (MQ=255)                                                                atgtcTTTTAAATCCGCCTGTGGATTAACTACCCAGTTGAAATCCCCAATATATCAGTGATAACGACATCTCGTGTTTGCATTATTTAGTTCCATGCAGAGAGGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGAATCAAGCGCCAC‑ACAGCCTCAACAGG‑GACTGTTAACCGGTTAAGATAGTTCCTAAAATGCCGCTCGCGACGCCAACTGACCCTTTATAGTTGCtgt                                                                                                                                                                         <  2:727491/247‑1 (MQ=255)                                                                agttcTTTTAAACCCGCACGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCATAAAATGCCGCTCGCGTCGCACACTGGCACTTTATATTTGCtgt                                                                                                                                                                         >  1:765948/5‑251 (MQ=255)                                                                   tcTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCAAAATAATTATGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTgg                                                                                                                                                                       >  1:362319/2‑250 (MQ=255)                                                                       ttAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCACACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAAt                                                                                                                                                                  >  1:891368/1‑251 (MQ=255)                                                                             ccgccctTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑GCAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAg                                                                                                                                                            <  1:482647/251‑1 (MQ=255)                                                                                accgTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAAACGG‑GACTGTTAACAGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGCGTc                                                                                                                                                          >  2:104887/2‑250 (MQ=255)                                                                                 ccgTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCAt                                                                                                                                                        <  1:110583/251‑1 (MQ=255)                                                                                        ttaaCTCACCATCTGAAATTCACAATAATTAAGTGATATCTACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGCAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTTTGGAAAATAGGGAGTCATTATaaa                                                                                                                                                  >  2:710039/1‑250 (MQ=255)                                                                                              cACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGGTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCGGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCTCGTCGCAAACTGACACTTGATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGtga                                                                                                                                           <  1:489502/251‑1 (MQ=255)                                                                                               aCCAGCTGTAATTCATACAAATTCAGTGATATCGACAGAGCGTTGTTGCATTTTTGTGTAACATGCAGAGATGAATTGCCGAATTAACCATCACTCTTCTTTGCTTTTGCTCAAACTCGCTGGCATCAAGCGCCAC‑ACAGACGCAA‑ATGTGACTTTTAATCGTGTAAGATAAGTACTAAAATGCCGCTCGAGACGCAAACTGACACTTTATATTTGCTGTGGAAAATATTGAGTCATTTTAATACGGtgat                                                                                                                                          <  2:589684/251‑1 (MQ=255)                                                                                                 cAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGtgatga                                                                                                                                        <  1:425993/251‑1 (MQ=255)                                                                                                  aGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGAc                                                                                                                                       >  2:193809/1‑251 (MQ=255)                                                                                                     tGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGAt                                                                                                                                    >  1:158666/1‑251 (MQ=255)                                                                                                         aTTCACAATAATTAAGGGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATGTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAgg                                                                                                                                <  2:224965/251‑1 (MQ=255)                                                                                                           tCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAgggg                                                                                                                              >  1:356665/1‑250 (MQ=255)                                                                                                               aataatTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGAGATGACGATGAGGGAtttt                                                                                                                          >  1:356970/1‑251 (MQ=255)                                                                                                                   atTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTGCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCt                                                                                                                      <  1:188320/251‑1 (MQ=255)                                                                                                                               cGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTc                                                                                                           <  1:439555/250‑1 (MQ=255)                                                                                                                                gACAGCGCGTTTTTGCATTATTTTGTTAAATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCAt                                                                                                         >  2:620122/1‑251 (MQ=255)                                                                                                                                      gcgTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGTTGAGGGATTTTTTCTTACAGCTATTCATAACGtt                                                                                                   >  2:489088/1‑251 (MQ=255)                                                                                                                                        gTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTa                                                                                                  >  1:582778/1‑250 (MQ=255)                                                                                                                                                     tttGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACg                                                                                    >  1:532500/1‑251 (MQ=255)                                                                                                                                                     tttGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTCAATTTGCTGCGCACg                                                                                    >  1:878530/1‑251 (MQ=255)                                                                                                                                                      ttGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGt                                                                                   >  1:79315/1‑251 (MQ=255)                                                                                                                                                           acaagcagagattaaTTGCCGATTTAACAAGCCCTTTTCTTTGCTTTTGCGGAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCCCGCGTCGCAAACCGACACTTTAAAATTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGa                                                                               <  1:458376/237‑1 (MQ=255)                                                                                                                                                           aCATGCAGAGATGAATGTCCGATTTAACAAGCACTGTTCTTTGTTTTTGCGCAAACACGCTGGCATCAAGCGTCAC‑TCAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCTGCTCGCGTCGCAAACTGACACTGTATATTTGCTGTTTAAAATAGTGAGTCATTTTAAAACGGTTATGGCGATTAGGGATTTTTCCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGa                                                                               <  1:677588/250‑1 (MQ=255)                                                                                                                                                                cagagATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTaa                                                                          >  1:513878/5‑250 (MQ=255)                                                                                                                                                                 agagATGAATTGCCGATTTAACAAGCTCTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTaaaa                                                                        <  1:390998/248‑1 (MQ=255)                                                                                                                                                                  gagATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAAt                                                                       >  1:830923/3‑251 (MQ=255)                                                                                                                                                                         aTTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAaacaac                                                                >  1:499009/1‑251 (MQ=255)                                                                                                                                                                           tGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTAc                                                                                                                   <  1:466750/198‑1 (MQ=255)                                                                                                                                                                           tGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTAc                                                                                                                   >  2:466750/1‑198 (MQ=255)                                                                                                                                                                             ccGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAACCTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTg                                                            >  1:432902/1‑251 (MQ=255)                                                                                                                                                                             ccGATTTAACAAGCACTTTTCTTTGCTCTGGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGC‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGGCACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACACCGCTg                                                            <  1:612447/251‑1 (MQ=255)                                                                                                                                                                                   tAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAAAGCTGATAtt                                                       >  2:926320/1‑250 (MQ=255)                                                                                                                                                                                    aacaagcaCTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGTTAtt                                                       >  2:50669/1‑249 (MQ=255)                                                                                                                                                                                               ttCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCTCCAC‑ACAGACTCAACAGG‑GACTGTTTACCAGGGAAAATATGTCCTATAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGTAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCTTAAACAt                                           <  2:106294/250‑1 (MQ=255)                                                                                                                                                                                                 cttggCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTAATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCgg                                        >  2:516761/5‑251 (MQ=255)                                                                                                                                                                                                 cTTTGCTTTTGCGCAAACACGCTGGCATCAAGCGCCAC‑CCAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGGGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAACCAACACTGATATTAGCCGTAAACATCgg                                        >  2:925348/1‑251 (MQ=255)                                                                                                                                                                                                     gctgttGCGCAAACCCGCTGGCAGCAAGCGCCAC‑ACAGACGCAAAAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCCCGCGTCGCAAACTGACACTTGATAGTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGTACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGttt                                    <  2:356970/247‑1 (MQ=255)                                                                                                                                                                                                      cTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAACATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTCTTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGGTATTAGCCGTAAACATCGGGtttt                                   >  1:632376/1‑251 (MQ=255)                                                                                                                                                                                                         ttGCGCAAACCCGCTGGCATCCAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGCTATTAGCCGTAAACATCGGGtttttt                                 >  2:162913/1‑250 (MQ=255)                                                                                                                                                                                                         ttGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGtttttt                                 >  1:749672/1‑250 (MQ=255)                                                                                                                                                                                                             gcAAACCCGCTGGCATCAAGCGC‑ACTACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCACACTGACACTTTATAGTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGTTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTc                            <  1:325343/251‑1 (MQ=255)                                                                                                                                                                                                              cAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACCGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCg                           >  1:728214/1‑251 (MQ=255)                                                                                                                                                                                                              cAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCg                           >  1:169558/1‑251 (MQ=255)                                                                                                                                                                                                                aaCCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTACTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGt                         <  1:358048/251‑1 (MQ=255)                                                                                                                                                                                                                aaCCCGCTGGCATCAAGCGCCAA‑AAAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACGCTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGCTGACGTTGAGGGATTTTTTCTTACAGCTATTCATAAAGTTAATTTGCTGCGCACGTTGGACGTAAATTAAACAACGCTGATATTAGCCGTTACCATCGGGttttttctctcgg                          >  2:134169/1‑243 (MQ=255)                                                                                                                                                                                                                 aCCCGCTGGCATCAAGCGCCAC‑ACAGACGCAACCGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGt                         >  1:269217/1‑250 (MQ=255)                                                                                                                                                                                                                       tGGCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGAGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCAGCGCACGTTGGACTAAAAATAAACAACACTGATATTAGCCGTAAACATCGGGATTTTTACATCGGTATGCCtt                  >  1:617752/1‑251 (MQ=255)                                                                                                                                                                                                                        ggCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTtg                 >  2:265494/1‑251 (MQ=255)                                                                                                                                                                                                                         gCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTtgt                <  1:905929/251‑1 (MQ=255)                                                                                                                                                                                                                         gCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTtgt                <  1:805962/251‑1 (MQ=255)                                                                                                                                                                                                                         gCATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTtgt                <  1:748058/251‑1 (MQ=255)                                                                                                                                                                                                                          cATCAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTtgtg               <  1:738021/251‑1 (MQ=255)                                                                                                                                                                                                                             cAAGCGCCAC‑ACAGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGAc             <  2:229572/250‑1 (MQ=255)                                                                                                                                                                                                                                      ccccaGACGCAACAGG‑GTCTTTTAACCGGGGAAGATATGTCCTAAAAGGCCGCTCGCGTCGCAAACTGACCCTTTATATTGGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGAGTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTGTTACCTCGGTATGCCTTGTGACTGGCTTGAc    <  2:891368/248‑1 (MQ=255)                                                                                                                                                                                                                                        acaGACGCAACAGG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTTATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACAGCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACaa  <  2:286547/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     GAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCAC‑ACAGACGTAACAAG‑GACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAA  >  NC_000913/89079‑89559

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