Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,193,996	T→G	D70A (GAT→GCT)	thiE ←	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,193,996	0	T	G	100.0%	163.1 / NA	56	D70A (GAT→GCT)	thiE	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (28/28);  total (28/28)

ATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCT  >  NC_000913/4193749‑4194240                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     aTGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGt                                                                                                                                                                                                                                                   >  1:206997/1‑251 (MQ=255)  tGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGtt                                                                                                                                                                                                                                                  >  1:73276/1‑251 (MQ=255)   gCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGTTTTGCGTCGGGAACACAGGTCCCAGCGGAATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTg                                                                                                                                                                                                                                                 <  2:204370/251‑1 (MQ=255)     cGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGTAACACATGTCCCAGCGCGATAGAAGAGGGGCGTGCTGCCAGCGCGACTTCGATGTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGTCGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGAt                                                                                                                                                                                                                                               <  2:178988/251‑1 (MQ=255)      gggCCAGCTGTTCCAGCCACTGCGGTGCAGAAGGTATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGCCACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCACATGGCCGC‑CATCCGCCTGATGCTTGATCACCAGCCGCCAGTACGCGTTGAt                                                                                                                                                                                                                                               >  2:215252/3‑250 (MQ=255)                 tCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGACCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATcgcg                                                                                                                                                                                                                                    <  1:140700/250‑1 (MQ=255)                       cccTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTACCAGCGCGATATAAGAGGAGCGCCCTGCCAACGCGAAGTCGATTTCCATATCGTCATGTGTCGACACGCCAAGCCGCAGGCCTGCCCCGCGGCTGCCATTGAGATCGGTG‑GCTTGCAAGTCTTCCTGCACCACCTCGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTTAGCGTTGCTAAACAATCGCGCGTtat                                                                                                                                                                                                                              >  1:25020/1‑250 (MQ=255)                                gCAGAAGGAATCTGTGTGGTTTGCGTCGGGAACAGATGGCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCTATTGCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCGGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCggcggc                                                                                                                                                                                                                     <  2:261045/250‑1 (MQ=255)                                 cAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGccc                                                                                                                                                                                                                   >  2:50314/1‑251 (MQ=255)                                          tCTGTTTGGTTAGCGTCGGGACCCCCTGTGCCAGCGCGATATAAGTGGGGCGTGCTGCCAGAGCGACGTCGATTTCCATATCGTCATGCGTCGAAACGCACAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATAAGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGTAGGCCCacagcattt                                                                                                                                                                                                          >  2:170043/1‑243 (MQ=255)                                               ttggtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccgc                                                                                                                                                                                                     >  1:153940/1‑251 (MQ=255)                                                  gtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGGGGGGCGTGCTGCCAGCGCGACGTCGATTCCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCTTTATAGCGGCGGTCCAGCGCAATTGCCGCCAc                                                                                                                                                                                                  <  1:79224/251‑1 (MQ=255)                                                     tgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGAc                                                                                                                                                                                               >  2:129145/1‑251 (MQ=255)                                                        ggcGGGAACACGTGCCCCAGCGCGATATAAGAGGGGCGGGCTGCCAGCGCGACGTCGATTCCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCGGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATc                                                                                                                                                                                            <  1:170907/249‑1 (MQ=255)                                                         ccgggaccgctggTCCCAGCGCGATATAAGGGGGGCTTGCGCCCAGCGCGACGTCGTTTTCAATATCGCGAGGGGTCTAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGAGCGGGG‑GCTTGCAAATCTTCCTACCCCAAATGGACGC‑CATACGCCTGATGCTTGAGCGCCAGCCGCCAGTAAGCGTTGATAAACAAGCGTGCGTTATAGCGGCGGTCCAGCGCAATTGCCGCGACGACATCg                                                                                                                                                                                           <  2:306976/239‑1 (MQ=255)                                                           gggAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTGGAAACGCCCGGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCGGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGc                                                                                                                                                                                         <  1:83933/251‑1 (MQ=255)                                                           gggAACACATGACCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGACGCCACGACATCGGc                                                                                                                                                                                         >  1:71411/1‑251 (MQ=255)                                                                     gTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGACATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAACTGCCGCCACGACATCCGCTTCCACctct                                                                                                                                                                               >  2:86745/1‑251 (MQ=255)                                                                           gcgcGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTCAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCCTcg                                                                                                                                                                         >  1:162501/1‑251 (MQ=255)                                                                              cgATATAAGAGGGGCTTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACAACGGCTTCCCCCTCTTCATCgcgc                                                                                                                                                                      >  2:11382/1‑251 (MQ=255)                                                                                  ataAGAGGGGCGTGCTGCCAGCGCGACGTCGATTCCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGa                                                                                                                                                                   <  1:210284/250‑1 (MQ=255)                                                                                            cGTGCTGCCAGCGCGACGTCGCTTTCAATATCGTCATGTGTCGAACCGCACAGCCGCAGGACTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGAGGCCCAGCGCAATTGCCGCCCCGCAATCGACTTCCACCTCTTCCTCGCGCCGATCTTTGATGCg                                                                                                                                                        >  2:61114/1‑251 (MQ=255)                                                                                                  gcCATCGCGGCGGCGATTTCCATATCGCAATGTGGCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCt                                                                                                                                                   <  1:225695/250‑1 (MQ=255)                                                                                                   ccAGCGCGACGTCGATTTCCATATCGTAATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTgg                                                                                                                                                 <  1:164874/251‑1 (MQ=255)                                                                                                       cgcgACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGgag                                                                                                                                               >  1:283521/1‑249 (MQ=255)                                                                                                       cgcgACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGGTATAGCGGCGGCCCAGCGCAATTGCCGCCACGGCATCGGCTTCCACCCCTTCATCGCGCCGCTCTTTGAATCGTAGCTGgaggg                                                                                                                                             >  2:186388/1‑249 (MQ=255)                                                                                                           gcGTCGATTTCCATATCGTCATGTGTCGAAACTCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACCC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAgtacg                                                                                                                                         <  1:57453/250‑1 (MQ=255)                                                                                                               caaTTTCCATATCGCCATGTGTCGAAACGCCCAGCCGCAGGCTGGCCCCGCGGATGGCATGGAGATGGGTG‑GCTTTAAAATCTTCCGGCCCCAAATGGACGC‑CATTCGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCTCCACGACATCGGCTGCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAgtacgtac                                                                                                                                      <  1:307297/248‑1 (MQ=255)                                                                                                                    tCCATATCGTCATGTGTCGAACCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGACTg                                                                                                                                 <  1:95994/250‑1 (MQ=255)                                                                                                                       acttcgacagttgtagaaccGCCCCGCCGCAGGCCCGCCGCGCGGATGGCTTTGAGATCG‑TGCGCTTGTAAGTCTTCATGCCCCAAATGGACGC‑CATACGTCTGATGTTTGACCCCGAGCCGCCAGGAAGCGTTGAGAACCAATCGCGCGGTATAACGGCGGCCCAGCGCAATTGCCGCCACGACATCTGCTTCCGCCTCTTCATCGCTCCGATCGTTGATGCGTAGCTGGAGATTACGTACGCCTGCATc                                                                                                                             <  2:205673/232‑1 (MQ=255)                                                                                                                           cGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCCGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGGACGCCTGCCTCCAAc                                                                                                                         >  1:187065/1‑251 (MQ=255)                                                                                                                                 gtgtCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGAGTTGATAAACAATCGCGCGTTATCCCGGCGGCCCATCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATATTTTATGCGGAGCTGGGTAGGACGGACGCCTgactcacacagaagg                                                                                                                   >  2:61766/1‑237 (MQ=255)                                                                                                                                       aaaCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATc                                                                                                             >  1:268753/1‑251 (MQ=255)                                                                                                                                           gCCCAGCCGCAGGCCTGCCGCGCGGATGGCCTTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACCTCGGCTTCCACCTCTTCCTCGCGCCACTCTTTTATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACCGACGATCGATCCACt                                                                                                         >  2:132430/1‑251 (MQ=255)                                                                                                                                               ccccGCAGGCCGGCCGCGCGGTGGGCATTGAGATCGGTG‑CCTTGCAAATCTTCCTGCCCCAAATGGTCGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTAc                                                                                                     <  2:57767/249‑1 (MQ=255)                                                                                                                                                  cGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTAAGCCTGCATCCAACAGACGTTCGATCCACTGTACGCt                                                                                                  >  2:158458/1‑251 (MQ=255)                                                                                                                                                        ccggCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTCCCTGCCCCAAATGGACGC‑CATCCGCCTGAGGCTGGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCTGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGACCCACTGTACGCTGTccac                                                                                            <  1:42407/248‑1 (MQ=255)                                                                                                                                                        ccTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGGGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTccac                                                                                            >  1:147204/1‑251 (MQ=255)                                                                                                                                                              gcgcgGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCATGCCCAACATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGAAATTGCCGCAACGCAATCGGCTTCCTCCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGGGGGGACGTACCACTGCATCCAACAGACGTTCGATCCACTGGCCGCTGTccacgcccgg                                                                                      >  2:8823/1‑245 (MQ=255)                                                                                                                                                                  ggATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTAc                                                                                  >  2:108874/1‑251 (MQ=255)                                                                                                                                                                  cggtgtgctttgagagggggggtatGCAAATCTTCCTGCCCCAAAGGGACGCACA‑CCGCCGGTGGCTTGACCGCCGGCCGCCAGTAAGCGTGGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCAGCTCTTCACCGCGCCGATCTTTAATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTACGATCCGCTTGACGCTGTCCAGCACCGGGTa                                                                                   <  2:68761/227‑1 (MQ=255)                                                                                                                                                                    aTGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAg                                                                                <  1:129145/251‑1 (MQ=255)                                                                                                                                                                                         ttGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTAc                                                              >  1:208324/1‑249 (MQ=255)                                                                                                                                                                                         ttGCAAATCTCCCTGCCCAAAATGGACCC‑TATCCCCCTGATGCTTGTACGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGTGTACAGTCCTGAACGAAAAGGTACa                                                             <  2:231870/250‑1 (MQ=255)                                                                                                                                                                                          tGCAAATCTTCCTGCCCCAAAGGGCCGC‑AATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACag                                                            <  1:253148/250‑1 (MQ=255)                                                                                                                                                                                            cAAATTTTCCCGCCCCAAATGGACCC‑CATACGCCTGATGCTTTACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCTTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACATTCCTGAACGAAAAGGTACaggag                                                         <  2:278282/251‑1 (MQ=255)                                                                                                                                                                                                   tcctccCCCAAACGGCCGC‑CATGCTCCTGAGGCTTGATCGCCAGCCGACAGTAAGCGTTGATAAACAATCGGGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCTCCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGGGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGaa                                                      <  1:246742/247‑1 (MQ=255)                                                                                                                                                                                                   tCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATc                                                  >  1:140946/1‑251 (MQ=255)                                                                                                                                                                                                                ggACGC‑CACACGACTGATGCTTGATCCCCAGCCGCCCCTAAGCGCTGATAAACAATCGCGCGTTATCGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGGTTCCACCTCTTCAGCGCGCCGATCTTTGGTTCGCAGCTGGAGAGTACGTAAGCCTGCATCCAACAGACGTTCGATCCCCTGTACGCTGTCCACCACCGGGTACAGTACTGGACGCACAGGTACAGGAGGAAAATCAGGCTGATACATc                                     >  2:62175/1‑251 (MQ=255)                                                                                                                                                                                                                        aaccccctgatttgtgctcccccgccgccAGTAAGCGTTGATAAACAATCCCGCGTTATAGCGGCGGCCCAGCGTAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCTCTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtc                              <  2:229740/228‑1 (MQ=255)                                                                                                                                                                                                                          aCGCCTGATGCTTGATCGCCAGCCGCCAGGAAGCGTTGATAAACAATCGCGCGTTATAGCGACGGCCCAGCGCAATTACCGCCACGACATCGGCTTCCAACTCTTCATCGCGCCGCTCTTTGATGCGTAGCTGGAGAGTACGTACGCCCGCATCCAACAGACGTTCGATACACTGTACGCTGTCCACCACCGGGTACAGCCCTGAACGAAAAGGTAAAGGAGGAAAATCAGACTGATACATCCCGCAtcct                            >  2:191257/1‑251 (MQ=255)                                                                                                                                                                                                                                gATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCtt                        <  1:270409/249‑1 (MQ=255)                                                                                                                                                                                                                                gATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCtt                        <  1:306066/249‑1 (MQ=255)                                                                                                                                                                                                                                         ccGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCGCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGa             <  2:49681/250‑1 (MQ=255)                                                                                                                                                                                                                                          cGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGAt            <  1:215252/251‑1 (MQ=255)                                                                                                                                                                                                                                                    cAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTAGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTACTTACGCAGGGAGATTTCTCCGccc  >  2:272549/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     ATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTG‑GCTTGCAAATCTTCCTGCCCCAAATGGACGC‑CATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCT  >  NC_000913/4193749‑4194240

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