Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,630,400	A→G	N41N (AAT→AAC)	ettA ←	energy‑dependent translational throttle protein EttA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,630,400	0	A	G	100.0%	164.5 / NA	56	N41N (AAT→AAC)	ettA	energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/22);  total (34/22)

AGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACG  >  NC_000913/4630154‑4630646                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       cgcggcaagtttgaaaaaaTCGGCATCCGGATCGGCGTACAGCGCATCCACTTCCTCCAGGCGTTCCGGGGCGTTAACCACTTCAGAAACCTCTTCTTCAAGGGACTCACGCGCGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTGGATGTCGGG‑CTGCGGACGCGCTCCACCTTCGATGT‑CTTTAACAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCa                                                                                                                                                                                                                                                    <  2:188267/236‑1 (MQ=255)   cgtccAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCCGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCACTGCCCGCCATAATGCTCAGCAGGGTAGACTTACCCGCGCCGTTccg                                                                                                                                                                                                                                                   >  2:161929/4‑248 (MQ=255)            tgccaaaaccgccaccgggaccGGCGTACAGCGCATACACTCCACCCAGGCGTTTCAGGGCGTAAACCACTTCAGAAACCGCTTCTTCAATGGCCTCCCGCACGGTGTTTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGCCGCGCTTCACCTGCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCTCCGTTCAGACCCAGGa                                                                                                                                                                                                                                           <  1:505387/228‑1 (MQ=255)            tGTCACAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGACCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCCGACCCAGGac                                                                                                                                                                                                                                          >  2:625231/1‑250 (MQ=255)                     cGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAAttttt                                                                                                                                                                                                                                >  1:559731/1‑251 (MQ=255)                          cccGGATCGGCGTCCAGCGCCTACCTTTCCCCCAGGCGTTTCAGGGCGTTAACCACTTCAGAACCCTCTTCTTCAATGGACCCCCCCCCGGTGTGTTCCGGTTCCAGCTGCGGTTCCTCCGCCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑TTTTATCACTGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCTCGCCGTTCAGACCCAGGACACCAATTTTTGcccc                                                                                                                                                                                                                           <  1:178694/250‑1 (MQ=255)                               aTCGGCGTACAGCGCATACACTTCATCCATGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCATGCGGCAGGTAACCAATATTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCCATGACCGCCATAAAGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGAACCAGGGCACCAATTTTTgagatcggc                                                                                                                                                                                                                       >  2:102333/1‑242 (MQ=255)                                    cGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGgaagaa                                                                                                                                                                                                                  >  2:476388/1‑250 (MQ=255)                                    cGTACAGCGCATACACTTCATACAGACGGTTCAGGGCGTTAACCAATTCAATAACCGCTTCTTCATTGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTTTGG‑ATGCGGACGCGCTTCACCTTCGCAGT‑CTTTATCAATGACCGCCATAATGCGCCGCAGGGGAGACTTACCCGCGCCGTTCATACACAGGACAACACTTTTTGTCCCAGGgcagaa                                                                                                                                                                                                                  >  2:42303/1‑250 (MQ=255)                                      tACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGCTGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGCCACCAATTTTTGCCCCAGGGAAGCAACt                                                                                                                                                                                                               >  2:164643/1‑251 (MQ=255)                                              aTACACTTCATCCAGGCGTTTCAGGGCGTTCACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCagag                                                                                                                                                                                                          >  2:86699/1‑248 (MQ=255)                                                      cATCCAGGCTTTTCAGGGCGTTACCCACTTCAAAAACCGCTCCTTAAATGGACTCCCGCACGGTGTGTTCCGGGTCCAGCGGCGGTCCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGCCGCGCTTCCCCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTc                                                                                                                                                                                               <  1:457412/251‑1 (MQ=255)                                                          cAGGCGTATCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGAATTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCaaa                                                                                                                                                                                            >  2:572452/1‑250 (MQ=255)                                                             gCGTTTCAGGGCGTTAACCCCTTCAGAAACCGCTTCTTAAATGGACTCACGCGCGGTGTGTGCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATAAAGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAatat                                                                                                                                                                                        <  1:473122/251‑1 (MQ=255)                                                                    aGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGttt                                                                                                                                                                                 >  2:113107/1‑251 (MQ=255)                                                                     gggCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGtttt                                                                                                                                                                                >  2:213968/1‑250 (MQ=255)                                                                       gCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGCTTTTTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATACCCCCCATAATGCGCAGCAGGGGAGACTTGCCCGCGCCGTTCAGACCCCATACACCAATTTTTGCCCCAGGGAAGAAAATCAGAGAGATGTTTTTCAAAATATGGAGTTTcgg                                                                                                                                                                              >  2:336752/1‑251 (MQ=255)                                                                         gTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑ATGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCACGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTcggcg                                                                                                                                                                            >  2:96131/1‑251 (MQ=255)                                                                              cccctgcagAAACCGTTTTTTCAATGTACTCCCGCACGGTGGGTTCCGGGTTCAGCGCCGGTCCCGGCGGCAGGTAACCAATTTTGATGTCTGG‑CGGCGGACGCGCTCCCCCTTCGATGT‑CTTTACCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGaaca                                                                                                                                                                       <  2:488999/243‑1 (MQ=255)                                                                                aCTTCAGAAACCGTTTCTTCAAGGGCCCCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCCCCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGaacaac                                                                                                                                                                     <  1:493475/251‑1 (MQ=255)                                                                                   tcagAAACCGCTTCTTCAATGGACTCACGCACGGCGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTCACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACCCCAATTTTTGCCCCCGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACttt                                                                                                                                                                  >  1:231638/4‑251 (MQ=255)                                                                                     agAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAAAAACTTTGc                                                                                                                                                                >  1:399653/2‑251 (MQ=255)                                                                                          ccgttttttcAATGGCCTCACGCTCGTTGTGTTCGGGGTTCGCCGGCGGTTCCGGCGGCAGTTAACCACT‑TGGGTTGCTGGGCGGCGGACGCGCTTCACCTGCGTTGT‑CTTTACAAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGTCACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATTACGTTTCGGCGGAACAACTTTGCCGaca                                                                                                                                                           <  2:429461/244‑1 (MQ=255)                                                                                           cGCTTCTTATATGGACTCACTCTGGGTGTGTTCCGGGTTCAGCTGCGGTTCCCGCGGCAGGTAACCATGCTTTATGTCTGG‑CTGCGCACGCGCTCCCCCCTCG‑TGTTCTTTCTAAATGTCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGAGGTTTTGCAAAAGATGACGTTTCGGCGGAACAACTTTGCCGacac                                                                                                                                                          <  2:248148/251‑1 (MQ=255)                                                                                               tcttcaggggccccccgcacgGTGTGTTCGGGTTTCAGCTGCGTGTCCTGCGGCAGGTAACCAATCTTGATGTCGGG‑CTGCGGACGCGCTCCACCTTCGATTT‑CTTTACCAGGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGAt                                                                                                                                                       <  1:86699/235‑1 (MQ=255)                                                                                                  tcAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCAt                                                                                                                                                   >  1:23611/1‑251 (MQ=255)                                                                                                                           cggggTTCAGGTGTGGTCCCTGCGGCAGGTAACAAATCTTGAGGTCGGG‑CGGGGGACGCGCTTCTCCTTCGATGT‑TTTTAAAAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTg                                                                                                                          <  1:178578/249‑1 (MQ=255)                                                                                                                            cGGGTTCAGCTTCGGTTCCTGCGGCAGGTAACCAATCTTGGTGTCTGG‑CTGCGGACTCGCTTCCCCTTCGATGT‑CTTTCTCAATTCCCGCCATAATGCGCAGCAGGGTAGACTTACCCCCGCCGTTCAGACCCAGGCCACCAATTTTTGCCCCAGGGAAGAAACTCCGAGAGATGTTTTTCAAAATATACCGTTTCGGCCGAACACCTTTGCCGACACGATGCATGGTATAAACGAATTTAGCCACGtttg                                                                                                                         >  2:42/1‑249 (MQ=255)                                                                                                                             gggTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACg                                                                                                                                         >  1:548369/1‑234 (MQ=255)                                                                                                                                ttCAGTTGGGGTTCCTGCGGCGTGTAACCAATCTTTATGTCTGG‑CTGCGGACGCGCTTCACCTGCGGTGT‑CTTTCTCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACt                                                                                                                      <  1:249573/250‑1 (MQ=255)                                                                                                                                    gCTGCGGTTCCGGCGGAAGGTAACCAATCTGGATGTCTGG‑CGGCGGCCGCGCTCCCCCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGcc                                                                                                                 <  1:353146/251‑1 (MQ=255)                                                                                                                                       gCGGTTCCGGCGGCAGGTACCCAATCTTGTTGTCGGG‑CGGCGGCCGCGCTCCCCCTGCGATGT‑TTTTACCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCct                                                                                                                <  1:585617/249‑1 (MQ=255)                                                                                                                                        cGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCCCCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCCTGGTATAAACGTATTGAGCCACGTTGGACTTCGCctctt                                                                                                             >  1:822/1‑250 (MQ=255)                                                                                                                                            tCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGtt                                                                                                         >  1:321708/1‑251 (MQ=255)                                                                                                                                                gCGGCAGGTACCCAATCTTGATGTCTGG‑CTGCGGCCGCGCTTCCCCTGCTGTGT‑CTTAATAAAGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTa                                                                                                       <  1:444012/249‑1 (MQ=255)                                                                                                                                                       gtAACCAATCTTGATGTCTGG‑CTGCGGCCGCGCTTCACCTTCGATGT‑CTTTATCAATACCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCCGAGAGATGTTTTTCAAAATATGACGATTCGGCGGAACAAATTTGCCGACACGAAGCATGGTATTAACGAATTGAGCCACGTTGGCCGTCGCCTCTATGTTTATCGTGAtg                                                                                               >  2:378275/2‑249 (MQ=255)                                                                                                                                                                 ttGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTc                                                                                     <  1:188495/250‑1 (MQ=255)                                                                                                                                                                  tGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                                                                                                                  >  2:427807/1‑140 (MQ=255)                                                                                                                                                                           ggcggCGGACCCGCTTCCCCTGCGTTGT‑TTTTCACAAGGCCCGCCATAATGCGCAGCAGGGGAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAg                                                                           <  1:131637/247‑1 (MQ=255)                                                                                                                                                                           g‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAAtt                                                                                                                                                                                                                                   >  1:168207/1‑98 (MQ=255)                                                                                                                                                                                  gACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGcc                                                                    >  1:100247/1‑251 (MQ=255)                                                                                                                                                                                          cctcctccGATGT‑CTTAACAAATCCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACTCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGGTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTccc                                                            <  1:425866/244‑1 (MQ=255)                                                                                                                                                                                           cgacttgtgatacgttaaccAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATTTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTcc                                                             <  2:3616/231‑1 (MQ=255)                                                                                                                                                                                           ccccTTCGATGT‑CTTTACCACTGCCCGCCATAATGCGCAGCAGGGTAGCCTTACCCTCGCCGTACAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGATATGTTTTTCAAACTACGACGTTTCGGCGGAACACCTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTTATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTccc                                                            <  1:494109/248‑1 (MQ=255)                                                                                                                                                                                                  gATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAAt                                                    >  1:609623/1‑251 (MQ=255)                                                                                                                                                                                                    tGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATcc                                                  >  2:495373/1‑251 (MQ=255)                                                                                                                                                                                                         tttATCAATGACCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGCTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCCCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAgt                                              >  2:157987/1‑250 (MQ=255)                                                                                                                                                                                                          ttatccctGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAgtt                                             >  1:483778/8‑249 (MQ=255)                                                                                                                                                                                                             tCAATGCCCGCCATAATGCACAGCCGGGTAGCCTTAACCGCGCCGTTCAGACCCAGGACACCACTTTTTGCCCCAGGGAAGAAACTCAGCGAGATGTTTTTCAACATATGACGTTTCGGCGGAACAACTTTGACGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGCAGTGTAGCCTTTTTCCCCGCCTACTCCCAgttggt                                          >  2:311568/1‑246 (MQ=255)                                                                                                                                                                                                                  ccccGCCAGAATGCGCATCAGGGTAGACTTACCCGCGCCGTTCAGGCCCAGGACACCAATTTTGGCCCCATGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACGTTGCCGACACGATTCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTTGATGTGTATCGTGAGGATGAGTTTTCAAAGGCGAAGTTTAGCCTTGGTCCCCGCCTAATCCCAgttggccgatc                                     <  2:623969/250‑11 (MQ=255)                                                                                                                                                                                                                    ccGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGCAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCAc                                   >  1:33838/1‑251 (MQ=255)                                                                                                                                                                                                                                        gtagACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGCTGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACATTCGCAgtaaaagtaaaa                >  1:256425/4‑250 (MQ=255)                                                                                                                                                                                                                                           gACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGGGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGTTGCATGGTATAACCGAATTGAGCCCCGTTGGACTTCGCCTCTATGTTTATCGTGATGATGGGTTTTCAAAGGCGAGGTTTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAgt            >  2:340562/1‑251 (MQ=255)                                                                                                                                                                                                                                               tACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAgtgt          <  1:60654/249‑1 (MQ=255)                                                                                                                                                                                                                                                aCCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCg       >  1:70702/1‑251 (MQ=255)                                                                                                                                                                                                                                                     cgcCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTAAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTACCCGCCTAATCCCAGTTGGTCGATCACCTTAGCAGTAAAAGTAAAAAAGTGTCCGTAACg  >  1:403109/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       AGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACG  >  NC_000913/4630154‑4630646

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