Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,322,325	G→A	L208F (CTC→TTC)	trpE ←	anthranilate synthase subunit TrpE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,322,325	0	G	A	100.0%	149.6 / NA	47	L208F (CTC→TTC)	trpE	anthranilate synthase subunit TrpE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/27);  total (20/27)

GATTACTCTTTTTCAGCACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTG‑CGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAA  >  NC_000913/1322082‑1322562                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 gATTACTCTTTTTCAGCACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACgttgt                                                                                                                                                                                                                                           >  1:200631/1‑251 (MQ=255)        ctTTTTCAGCACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTttctt                                                                                                                                                                                                                                    >  2:458135/1‑251 (MQ=255)           tttCAGCACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTttcttct                                                                                                                                                                                                                                  <  2:200631/250‑1 (MQ=255)            ttCAGCACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTttcttctt                                                                                                                                                                                                                                 <  1:30551/250‑1 (MQ=255)             tCAGCACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAAGCGAGCAGTGAAACGTTGTTTttcttcttc                                                                                                                                                                                                                                <  1:423153/250‑1 (MQ=255)                      aaTAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACAGTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGTCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGGTCAGGCGAGCAGTGAAACGTTGTTTTTTTTCTTCATT‑CGGAGc                                                                                                                                                                                                                      <  2:542377/250‑1 (MQ=255)                        tAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACAACTGGCAGCGGCGGCGCGGCTTCGGACAGTTGCTGACGTAGTTCGTTCAGGCGAGCCGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCac                                                                                                                                                                                                                    >  2:366709/1‑249 (MQ=255)                                 aGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGc                                                                                                                                                                                                          >  1:295160/1‑251 (MQ=255)                                    ggTGACGGGCAGTGCATAGGTAGAAGGCGAGATGGCCCCACCTGGAAAATTACTCCAGCGCGAATCGCTTTTTGCATCAAACGCTCTACGCCACCGAACTCATCATCGCTCTGAGTACATTCACAGCGCATATGCCGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCGGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCGTCATT‑CGGAGCAAACAGGCTGGCCTg                                                                                                                                                                                                       <  2:478789/251‑1 (MQ=255)                                       gACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGaa                                                                                                                                                                                                     <  2:520217/250‑1 (MQ=255)                                              agggcagAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTg                                                                                                                                                                                             >  1:37902/1‑251 (MQ=255)                                                  cagAGAGAAACGGCGAGATGGCACTACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCGCTACGCCACCGAACTCTTCATCGCTCTGATTACAGTCACAACGCATATGCGGCACTGAAACCACTGGCATCGGCGGCGCGGCTTCGGTCAGTTGCTGTCGTAGTTCGTTCAGGCGAGCAGTGAAACTTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCttt                                                                                                                                                                                         <  1:431471/251‑1 (MQ=255)                                                    gagagaAACGGCGAGATGGCCCCACCTAGAAAATTTCTCCAGCGCGAATCGCTTGTTGCAACAAACGCACTACGCCACCGAACTCTGCAACGCTCTGATTACATTCACAACGCATATGCTGCACGGAAACCACTGGCAGCGGCGGAGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGATGGCCTTAATACGGGTGCttttt                                                                                                                                                                                       <  2:103066/251‑1 (MQ=255)                                                     agagaAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCtttttt                                                                                                                                                                                      <  1:220323/251‑1 (MQ=255)                                                     agagaAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCtttttt                                                                                                                                                                                      <  1:175236/251‑1 (MQ=255)                                                      gagaAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCttttttt                                                                                                                                                                                     <  2:358724/251‑1 (MQ=255)                                                       agaAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTc                                                                                                                                                                                    <  1:458135/251‑1 (MQ=255)                                                               cGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTc                                                                                                                                                                           >  1:447871/1‑251 (MQ=255)                                                                gagaTGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCa                                                                                                                                                                          <  1:105955/251‑1 (MQ=255)                                                                 agaTGGTACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCa                                                                                                                                                                          >  1:144962/1‑250 (MQ=255)                                                                                gAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAg                                                                                                                                                                                                                       <  2:63152/191‑1 (MQ=255)                                                                                gAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAg                                                                                                                                                                                                                       >  1:63152/1‑191 (MQ=255)                                                                                 aaaaTTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTc                                                                                                                                                         >  2:472223/1‑251 (MQ=255)                                                                                  aaaTTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCAACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGGAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCt                                                                                                                                                                                           >  1:336714/1‑217 (MQ=255)                                                                                  aaaTTTCTCCAGCGCGAATCGCTTTTGGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCGCTGGCAGCGGCGGCGCGGCGTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTGCTTCGTCATT‑CTGAGCAAACAGGCTTGCCTGAATACGGGTGCt                                                                                                                                                                                           <  2:336714/217‑1 (MQ=255)                                                                                    aTTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGc                                                                                                                                                      >  2:197050/1‑251 (MQ=255)                                                                                          ccAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCgag                                                                                                                                                   <  2:406095/248‑1 (MQ=255)                                                                                                cgAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCCTCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAAc                                                                                                                                           >  1:458445/1‑250 (MQ=255)                                                                                                            tGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGTAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGGTCGTTCAGGCGAGCAGTGAAACGGTGTTTTGCTTCTTCATT‑CGTAGCAAACAGGCTGGCCTGAATTCTGGTGCTTTATTTCTG‑ATTGTCAATCACCAGCAGCGTTTCAGCGATATAAAAACAGAAATCAGGGc                                                                                                                               <  1:487730/250‑1 (MQ=255)                                                                                                             gCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTGTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGGGCGTGTTTTCTG‑ATGGGCAATCACCAGCAGCGGTTCAGCTAGATAATAACAGAAATTAGGGCa                                                                                                                              <  2:205411/250‑1 (MQ=255)                                                                                                                 caaACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAt                                                                                                                         <  1:84673/251‑1 (MQ=255)                                                                                                                      gcactaggcccgcGAACGCTTCAGCGCTCTGATTACATTCACAATGCCTATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGACGTCGGTCAGTTGCTGACGTAGTTCGTTCGGGCGGGCAGTGAAACGTTGGTTTTCTTCTTCATG‑CTGGGCAAACAGGCTGGCCTGAAGACGGGTGATTTGTGTCTG‑ATGGTGAGTCACGATCGGCGTGTCAGCGAGATATAAACAGATATCAGGGCCTTTATTTTcc                                                                                                                    <  2:540943/239‑1 (MQ=255)                                                                                                                      ccACTACGCCACCGAACTCTTCATCGCTCTGATTACAGTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGTCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCATGGCAGTTATTTTcc                                                                                                                    <  1:433439/250‑1 (MQ=255)                                                                                                                       cACTACGCCGCCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCg                                                                                                                   >  2:401538/1‑251 (MQ=255)                                                                                                                         cTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGCTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCt                                                                                                                 >  1:486858/1‑251 (MQ=255)                                                                                                                                         ttCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTG‑CGGTaaa                                                                                                 >  2:426450/1‑250 (MQ=255)                                                                                                                                                      ttACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTG‑CGGTAAATCTTCAAATCCCGc                                                                                   >  2:248815/1‑251 (MQ=255)                                                                                                                                                                             cACGGAAACCACTGGCGGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTGCAGGCGAGCAGTAAACCGTTGTTTTTCTTCTTCATG‑CGGAGCAAACAGGCTGGCCTGAATGCGGGTGTTTTTTTTCTG‑ATGGCAAATCACCATCAGCGTTTCAGCGAGAAAAAAACAGAAACCAGGGCAGTTATTTCCCGCTGACAGTGG‑CGGTAAACCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGgc                                                            <  1:472223/251‑1 (MQ=255)                                                                                                                                                                              aCGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTG‑CGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGgcc                                                           >  1:242662/1‑251 (MQ=255)                                                                                                                                                                                          ggCAGCGGCGGCGCGGCTTCGGGCCGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGGTG‑CGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAAcc                                                >  2:374873/1‑249 (MQ=255)                                                                                                                                                                                                  gggcgcggcTTCGGTCATTTGGTGGTGTAGTTCTTGCGGGCGAGCAGTGAAACTTTGTTTTTTTTTTTCATT‑GGGAGCAAACAGTCTGCCCTAAATACGGTGGTTTTTTTTCTG‑AGGGTCAATCACCATCAGCTTTTCAGCGAGAAAAAAACAAAAATCAGGGCAGTTATTTGCCGCTGACAGTTG‑CGGTAAATCTTCAAACCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTtctc                                       <  1:426450/250‑1 (MQ=255)                                                                                                                                                                                                  cggcgcggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGCGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTG‑CGGTAAATCTTCAAATCCCGCCACAATGTCATAAGAGAACAGTCCGCCGAAGAACATGGCTtctc                                       >  2:447119/1‑251 (MQ=255)                                                                                                                                                                                                   ggcgcggcTTCGGTCAGTGGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACTTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGAAAAAAACAGAAACCAGGGCAGTTATTTCCCGCTGACAGTGG‑CGGTAATTCTTCAAACCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCg                                      <  1:197050/251‑1 (MQ=255)                                                                                                                                                                                                                  aGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTG‑CGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTAc                       >  1:86874/1‑251 (MQ=255)                                                                                                                                                                                                                              aGTTCGTTCAGGCGAGCAGGGAAGCGTTTTTTTTCTTCTTCATT‑CGGAGCAAACGGGCTGGCCTGAATACGGGTGTTTTTTTTCTG‑AGGGTCAATCCCCATCAGCTTTTCAGCGAGAAAAAAACAGAAATCAGGGCAGTTTTTTCCCGCGGACATTTG‑CGGAAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGAt            <  1:248815/250‑1 (MQ=255)                                                                                                                                                                                                                                  tgttcggtagatcAGTGAAATGTTTTTTTTCTTCTT‑ATTGCGGAGCAAACAGGTTGCCCTGAATACGGTGGCTTTTTTT‑TGTAGGGTAAACCACCATCAGCGTTTAAGGAAGATAAAAACAGAAATCAGGGCTGTTTTTTCCCGCTGTC‑GTTGGCGGTAAATTTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTg        <  1:248000/238‑1 (MQ=255)                                                                                                                                                                                                                                       gggCGAGCAGTGAAACTTTTTTTTTCTTCTTCATT‑GGGAGAAAACAGGCTGGCCTGAATCCGGGTTTTTTTTTTCTG‑TTGTTCAATCCCCATCAGCGTTCCAGCGAGATAAAAACAAAAACCAGGGCAGTTATTTTCCGCTGACAGTTG‑CGTAAAATCTTAAAACCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATaa  <  1:401538/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 GATTACTCTTTTTCAGCACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATT‑CGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTG‑ATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTG‑CGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAA  >  NC_000913/1322082‑1322562

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