Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,912,074	G→A	S532F (TCC→TTC)	glmS ←	L‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,912,074	0	G	A	100.0%	152.2 / NA	43	S532F (TCC→TTC)	glmS	L‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/17);  total (26/17)

ATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCCTGA  >  NC_000913/3911834‑3912319                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      aTTTATTAATCAACCGCAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGttt                                                                                                                                                                                                                                             >  1:83321/1‑251 (MQ=255)               cGTAACCGCTTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCgt                                                                                                                                                                                                                               >  1:138176/1‑251 (MQ=255)                gTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCgtt                                                                                                                                                                                                                              <  1:46271/251‑1 (MQ=255)                                 gTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGcaacaacaa                                                                                                                                                                                                             <  1:15973/251‑1 (MQ=255)                                      gcgGCTGGTCAAAGTCGGTGCCTTTGATCAGCGGGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGGTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAAcc                                                                                                                                                                                                        <  2:245782/251‑1 (MQ=255)                                                      ggTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCa                                                                                                                                                                                         >  1:215553/1‑250 (MQ=255)                                                                  cAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCg                                                                                                                                                                            >  2:160167/1‑251 (MQ=255)                                                                   agagcgACAGGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTAACGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCgg                                                                                                                                                                           <  2:299909/248‑1 (MQ=255)                                                                                                 ggAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGc                                                                                                                                             >  2:295670/1‑251 (MQ=255)                                                                                                     cGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTc                                                                                                                                          >  1:331862/1‑250 (MQ=255)                                                                                                         gtAGAAGATCGGTGCAAGCACCTCTGCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGt                                                                                                                                     <  2:27797/251‑1 (MQ=255)                                                                                                           agaagaTCGGTGCAATCACCTCTTCCACATGCGGTATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTg                                                                                                                                    <  1:77758/250‑1 (MQ=255)                                                                                                           agaagaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTg                                                                                                                                    <  1:141184/250‑1 (MQ=255)                                                                                                           agaagaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTg                                                                                                                                    >  2:8392/1‑250 (MQ=255)                                                                                                            gaagaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGa                                                                                                                                   >  2:140269/1‑250 (MQ=255)                                                                                                                   ggTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCata                                                                                                                                                                                                  >  1:289580/1‑180 (MQ=255)                                                                                                                   ggTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACATACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCata                                                                                                                                                                                                  <  2:289580/180‑1 (MQ=255)                                                                                                                     tGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAag                                                                                                                            >  2:71222/1‑248 (MQ=255)                                                                                                                                       aTGCGGCATCTCGATGATGCGCATGTTATCGCTACTTACAAAACCCGCCTCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCCATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATg                                                                                                       >  2:40413/1‑251 (MQ=255)                                                                                                                                          cGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATgcgc                                                                                                    >  2:127324/1‑251 (MQ=255)                                                                                                                                                ctcGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTcc                                                                                               >  1:148362/1‑250 (MQ=255)                                                                                                                                                 tcGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCa                                                                                                          >  1:326069/1‑238 (MQ=255)                                                                                                                                                 tcGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCa                                                                                                          <  2:326069/238‑1 (MQ=255)                                                                                                                                                       atgggCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGAt                                                                                       <  1:199608/247‑1 (MQ=255)                                                                                                                                                          tgCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTgg                                                                                    >  1:287260/1‑251 (MQ=255)                                                                                                                                                                  aTCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTACGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGAt                                                                            >  2:340657/1‑251 (MQ=255)                                                                                                                                                                  aTCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGa                                                                             <  2:364207/250‑1 (MQ=255)                                                                                                                                                                           tACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACgg                                                                    <  2:48272/250‑1 (MQ=255)                                                                                                                                                                               aaaCCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGTATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGccc                                                                 >  2:291258/1‑249 (MQ=255)                                                                                                                                                                                   ccGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAAc                                                           <  1:390750/251‑1 (MQ=255)                                                                                                                                                                                   ccGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAAc                                                           >  2:209064/1‑251 (MQ=255)                                                                                                                                                                                    cGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAAc                                                           >  2:132348/1‑250 (MQ=255)                                                                                                                                                                                      cATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAgc                                                        >  2:397005/1‑251 (MQ=255)                                                                                                                                                                                       aTCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAg                                                                                                                                                                                    >  1:8675/1‑126 (MQ=255)                                                                                                                                                                                       aTCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAg                                                                                                                                                                                    <  2:8675/126‑1 (MQ=255)                                                                                                                                                                                                                tGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAAtcttc                              >  1:106481/1‑251 (MQ=255)                                                                                                                                                                                                                      ccACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAg                        <  2:350026/251‑1 (MQ=255)                                                                                                                                                                                                                        aCGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAgcg                      >  2:186135/1‑251 (MQ=255)                                                                                                                                                                                                                                 aaCTTCTTCAGTGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGc             >  1:170006/1‑251 (MQ=255)                                                                                                                                                                                                                                 aaCTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGaaa                                   <  2:28738/229‑1 (MQ=255)                                                                                                                                                                                                                                 aaCTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGaaa                                   >  1:28738/1‑229 (MQ=255)                                                                                                                                                                                                                                      cttcAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTg        >  2:344897/1‑251 (MQ=255)                                                                                                                                                                                                                                             gTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCCTga  <  1:186135/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      ATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCCTGA  >  NC_000913/3911834‑3912319

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