Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,592,188	C→T	H178Y (CAT→TAT)	tsr →	methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,592,188	0	C	T	100.0%	141.7 / NA	44	H178Y (CAT→TAT)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/23);  total (21/23)

TTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCG‑ATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCAT  >  NC_000913/4591944‑4592430                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      ttCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGat                                                                                                                                                                                                                                                  >  1:176451/1‑250 (MQ=255)    gCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCg                                                                                                                                                                                                                                              <  1:134409/251‑1 (MQ=255)                ggAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGAt                                                                                                                                                                                                                                   >  1:450426/1‑250 (MQ=255)                     cgtgcggggaagATTACGAAGCGTGGCCGCTTGACCCGCGCCAGATCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAAt                                                                                                                                                                                                                             <  2:400388/240‑1 (MQ=255)                           gggcagtttACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGcctc                                                                                                                                                                                                                        <  1:449140/243‑1 (MQ=255)                           gggCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAAt                                                                                                                                                                                                                             >  1:400388/1‑245 (MQ=255)                                               ccGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCCACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGCAATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGAt                                                                                                                                                                                                   >  1:34551/1‑251 (MQ=255)                                                                 agCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTgatgat                                                                                                                                                                                 <  1:142926/251‑1 (MQ=255)                                                                  cccccgccgCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtc                                                                                                                                                                                <  2:450426/248‑1 (MQ=255)                                                                     ccgccgCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGACTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCTTGATGAtcgtc                                                                                                                                                                             >  1:342616/1‑251 (MQ=255)                                                                         cgcAGCGGCAGAGAACAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTAGGATATCGCCGTAAGCGATAACAATGCCTCCTACAGCCAGGGGGTGGGGGTTCTGGTGGGCGGGCTGATGGTCGTAc                                                                                                                                                                         >  2:10024/1‑251 (MQ=255)                                                                                  agagATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTAAAGCCAGGCGATTGGGATTCTGGGGGGCGAGATGATCGTCGTACTGGCGGTc                                                                                                                                                                 >  2:495881/1‑250 (MQ=255)                                                                                     gaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTAc                                                                                                                                                                                                                    >  1:177647/1‑196 (MQ=255)                                                                                     gaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTAc                                                                                                                                                                                                                    <  2:177647/196‑1 (MQ=255)                                                                                                  tACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCAGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGAGGATTCTGGTGGGCGAGGTGATCGTCGTACTGGCGGGAATCTTTGCCG‑TCTgggt                                                                                                                                               >  1:241026/1‑249 (MQ=255)                                                                                                       tatTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCCTCTTCCCCG‑GCTGGGTCgggt                                                                                                                                          >  1:127996/1‑249 (MQ=255)                                                                                                        atTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGCTAACAATGCCTCCTACAGCCAGGCGATGTGAATTCTGGTGGGCGTGATGCTCGTCGTACTGGCgggacttttctccggcgcgtgggtggg                                                                                                                                          >  2:188815/1‑226 (MQ=255)                                                                                                           tATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATATTCTCCG‑GCTGGGTCGGTAttt                                                                                                                                       >  1:415850/1‑249 (MQ=255)                                                                                                              cacaATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCGGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGCTTCTGGTGGGCGTGATGATCGTCGTACTTGCGGTCATCTTCGCCGCTCGGGTTGGGT‑TTAAAGc                                                                                                                                   >  1:40546/1‑251 (MQ=255)                                                                                                                    gcgcTGGCGGAGCTGATCCAACTAATAGGTGCAGGCAAAATCAACGAGGTCTTTGATCAGCCGGCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCGTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCt                                                                                                                             <  2:75260/251‑1 (MQ=255)                                                                                                                                    tCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGCATCTTCGCCG‑CCTGGGTTGGTATTAAAGCCTCTCTGGTAGCG‑cccataagt                                                                                                            >  1:439129/1‑244 (MQ=255)                                                                                                                                     aacacttttAGGTGCAGGCAAAAGCAACGAGTTCTGTTATCAGCAGCCCCAGGGATATCCGGACGGTTTCGATAATCAGTATGTGGCTTACATGGATCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCCATGTTGAGTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCTG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGTATc                                                                                                           <  2:395403/244‑1 (MQ=255)                                                                                                                                           attAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCt                                                                                                       <  2:11978/248‑1 (MQ=255)                                                                                                                                              aGGTGCAGGCAAAAACAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGAtt                                                                                                   <  2:127996/250‑1 (MQ=255)                                                                                                                                                 ggCAGGAAAAATCAACTAGTTTGTTGATCAGCCGACCCAGGGATATCAGGACGGTTGCGAGAAGCATTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGTATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTTGTAGCG‑CCAATGAATCGCCTGATTGACa                                                                                               <  2:370718/250‑1 (MQ=255)                                                                                                                                                   cAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTCCCG‑TCTTGTTCGGTATTAAAGCCTCGCTgggtgcggcactgaatagcccgttt                                                                                                    >  2:405680/1‑222 (MQ=255)                                                                                                                                                   cAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGAtt                                                                                                   <  1:405680/245‑1 (MQ=255)                                                                                                                                                      gCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGAAGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATTTTCCCCG‑GCTGGTTCGGGATTAAAACCTCACTGGTAGCGACCAA‑GAAACGCCCGGTTGACCGCAtt                                                                                          >  1:146031/1‑251 (MQ=255)                                                                                                                                                                  aggtCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACTATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCa                                                                              <  2:415850/248‑1 (MQ=255)                                                                                                                                                                  aGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCa                                                                              <  2:188371/251‑1 (MQ=255)                                                                                                                                                                        tggATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCGTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAggc                                                                           <  2:146031/246‑1 (MQ=255)                                                                                                                                                                          gATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATc                                                                                                                                                              >  1:54976/1‑165 (MQ=255)                                                                                                                                                                          gATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATc                                                                                                                                                              <  2:54976/165‑1 (MQ=255)                                                                                                                                                                                         ggATATCATGACGGTTGCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCTGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCG‑At                                                      <  2:485162/251‑1 (MQ=255)                                                                                                                                                                                         ggATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCATTCGTCat                                                                                    <  1:445818/222‑1 (MQ=255)                                                                                                                                                                                                   gcGGTTTCGAGAAGCAGTATGTGGCCTACATGGAGCAAAACGATTGGCTCTATGATTTCGCCGTCAGCGATAACAATGCCTCCTCCAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAGCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCATGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCG‑ATTGAGGTGGAt                                            <  2:263690/250‑1 (MQ=255)                                                                                                                                                                                                      gTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGGTTTTAAGCCTCGCTGGTAACG‑GCAATGAATCGCCTGATTTAAAGGATTCTTCATATTTCAGGCGGCGGTCTGGTGAAAACG‑ATTGAggggggtggg                                         >  1:391092/1‑243 (MQ=255)                                                                                                                                                                                                         tCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGACGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCAGTCGGCATATTGCAGGCGGCGATCTGGTGAAACCG‑ATTGAGGTGGATGGctc                                       <  1:24350/250‑1 (MQ=255)                                                                                                                                                                                                             gaAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑GCTGGTTTGGTATTTAAACCTCCGTGGTAGCG‑CCCATGAATCGCCTGATTGAAAGCCTTCGTCCTATTTCAGGCGGCGATCTGGTGGAACCGCATTG‑GGGGGGTGGCTCTAATg                                  >  1:169961/1‑251 (MQ=255)                                                                                                                                                                                                                    tATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATTTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATAGCCTGCTTGGCAGCATTCGTCATATTTCAAGCGGCCATCTGGGGTAAACG‑ATTTGGGTGGGTGGGTCTAATTGGATggg                           >  1:454793/1‑251 (MQ=255)                                                                                                                                                                                                                    aatttgGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACCATGCCGCCTACAGCCAGGCGAGGTGGATTCTGGTTGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCATGATTTACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCG‑ATTGAGGTGGATGGCTCTAATGAGATggg                           <  2:455289/247‑1 (MQ=255)                                                                                                                                                                                                                              aCATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCG‑ATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGc                  <  1:405989/250‑1 (MQ=255)                                                                                                                                                                                                                                      cAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGGTTGGGATTAAAAACTCGCTTGTAGCG‑CCAATGAATCGCCCGTTTGGCAGCATTCGTCCTATTGCAGGCGGCGGTCTGGTGAAACCG‑GTTGAGGGGGGAGGCTCTAATGAGATgggggcacaggcagggaggtt         >  1:357649/1‑234 (MQ=255)                                                                                                                                                                                                                                             aTCGGCGCTATGATATCGCCGGCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCG‑ATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCat  <  2:36683/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      TTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCG‑CCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCG‑ATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCAT  >  NC_000913/4591944‑4592430

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