Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,564,294	C→T	K480K (AAG→AAA)	dosP ←	oxygen‑sensing c‑di‑GMP phosphodiesterase DosP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,564,294	0	C	T	100.0%	188.0 / NA	57	K480K (AAG→AAA)	dosP	oxygen‑sensing c‑di‑GMP phosphodiesterase DosP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/28);  total (29/28)

TAAACCAGTTTCAGTTGGTTATTG‑CTAATCGCTTCTTTCAGCGCTGCGCCTAAAACCAAACGCTCTTTT‑ACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGCTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTG  >  NC_000913/1564052‑1564537                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    aaaaCCAGTTTCAGTTGGTTATTG‑TTAATCGCTTCTTTCAGCGCGGCGACTAAAACCAAACGCTCTTTT‑CCCATTTCGTTCAGCGCCGGGCTGAAAAATTGCCGGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGAc                                                                                                                                                                                                                                             <  2:169172/250‑1 (MQ=255)            aaGTTGGTT‑TTGTCTAATCCCTTCTTTCAGCCCTCCCCCTAAAACCAAAC‑CTCTTTTTACCATTCCTTTCATCCCCGGCCGGAAGAACCGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTa                                                                                                                                                                                                                                  <  2:502828/250‑1 (MQ=255)                 ggTTATTG‑CTAATCGCTTCTTTCAGCGCTGCGCCTAAAACCAAACGCTCTTTT‑ACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGCAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTc                                                                                                                                                                                                                              >  2:475465/1‑250 (MQ=255)                     aTTG‑CTAATCGCTTCTTTCAGCGCTGCGCCTAAAACCAAACGCTCTTTT‑ACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCgg                                                                                                                                                                                                                         >  1:69371/1‑251 (MQ=255)                                   tcttccaccccTGCGCCTAAACCCAAACCCTCTTTT‑CCCATTTCGTCCATCCCCGGTCTGAAAACCTCCCAGCCTTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTa                                                                                                                                                                                                              <  2:28033/239‑1 (MQ=255)                                       tcagggcggcgcTTAAAACAAAACTGTCTTTC‑ACCATTTCGTTCATCGCCGGGGGGAAAAACGCCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGtt                                                                                                                                                                                                         <  2:217854/242‑1 (MQ=255)                                             cTGCGCCTAAAACCAAACGCTCTTTT‑ACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGAcg                                                                                                                                                                                                  >  2:148464/1‑251 (MQ=255)                                               gcgcCAAAACCCAAACGCTCTTTT‑ACCTTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGTCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCTCATTCCGTATCTCATCGGCGATTTGGGTAATGTTACTGAcgtc                                                                                                                                                                                                <  1:558240/251‑1 (MQ=255)                                                    tAAAACCAAACGCTCTTTT‑ACCATTTCGTTCACCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGtt                                                                                                                                                                                             <  1:360766/249‑1 (MQ=255)                                                                ctctTTT‑ACCATTTCGTTCATCGCCTGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTACGAATATAATCCCTTGCATTGTGAGCAGTGGAGAGCACGTAATCGCGGTTTTTACCCTGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAATGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACg                                                                                                                                                                                >  2:24559/1‑250 (MQ=255)                                                                ctctTTT‑ACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACgg                                                                                                                                                                               >  1:138160/1‑250 (MQ=255)                                                                 tctTTT‑ACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACg                                                                                                                                                                                >  1:140010/1‑249 (MQ=255)                                                                          cattcctttCAACGCGGGGCTGAAAAACGCCCGGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTGCCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCTTAGCTCATCGGCGATGTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACt                                                                                                                                                                      <  2:237958/244‑1 (MQ=255)                                                                                     ccgctggggtGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAAACCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATTCTCAAGGTAAGCGGGAAGGGTTTATCGTTAATCATTATCTGTTTGCTGACCACATTCCGTAGCTCATCGGCTATTTGTGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTc                                                                                                                                                            <  1:416946/241‑1 (MQ=255)                                                                                       gCCGGGCTGAAGAACTGCCAGCCGTTACAGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAAATGCGTACCtt                                                                                                                                                             >  1:547505/1‑247 (MQ=255)                                                                                            gCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGAc                                                                                                                                                    >  1:397171/1‑251 (MQ=255)                                                                                              tGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAGGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGGTTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAAATGCGTAACTTCGATAAGAccg                                                                                                                                                  >  1:93612/1‑249 (MQ=255)                                                                                              tGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGGCAAACTGCGTACCTTCGATACGAc                                                                                                                                                    >  2:81806/1‑249 (MQ=255)                                                                                              tGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGAc                                                                                                                                                    >  2:90327/1‑249 (MQ=255)                                                                                                aagaccggccatcctttACCGCCATTCTTGCGAATATCATCCATTTCATTGTGTGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATTCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTGTGCTGACCACATACCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAAGGCAAACTGCGTACGTTCGATACGACaga                                                                                                                                                 <  2:222309/235‑1 (MQ=255)                                                                                                       gCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTg                                                                                                                                          >  2:170563/1‑250 (MQ=255)                                                                                                         cAGCCGTTACCGCAATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAACCAGGTCGTAGCTGATGCCAATACTAAAGGAAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGAAAAACTGCGTACCTTCGATACGACAGAGATACTNaac                                                                                                                                       >  2:608353/1‑249 (MQ=255)                                                                                                         aaGCCGTTACCGCCATTCTTGCTAATATAATCCATGGCATTGCGAGCAGTGGAGAGTAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTTCTGATGTCGTTTTCGAGGCTCACGAGGACAAACTGCGGACCTTCGATACGACAGAGATCCTGATc                                                                                                                                       <  2:167644/250‑1 (MQ=255)                                                                                                           gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATAGGCGATTTGGGTAATGTTACTGACGACGTTTTCTAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCg                                                                                                                                     >  1:497725/1‑251 (MQ=255)                                                                                                             cGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACAGCGTACCTTCGATCCGACAGAGATACTGATCCGGt                                                                                                                                   >  2:369802/1‑251 (MQ=255)                                                                                                                   cGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGTTTTTTACCCAGGTCGTAGCTGATGCCAATGCTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATGCGGATCCGGTTTGAGt                                                                                                                             <  2:616349/251‑1 (MQ=255)                                                                                                                    gCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGtt                                                                                                                            >  1:652146/1‑251 (MQ=255)                                                                                                                       aTTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGttttt                                                                                                                         >  1:636641/1‑251 (MQ=255)                                                                                                                         tcctGCGAATATAATCCCTTGCATTGTGAGCAGTGGAGAGCCCGTAATCGCGGTTTTTACCCAGGTAGTAGCTGATGCCACTACTCAAGGAAAGCGGGAAGGGTTTATCGTCAATCAGTATGGGTTTGCTGAACACATTCCGAAGCTCATCGGCGATTTGGGTAATGTTACTGACGACGTTTTCGAGGCGCACGAGGAAAAACTGCGTACCTTCGGTACGACAGCGATACTGCTACGGTTTGAGtttgtc                                                                                                                        >  2:610261/4‑247 (MQ=255)                                                                                                                              cGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGaaa                                                                                                                  <  1:14382/251‑1 (MQ=255)                                                                                                                              cGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGaaa                                                                                                                  <  1:442600/251‑1 (MQ=255)                                                                                                                                aaTATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATGTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGc                                                                                                                <  2:238663/251‑1 (MQ=255)                                                                                                                                   ataATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGa                                                                                                              <  2:220397/250‑1 (MQ=255)                                                                                                                                                             agagCAAGTAATCGCGGCTTTTACACAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGAGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTAACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATACGGTGTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCg                                                                                   >  2:423682/1‑251 (MQ=255)                                                                                                                                                                  aaGTAATCGCGGCTTTTACCCAGGTCGTAGATGATGCCAATACTCAAGGTAAGCGGGCAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGCAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGCTTTTCACGAAAACGCTTGCCCACTTCCAGCAATGCTTGATCGcccc                                                                               >  2:636609/1‑250 (MQ=255)                                                                                                                                                                         cgcagTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGTCCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAg                                                                       <  1:24834/247‑1 (MQ=255)                                                                                                                                                                          gcgGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTATTGACGTCGTTTTCGAGGCCCACGGGGAAAAACGCCGTACCTTCGATACGCAAGAGAAACTGTTCGGGTTTGAGTTTTCAACGAAAGCGTTTGACCACTTCCACCAATGCCTGATCGCCCCACACATacc                                                                      >  2:575876/1‑249 (MQ=255)                                                                                                                                                                            ggTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGTCCACTTCCAGCAATGCCTGATCGGCCCACGCATAGccc                                                                    >  2:424895/1‑250 (MQ=255)                                                                                                                                                                               ttttACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAg                                                                  <  1:81806/250‑1 (MQ=255)                                                                                                                                                                                   aCCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCt                                                                                        >  2:200823/1‑224 (MQ=255)                                                                                                                                                                                   aCCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCt                                                                                        <  1:200823/224‑1 (MQ=255)                                                                                                                                                                                       aGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACg                                                                                                                                                      <  1:265565/158‑1 (MQ=255)                                                                                                                                                                                       aGGTCGTAGCTGATGCCAAAACTCAAGGTAAGCGGGAAGGGTTTATCGTCACTCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGAACCTTCGATACg                                                                                                                                                      >  2:265565/1‑158 (MQ=255)                                                                                                                                                                                         gTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCCAGGCTatcaatc                                                        >  2:294448/1‑250 (MQ=255)                                                                                                                                                                                                 gggtGCCCATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCTGTTTGCTGACCACATTGCGTAGCTCATCTGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAAGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCAGTTTGAGTTTTTCCCGCAAGCGATTGATCACTCCCAGCAATGCCTGATGCGCCCACGCATAGCCAAGGTTATCAATCGCATCCTGa                                               <  2:544372/248‑1 (MQ=255)                                                                                                                                                                                                  gATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGaa                                              <  2:448788/251‑1 (MQ=255)                                                                                                                                                                                                   aTGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAat                                             <  2:164824/251‑1 (MQ=255)                                                                                                                                                                                                     gCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAata                                            <  1:440216/250‑1 (MQ=255)                                                                                                                                                                                                      ccAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATg                                          >  1:158512/1‑251 (MQ=255)                                                                                                                                                                                                       cAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATgg                                         >  2:329807/1‑251 (MQ=255)                                                                                                                                                                                                              cAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAacac                                  <  1:165477/251‑1 (MQ=255)                                                                                                                                                                                                                     aGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATgag                           <  1:445790/251‑1 (MQ=255)                                                                                                                                                                                                                     aGCGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATgag                           <  1:20811/251‑1 (MQ=255)                                                                                                                                                                                                                       cGGGAAGGGTTTATCGTCAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGAt                         <  1:525951/251‑1 (MQ=255)                                                                                                                                                                                                                                        cAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGAc         <  1:646838/250‑1 (MQ=255)                                                                                                                                                                                                                                        cAATCATTATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACg        >  2:571628/1‑251 (MQ=255)                                                                                                                                                                                                                                              ttATCGGTTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGACCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTg  >  1:628765/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    TAAACCAGTTTCAGTTGGTTATTG‑CTAATCGCTTCTTTCAGCGCTGCGCCTAAAACCAAACGCTCTTTT‑ACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAAGGGTTTATCGTCAATCATTATCGGCTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGGGTAATGTTACTGACGTCGTTTTCGAGGCTCACGAGGACAAACTGCGTACCTTCGATACGACAGAGATACTGATCCGGTTTGAGTTTTTCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTG  >  NC_000913/1564052‑1564537

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