Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,860	T→A	A508A (GCT→GCA)	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,860	0	T	A	94.7%	102.4 / ‑3.0	39	A508A (GCT→GCA)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (15/21);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

GTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGA  >  NC_000913/1618‑2096                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             gTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCAc                                                                                                                                                                                                                                      <  1:218887/251‑1 (MQ=255)     ccATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAAt                                                                                                                                                                                                                                  >  1:244003/1‑251 (MQ=255)      cATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGAGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGAGTTGGCGGTGCGCTGCTGGGGCAACTGTAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAAATAGAAGGCACTGCGCACCAATg                                                                                                                                                                                                                                 >  1:102613/1‑251 (MQ=255)              aGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGACGAATAAACATATCGACTTACGTGTATGCGGTGTTGCCAACTCGAATGCACAGCTCACCAATGTGCATGGc                                                                                                                                                                                                                         >  2:51367/1‑251 (MQ=255)               gggATCGTCTGAACTCTCAAGCCCGGGCGTGGTAAATAACGAGGATGCGACCCCTGGCGTGCGCGTTACGCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGAGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGGGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGTCTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTTACAAATGTGCATGGcc                                                                                                                                                                                                                        <  2:360943/251‑1 (MQ=255)                 gctcttctGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCAGCAGATGCTGTTCAATAACGATCAGGTTATCGCTGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGGGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCGGGCTGAAGAATAAACAAATCGACTTACTGGTCTGCGGTGTTGCCAAATCGAAGGCACTGCTCACCAATGtgacgggcct                                                                                                                                                                                                                       >  2:239662/3‑241 (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                                                                                                                                                                                <  2:328721/107‑1 (MQ=255)                                                                                                                                                                                                       aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCaa                                                                                                                                                                                >  1:328721/1‑107 (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)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             GTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGA  >  NC_000913/1618‑2096

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