Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,875	A→G	V513V (GTA→GTG)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,875	0	A	G	95.0%	114.8 / ‑2.6	40	V513V (GTA→GTG)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (14/24);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-01

AGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTAT  >  NC_000913/1631‑2118                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    aGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGACGAATAAACATATCGACTTACGTGTATGCGGTGTTGCCAACTCGAATGCACAGCTCACCAATGTGCATGGc                                                                                                                                                                                                                                               >  2:51367/1‑251 (MQ=255)  gggATCGTCTGAACTCTCAAGCCCGGGCGTGGTAAATAACGAGGATGCGACCCCTGGCGTGCGCGTTACGCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGAGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGGGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGTCTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTTACAAATGTGCATGGcc                                                                                                                                                                                                                                              <  2:360943/251‑1 (MQ=255)                         ctGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGga                                                                                                                                                                                                                        >  1:163939/1‑250 (MQ=255)                           gTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGgaag                                                                                                                                                                                                                      <  1:333226/250‑1 (MQ=255)                                   aaaTAACGATGATGCGACCACTGGCGTGCGCGTTAATCATCAGATGCTGTTCAATACCGCTCAGGTTATCGAAGTGATTGTGATTGGCGTCGGTGGCGTAGGCGGTGCGCTGCTGGGGCAAATGAAGCGTCAAAAAAGCTGGCGGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGACAGAGCTGCTCACCCATGTGCATGGCCTAAATCTGGAAAAATGGCAGGAACAACTGtcgc                                                                                                                                                                                                             >  2:328530/1‑251 (MQ=255)                                        aCGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGcc                                                                                                                                                                                                        <  1:108430/251‑1 (MQ=255)                                                gCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCTCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGcc                                                                                                                                                                                                <  1:276576/251‑1 (MQ=255)                                                       cTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGTAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAATGCACTGCTCACCACTGTGCATGGCCTATATCTGGCAAACTGGAAGGAAGAACTGGCGCAAGCCAACGAGCCGTTTAAt                                                                                                                                                                                         >  1:14359/1‑251 (MQ=255)                                                         ggCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAtct                                                                                                                                                                                       <  1:350159/251‑1 (MQ=255)                                                           cGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCg                                                                                                                                                                                     <  1:141237/251‑1 (MQ=255)                                                                   ttACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGAAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCtt                                                                                                                                                                              >  1:126605/1‑250 (MQ=255)                                                                   ttACTCATCAGATGCTGCTCCATACCAATCAGGTTATCGAAGTGCTTATCATTGGCGTCGGTGGCGTTGACCGTGCGCTCCTCGAGCAGCTGAACCGTCAACAAAGATGGCTGAAGAATAAACATGTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGAACTGCTCACAAACGTGCCTGGCCTAAATCTAGAAAACTGGCAGGGAAAACTGGCGCAAGCCAAGGAGACGATTAATCTCgagagcct                                                                                                                                                                              >  2:252619/1‑243 (MQ=255)                                                                        catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                                         <  1:280262/250‑1 (MQ=255)                                                                                   gtccAATCCCGATCAGGTTATCGAAGTGTTTGTGATTGGTGGCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCCGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCGCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAATAACGGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGaa                                                                                                                                                               <  1:215293/246‑1 (MQ=255)                                                                                      cAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAat                                                                                                                                                          <  1:137301/251‑1 (MQ=255)                                                                                                cAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgc                                                                                                                                                 <  2:244003/250‑1 (MQ=255)                                                                                                        cGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCg                                                                                                                                         <  1:308242/250‑1 (MQ=255)                                                                                                                         gCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACt                                                                                                                       <  1:239662/251‑1 (MQ=255)                                                                                                                           gTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAAAAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTc                                                                                                                     >  1:319574/1‑251 (MQ=255)                                                                                                                                 ggCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccagcc                                                                                                                <  1:352052/250‑1 (MQ=255)                                                                                                                                       ggCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGt                                                                                                         >  1:203298/1‑251 (MQ=255)                                                                                                                                               gcgtctgaaGCAACTTGAGCGTCAACAAAGCTGGCTGAAGAAGAAACATATCGACTTACGTGTCTGCTGTGTTTCCAACGCGAAGGCACGGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGATCCGTTTAAGCTCGGGCGCTTAATTCGCCTCTTGAAAGAATATCATCTGCTGAACCCGGTCAGTGTTGACTGCACTTCCAGCCAGGCATTGGCGGATc                                                                                                 <  2:371715/243‑1 (MQ=255)                                                                                                                                                                aGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTcttg                                                                                <  1:203420/251‑4 (MQ=255)                                                                                                                                                                              tggctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAcc                                                                                                                                           <  1:103482/178‑1 (MQ=255)                                                                                                                                                                              tggctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAcc                                                                                                                                           >  2:103482/1‑178 (MQ=255)                                                                                                                                                                               ggctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGGTGCCAACTCGAAGGCACTGCGCCCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCAc                                                                 <  2:246174/251‑1 (MQ=255)                                                                                                                                                                                         tAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGATGGTTTCCACGTTGTCACGc                                                       >  1:50645/1‑251 (MQ=255)                                                                                                                                                                                          aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCgg                                                                                                                                                                                    <  1:384194/125‑1 (MQ=255)                                                                                                                                                                                          aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCgg                                                                                                                                                                                    >  2:384194/1‑125 (MQ=255)                                                                                                                                                                                          aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCaa                                                                                                                                                                                                      >  1:328721/1‑107 (MQ=255)                                                                                                                                                                                          aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCaa                                                                                                                                                                                                      <  2:328721/107‑1 (MQ=255)                                                                                                                                                                                             cATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGa                                                    >  1:2116/1‑250 (MQ=255)                                                                                                                                                                                               tatCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGTTCACCAATGTGAATGGCCTAAATCTGGAACACTGGAAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCTTGCGCGGAGGGTTCCCCGTTTTCACGCCGAAcc                                                 >  1:130967/1‑250 (MQ=255)                                                                                                                                                                                                                  ggTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGATCTGGCGCAAGCCAAAGAGCCGTTTAAGCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGTCTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCATCACCtcgtc                              <  1:252619/251‑1 (MQ=255)                                                                                                                                                                                                                        gCCAACTCGAAGGGACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATGCGCCTCGTGAAAGAATATCATCTTCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCATTGGCGGATCAATATGCCTACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCATCACCTCGTCGATGGa                        <  1:184039/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggcaatgcgcACCAATGTGTAAGGCCTAAATCTGGAAAACTGGCTGGAAGAACGTGCGTAAGCCAAAGAGCCGTGTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATTATCTGCTGAACCCGGTCATTGTTGACTGCACTCCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTAGCGCTAAGGTTTCCACGTTGTCACGCCGAACAAAAAGTCAAACGCCTCGGCGATCGATGACTACCATc             <  2:165951/242‑1 (MQ=255)                                                                                                                                                                                                                                       ctGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGATTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCCACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATctgat         >  1:67864/1‑247 (MQ=255)                                                                                                                                                                                                                                         gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt         <  1:376086/249‑1 (MQ=255)                                                                                                                                                                                                                                         gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTGCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt         <  1:171287/249‑1 (MQ=255)                                                                                                                                                                                                                                              ccAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGtcat  <  1:159096/251‑4 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    AGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTAT  >  NC_000913/1631‑2118

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