Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,149,535	C→T	G244D (GGC→GAC)	eptC ←	phosphoethanolamine transferase EptC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,149,535	0	C	T	100.0%	146.6 / NA	50	G244D (GGC→GAC)	eptC	phosphoethanolamine transferase EptC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/24);  total (26/24)

GGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGCCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACGCGAGGCCAGGTTATCCAACGTTTTCTCAAACGGCTTGTTTTTGATAAACGTATTCATGGCGATCGGATGCAG  >  NC_000913/4149294‑4149773                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                ggTTTTTTAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCAcc                                                                                                                                                                                                                                         >  1:127018/1‑250 (MQ=255) ggTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGCTCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGCTTTCCATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTAGGATCGGTTTTATGCAGCGCATCCAGATCCGGCGTGGTTTCACGCGGATAACCGTAAAGCCTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCAcc                                                                                                                                                                                                                                         >  2:248659/1‑250 (MQ=255)   ttttATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGc                                                                                                                                                                                                                                      >  1:196257/1‑251 (MQ=255)     ttATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAc                                                                                                                                                                                                                                    >  1:241467/1‑251 (MQ=255)          aCCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCCGTTGCAGGATTTCAATGGTGTCCGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑Taa                                                                                                                                                                                                                                >  1:396814/1‑249 (MQ=255)                  gTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCACTCACCAGCAC‑TAAAGTGCGCGGt                                                                                                                                                                                                                      >  1:261792/1‑251 (MQ=255)                    ttCATCATGTTCATCAGCGACGGCTGCGTCAGATAAAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCACTGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGGCCACCCACCTCCGGCGCGGTTTCCCGCGGTTAACCGACCAG‑CTCCATGCGTCCGCGCGGGGTCGCCTGGTCATCCACCAGCAC‑TAAAGTGCGCGGTTc                                                                                                                                                                                                                    >  1:381943/1‑251 (MQ=255)                      catcatGTTCATCAGCGACGGCTGCGTCAGATACAGGTCCGGGTTCTTTTCATTGTCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCg                                                                                                                                                                                                                   <  1:383648/250‑1 (MQ=255)                               catcaGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTGTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGAt                                                                                                                                                                                                          <  2:347600/250‑1 (MQ=255)                                    gCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCt                                                                                                                                                                                                    >  1:201754/1‑251 (MQ=255)                                     cGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCtt                                                                                                                                                                                                   >  1:8186/1‑251 (MQ=255)                                     cGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCtt                                                                                                                                                                                                   <  1:344219/251‑1 (MQ=255)                                      gacggtttcGTCAGAGACAGATCCAGTTTCTTTTCATTGGCAAAGGTCAGCGACTTTTTCAGGATTTCAATGGGGGACTGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCTCCAGCCC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCttt                                                                                                                                                                                                  <  2:320756/243‑1 (MQ=255)                                         ggCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCGCCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGa                                                                                                                                                                                                >  1:103891/1‑250 (MQ=255)                                          gCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGaa                                                                                                                                                                                               >  1:297734/1‑250 (MQ=255)                                              cGTCAGCTACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTTTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGAACGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTA‑ATGCGTCCGCGCTGGGTCGACTCGTCAATAACCAGCAC‑TAAAGTGCGCGGTTCGTTAACCGATTCATCTTTGAAAtt                                                                                                                                                                                            >  2:263015/1‑249 (MQ=255)                                                     tACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGCTTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTgg                                                                                                                                                                                   >  1:221487/1‑251 (MQ=255)                                                            ccGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAg                                                                                                                                                                            <  1:136425/251‑1 (MQ=255)                                                             cGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTCCGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAgg                                                                                                                                                                           <  1:326220/251‑1 (MQ=255)                                                                         aTTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTttcct                                                                                                                                                               >  2:317221/1‑249 (MQ=255)                                                                           tGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTTGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTttcatt                                                                                                                                                              <  1:309194/250‑1 (MQ=255)                                                                                       gcgcCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATGCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTGCATTCAGTAACTTTGtt                                                                                                                                                 <  2:282668/251‑1 (MQ=255)                                                                                                     aTTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGCCTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTCTTTTCATTCAGTAACTTTGTTAGCGAGGTTAGTTg                                                                                                                                   >  2:20926/1‑251 (MQ=255)                                                                                                         cAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCTGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGa                                                                                                                               <  1:224039/251‑1 (MQ=255)                                                                                                                   cGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCGCCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGat                                                                                                                     <  2:103891/251‑1 (MQ=255)                                                                                                                         agaCGTAAATACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTACCTTTGTTAGCGAGTTTAGTTGCTGACGGTACTGATAATAGc                                                                                                               >  2:367606/1‑251 (MQ=255)                                                                                                                             gTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTTCCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGt                                                                                                           <  2:74643/251‑1 (MQ=255)                                                                                                                             gTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTCTATGCAGCGCATCCAGCTCCAGCGTGGTTTCACGCGAATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGGAACTTTGTTAGCGGGTTTAGTTGCTGAAGATACTGATAATAGCCGGt                                                                                                           >  2:45558/1‑251 (MQ=255)                                                                                                                              aacagaCGTTATTAAACACGGTCTGATTGGGATCGGTTTTATGCATCGCATGCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGGTGGTTCGACTCGTCAATCGCCTGTAC‑TAAAGTGCGCGGTGCGTTACCCGATTCATTTTTGAAATTATCCAGGGGCGGCAAGGCAGTATTTTCATTCAGTAACGATGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTc                                                                                                          <  2:102926/246‑1 (MQ=255)                                                                                                                               aaCTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCa                                                                                                         <  1:6523/251‑1 (MQ=255)                                                                                                                                  tACGTTATTCAAGACGGTCAGGTTCGGGTAGGTTTTCTGCAGCGCATCCAGATCCGGCGTGGATTCACGAGGTCAACCGTAAAGCCTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAACCACCAG‑ACATACAGTGCGCGGTTCGTTACCCGATTCATCTTTGCAATTAGCCCGTGGCGGCCAGGCATTATTTTCATTCAGTAACTTTGTTAGAGAGTTTAGTTACTGACGATAATGATAATAGCCGGACAgt                                                                                                       >  2:259294/1‑249 (MQ=255)                                                                                                                                      ttattaAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTg                                                                                                  >  1:108821/1‑251 (MQ=255)                                                                                                                                                cGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGAGCGGGTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCCTTCAGTAACTTTGTTAGCGAGTTTAGTTGATGACGATACTGATAATAGCCGGTCAGGAATTGCCACgtggcg                                                                                        >  2:308561/1‑246 (MQ=255)                                                                                                                                                          cGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGAATC‑ATGCGCCCGCGCTGGGCCGACTCGTCAATCACCAGCAC‑TAACGTGCGCGGTTCGTTCACCGCTTCCTCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCATTAACTTTGTTAGCGAGTTTAGTTGCTGCAGATACTGATAATATGCGGTAAGGAATTGCCACGGCGCGGCAGGGTcc                                                                               >  2:357156/1‑250 (MQ=255)                                                                                                                                                              gcggtttaagccAGCTCATCCCGCTCCGGCGTGGTTTCACGCGGAGAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCCC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTCGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGGTTCTGACCATACTGATAATAGCCGTTCAGTAATTGCCACGGTGCGGCAGGCTCCATAc                                                                           <  1:413106/239‑1 (MQ=255)                                                                                                                                                                 gTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACGCGa                                                                       <  1:247902/251‑1 (MQ=255)                                                                                                                                                                  ttttATGCAGCGCATCCCGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGAGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAACTTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAAATTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCAAGGTGCGGCAGGCTCCATAagcgcg                                                                      >  2:304971/1‑249 (MQ=255)                                                                                                                                                                     tATGCAGCGCATCCAGCTCCGGCGTTGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGTTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGTCAAGGCATTATTTTCAGTCAGTAACTTTGTTAGCGATTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAAGCGCCACGGTGCGGCAGGCTCCATACGCGAGGc                                                                    <  1:259294/250‑1 (MQ=255)                                                                                                                                                                     tATGCAGCGCATCCAGCTCCGGCGTGGTTTAACGCGGATAACAGTAAAGACAC‑ATGCGTCCGCGCCGGGCAGACTCGTCAATACCCAGCAA‑TACTGTGCGCGCATCGCTACCCGATTCATCTTTGAAATTACCAAGTGGCGGCAAAGCATTCTTTTCATTCAGTACCCTTGTTAGAGCGTTTAGTTGCTAACGATACTGATCATCACCGGTCAGGAATTGCCACGGTGCGGCAGGCCCCAAAAGCGAGGcc                                                                   >  1:415013/1‑251 (MQ=255)                                                                                                                                                                            cgcATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACGCGAGGCCAGGTTAt                                                            >  1:437035/1‑251 (MQ=255)                                                                                                                                                                                 ccAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACGCGAGGCCAGGTTATCCaa                                                        <  1:99004/250‑1 (MQ=255)                                                                                                                                                                                   aGCTCCGGCGTGGTTTCACGGGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGATTCGTCAATCACCAGCAC‑TAAAGTGCGCTGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTGGCTGACGATACTGATAATAGCCGGGCAGGAATTGCAACGGTGCGGCAGGCTCCATACGCGAGGCCAGGTTAGCCAACGt                                                     <  2:112612/251‑1 (MQ=255)                                                                                                                                                                                   aGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGCCTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACTCGAGGCCAGGTTATCCAACGt                                                     <  2:372814/251‑1 (MQ=255)                                                                                                                                                                                                   cACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCGCCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACGCGAGGCCAGGTTATCCAACGTTTTCTCAAACGGCTTg                                     <  1:374995/251‑1 (MQ=255)                                                                                                                                                                                                                 tccagtcTC‑ATGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTCCGTTACCCGATTCATCTTTGAAATAAGCCAGTGGCGGCAAGGCCTTATTTTCCTTCAGGAACTTTGTTAGCGTGTTTAGTTGCTGACGATCTTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGGTACATACGCGAAGGCAGGCTATCCAACGGTTTTTCCAACGGCTTGATTTTGAGAAAc                         >  2:377907/7‑249 (MQ=255)                                                                                                                                                                                                                           aTGCGTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTCTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGGTCCATACGCGAGGCCAGGTTATCCAACGTGTTCGCAAACGGCTTGTTTTTGATAAACGTATTCATGGc              <  1:172467/251‑1 (MQ=255)                                                                                                                                                                                                                               gtccgagatGGGTCGACTCATCCATCACCAGAAC‑TAAAGTGCGGGGTTCGTTACCCGAGTCATCTTTGAAATTAGCCAGTGGCGGCCAGGCCTTATTTTCATTCAGAAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGTAAATAGCCGGTTAGGAATTGCCACGGTGAGGCAGGGTCCAAACGAGAGGCCAGGTTAACCAACGTTTTCTCAAACGGATTGTTTTTGCTACACGTATTCATGGCggac          >  2:57970/9‑248 (MQ=255)                                                                                                                                                                                                                               gTCCGCGCTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCAGGCAGGAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACGCGAGGCCAGGTTATCCAACGTTTTCTCAAACGGCTTGTTGGTGATAAACGTATTCATGGCGATc          <  2:221487/251‑1 (MQ=255)                                                                                                                                                                                                                                    cgcTGGGTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACGCGAGGCCAGGTTATCCAACGTTTTCTCAAACGGCTTGTTTTTGAGAAACGTATTCATGGCGATCGGATg     <  1:171501/251‑1 (MQ=255)                                                                                                                                                                                                                                        gggTCGACTCGTCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACGCGAGGCCAGGTTATCCAACGTTTTCTCAAACGGCTTGTTTTTGATAAACGTATTCATGGCGATCGGATGCAg  <  1:248659/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                GGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTC‑ATGCGTCCGCGCTGGGTCGACTCGCCAATCACCAGCAC‑TAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTGATAATAGCCGGTCAGGAATTGCCACGGTGCGGCAGGCTCCATACGCGAGGCCAGGTTATCCAACGTTTTCTCAAACGGCTTGTTTTTGATAAACGTATTCATGGCGATCGGATGCAG  >  NC_000913/4149294‑4149773

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