Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	557,635	A→T	L36Q (CTG→CAG)	folD ←	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	557,635	0	A	T	100.0%	114.2 / NA	39	L36Q (CTG→CAG)	folD	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (20/19);  total (20/19)

GATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACT  >  NC_000913/557396‑557882                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        gATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATTAATTTGCGATGCAGGGTTACTACCCACCCGCACAACGGCCTGTCCTGGTGcc                                                                                                                                                                                                                                              >  2:156861/1‑251 (MQ=255)          gTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCg                                                                                                                                                                                                                                     >  1:185452/1‑251 (MQ=255)                     gtgaccattaccatccccGCCGTTAACGGCAGTTGAACCAGAATGCCATCGTTTGGGTTGGCGGCATGCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTGCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCATGGTTGCTGCCCACCATCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATAc                                                                                                                                                                                                                           <  2:205991/235‑1 (MQ=255)                         atgttctatcaaccgccgGTAACGGCAGTTGACACAGAATGCCTTCTGTGGTGTGGTCGGCATTCAGCGTACCGATAAGCCCCAGCATCGCTGCTTCGCTGGGGGTTGCCGGTAGGTCAGAATAGCGGGCGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGTCATAAATTTGCGATGCAGGGTTACTACCCACCAGCCCAACGGCCTGCCCTGGTGTCCGCCGTCCGGcataaataagcgcc                                                                                                                                                                                                                      <  2:406671/239‑14 (MQ=255)                            tATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTg                                                                                                                                                                                                                    >  1:384422/1‑250 (MQ=255)                                    cccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTg                                                                                                                                                                                                           <  1:239653/251‑1 (MQ=255)                                          ggTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCCCCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCACTACGCGCCTGAACTTTTTGAGCAAc                                                                                                                                                                                                     >  1:322310/1‑251 (MQ=255)                                                  cAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAgc                                                                                                                                                                                             <  1:312353/251‑1 (MQ=255)                                                    gTTGAACCAGAATGCCATCGATGGAGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGACATAAGAGCGGGAGGCGAACCCGACTTCT‑TCACAAGCCTTGAGTTTGCTTGCGACATAAGTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCATGTCCTGGGGCCCGCGGTCCGGCTGCAATAAGCGCCTGAACTTTTTGAGCAACTTCAGAgcgc                                                                                                                                                                                           >  1:100457/1‑251 (MQ=255)                                                        aaCCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCagaacgaccct                                                                                                                                                                                       >  1:453108/1‑243 (MQ=255)                                                            tgattGCCATCGACGCGGTTGTCGGCATTCAGCGACTCGCTAAGCTCCAACAGCCACGCTTCGCTGGAGGATTACGGGAGGTCTTAAGAGCGGGTGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTAGCTTGCGCAATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACAGCCTGTCCTGGTGACCGCTGTCCGGCTGCAATACGAGCCTGAAATTTTTGACCTAATTCAGACCGCAcacac                                                                                                                                                                                     >  2:52945/5‑245 (MQ=255)                                                             gAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCtgctgc                                                                                                                                                                                  <  1:218889/251‑1 (MQ=255)                                                             gAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCtgctgc                                                                                                                                                                                  >  1:375842/1‑251 (MQ=255)                                                                  ccATCGATGGTGATGTCGGCATTCAGCGAATCGTTATGATCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGAGAGGTCATCAGAGCGGGAGACGCACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATACCGGGTTACTACTCACCAGCACAACGCCCTGTCCTGGAGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTAGGCAACTTCACAGCGCACctactacgaaat                                                                                                                                                                             >  2:258298/1‑241 (MQ=255)                                                                        atgtggttgtCGGCATGCAGCGTATCGATACGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGtttt                                                                                                                                                                       <  2:124130/246‑1 (MQ=255)                                                                                 tCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCaat                                                                                                                                                              >  1:352945/1‑251 (MQ=255)                                                                                  cggcggtctgtgtacagatcagtcccAGCAGCTCGGCTTTGCTGGTGGGTTCCCGGAGGTAATCAGAGCTGGTGACGAATCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGTGACATAAATTTGCTATGCATGGTTCCTTCCCACCTGCACAACGGCCTGTTCTGGGGCTCGCAGTCCGGCTGCTATACGCGCTTGATCTTTTTTAGCAACTTCAGAGCGCACCTGCTGCTCCATCGTTTTACCGTCaata                                                                                                                                                             <  2:280632/227‑1 (MQ=255)                                                                                   ggCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGCCTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCaata                                                                                                                                                             >  1:160900/1‑250 (MQ=255)                                                                                       ttCAGCGGATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCt                                                                                                                                                         <  2:167826/250‑1 (MQ=255)                                                                                                              cagcTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTc                                                                                                                                 >  1:9984/1‑251 (MQ=255)                                                                                                                   ccGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATc                                                                                                                             <  1:81894/250‑1 (MQ=255)                                                                                                                   ccGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCt                                                                                                                            <  1:136904/251‑1 (MQ=255)                                                                                                                                   ttCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCACTACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGCGAGAGGATTCCATCTGTTACGTAGATCGacg                                                                                                            >  1:311254/1‑249 (MQ=255)                                                                                                                                     ccGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAggg                                                                                                          <  1:22101/251‑1 (MQ=255)                                                                                                                                      cGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAgggg                                                                                                         >  1:404503/1‑251 (MQ=255)                                                                                                                                           ggTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGAGACATAAATTTGCGATGAAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGCGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGCGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATgcg                                                                                                    >  1:200034/1‑251 (MQ=255)                                                                                                                                                         gagaCGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGa                                                                                      <  1:100781/251‑1 (MQ=255)                                                                                                                                                          agaCGAACCCGACGTCTGTC‑CAATCCTTGCGGTTGGTTGCGACACAAGTTTGCGATGGTGGGGTTCTCGCCGCCAGCACAACGGGCTGTCCGGGTGCCCTCAGTCCGGCTGCAATACGCAGCTGATCTTTTTGAGGAACTTCAGAGCGCACAGGCTGCGCAATCGTTTTACAGTGAATAATCATTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATTGAAGTGGATGCGTCTATTTTGTCAGaa                                                                                     <  2:37632/251‑1 (MQ=255)                                                                                                                                                               aaCCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCggg                                                                                >  1:39894/1‑251 (MQ=255)                                                                                                                                                                    aaCTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGACCTGTCCTGGTGCCCGCCGTCCGGCTACAATACGCGCCTGAACTTTTTGAGACACTTAAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAAAAATCTTTCCTGCCCTCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGAGCCTGTTTTGtaggaagctggacgcg                                                                           >  2:318165/2‑236 (MQ=255)                                                                                                                                                                            ccacaatcactgtggtggcttgcGACAAAAATTTGCGATGCAGGGTTACTACCCCCCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCAGCTGCAATACGCGCCTGAAACTTTTGCGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTTCCGTCAATAATCTTTGCTTCCATCAGAGATAGGACTCCATCTGTTCCGTAGATCGAAGGGGATGCGCCTACGTTGTCAGAAGCGGGGCTCGCTGTCa                                                                     <  2:413610/233‑1 (MQ=255)                                                                                                                                                                                 aGCCTTGCGATTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCCCAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCGGCTGCGCAATCGTGTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGTGTATGCGCCTATTTTGTCAGAAGCGGGTCGCGCTGTCTGgtttc                                                               <  2:200034/251‑1 (MQ=255)                                                                                                                                                                                   ccTTGCGTGTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgt                                                             <  1:439676/251‑1 (MQ=255)                                                                                                                                                                                           ttgcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGAtt                                                     >  1:153912/1‑251 (MQ=255)                                                                                                                                                                                                   gACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAAGGCCTGTCCTGGTTACCGCAGTCCCGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGAGCCTATTTTGTCAGAAGCGGGGCGCGCTGACCGGTTTAGTTTACGATTTATCGCGt                                             >  2:433086/1‑251 (MQ=255)                                                                                                                                                                                                        aaaTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAg                                                                                                            <  1:39090/183‑1 (MQ=255)                                                                                                                                                                                                        aaaTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAg                                                                                                            >  2:39090/1‑183 (MQ=255)                                                                                                                                                                                                                 aatcAGGGTTACTACCCACCCGCCCAACTGTCTGTTCCGGTGCCCGCAGGCCGGCTGCAATACGCGCTTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTg                               <  2:352945/248‑1 (MQ=255)                                                                                                                                                                                                                    cAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGa                            <  1:286885/251‑1 (MQ=255)                                                                                                                                                                                                                                              ggtaGTCATGGAGCCCGCAGGCTGTCTGCAATTCGCTCCAGAACTTGGTGGGCAACTTCAGAGCGCACCTGCTGCTCGATAGTTTTCCCGTCAATAATCTGTGCTGCCATCAGAGAGAGGATTCCATCTGCTACGTAGATCGAAGGGGATGCGCCTATTGTGTCAGAAGCGGTACGCGCTGTCAGGTTTCGTTTCAGACTTATCGCGTGCAGCGACCCCTTTCGAAGGTGAGGCGCACCGTCGCTGAGACt  <  2:100208/248‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        GATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCT‑TCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACT  >  NC_000913/557396‑557882

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