Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,630,385	G→A	S46S (TCC→TCT)	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,385	0	G	A	100.0%	155.7 / NA	53	S46S (TCC→TCT)	ettA	energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/22);  total (31/22)

TGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCAC  >  NC_000913/4630149‑4630613                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     tGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGcc                                                                                                                                                                                                                        >  1:321023/1‑251 (MQ=255)  gTTCGGCGGCCAGCTTGACAAAATCGGCCTCCGGAGCGCCTTACAGCCCATCCCCTTCACCCCGGCTTTTCAGGGCTTTAACCTCTCCAAAACCCGCTTCTTCAATGGACTCCCGCACGTTGTTTTCCGGGTTCATCTGCGGTTCCTGCGGCGGGTACCCATTCTTGATGTCTGG‑CTGCGGACGCGTTTCACCTCCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGcc                                                                                                                                                                                                                        <  1:24686/250‑1 (MQ=255)      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)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     TGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTTCACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCAC  >  NC_000913/4630149‑4630613

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