Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,258,484	T→C	Q389R (CAA→CGA)	tdcG ←	L‑serine deaminase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,258,484	0	T	C	100.0%	169.4 / NA	47	Q389R (CAA→CGA)	tdcG	L‑serine deaminase III
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/15);  total (32/15)

GGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGT  >  NC_000913/3258243‑3258727                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    ggTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTc                                                                                                                                                                                                                                            <  1:130740/251‑1 (MQ=255) ggTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTc                                                                                                                                                                                                                                            <  2:271842/251‑1 (MQ=255)     cacaTTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGc                                                                                                                                                                                                                                        <  1:117645/251‑1 (MQ=255)     cacaTTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGc                                                                                                                                                                                                                                        <  1:117653/251‑1 (MQ=255)       caTTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCa                                                                                                                                                                                                                                       >  1:317724/1‑250 (MQ=255)           tgtgCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCTGCAACCg                                                                                                                                                                                                                                   <  2:31299/250‑1 (MQ=255)            gtgCACCCAAGGATAAAAGCTGTCAGCAATGTCCGCCTCAGCCCACTTTAATGGCCAGACCTCCGCGTGATGTTTCGCGGTATTTATCGTCCCTATCTTTGTCGGTTTCATACATCGTCTCGATCACTTTCTCGAGTGTAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCACCCGGCGTTTACTGCTTTCACGGCCTTACTGGCATTACGTTCATTGCACGGGATTCGTACCTGTcagcacaccgt                                                                                                                                                                                                                                  >  2:170180/1‑240 (MQ=255)                aCCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGc                                                                                                                                                                                                                             >  2:108417/1‑251 (MQ=255)                     aGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTc                                                                                                                                                                                                                        >  1:259349/1‑251 (MQ=255)                        aTGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGc                                                                                                                                                                                                                     >  2:252053/1‑251 (MQ=255)                            aaGCTGACAGCAATGTCAGCCGCAGACCCCTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCa                                                                                                                                                                                                                  >  1:316633/1‑250 (MQ=255)                            aaGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCGTATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCa                                                                                                                                                                                                                  >  2:176519/1‑250 (MQ=255)                                      cAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTc                                                                                                                                                                                                       >  1:46000/1‑251 (MQ=255)                                       aaTGTCAGCCGCAGAACACTTTAATGGCCAGGCCTCAGCGTGATGTTTCGCGGTATTTATCATTCATATCTTTACCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGAACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTcc                                                                                                                                                                                                      >  1:342697/1‑251 (MQ=255)                                                   aGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCACGGTCTTTATCGTTCATATCTTTGCCGGGTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGTCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTAAGCCCAAGGTTATGCTCCATCGCGATTTcc                                                                                                                                                                                          >  1:162746/1‑251 (MQ=255)                                                    gACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccg                                                                                                                                                                                         <  2:208415/251‑1 (MQ=255)                                                     aCCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                                                                                                                                                                        >  1:303660/1‑251 (MQ=255)                                                                   cAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCt                                                                                                                                                                          <  2:192212/251‑1 (MQ=255)                                                                          ccGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCgg                                                                                                                                                                   <  1:10368/251‑1 (MQ=255)                                                                          ccGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCg                                                                                                                                                                    <  1:16235/250‑1 (MQ=255)                                                                           cgcgTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCgg                                                                                                                                                                   >  1:200162/1‑250 (MQ=255)                                                                              gTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACtt                                                                                                                                                               >  2:341884/1‑250 (MQ=255)                                                                              gTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGAATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCACGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGCCTg                                                                                                                                                               >  2:289217/1‑251 (MQ=255)                                                                                                            cTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc                                                                                                                                 >  1:137048/1‑251 (MQ=255)                                                                                                              ttGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccg                                                                                                                                <  2:282679/250‑1 (MQ=255)                                                                                                                  cGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATCCCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAtt                                                                                                                           <  2:32522/251‑1 (MQ=255)                                                                                                                     tttCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGTTGCCGAGGTGAGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTACTGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCACCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCAGCCCAGTCGTTCAGTTAACCCTGCCGCCGCCATTGa                                                                                                                         >  2:160379/1‑250 (MQ=255)                                                                                                                      ttCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAAc                                                                                                                       >  2:188542/1‑251 (MQ=255)                                                                                                                                 tctcGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGcc                                                                                                             >  1:53795/1‑250 (MQ=255)                                                                                                                                                              ggTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGcc                                                                               <  1:192746/251‑1 (MQ=255)                                                                                                                                                                   cGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCag                                                                             >  1:143253/1‑248 (MQ=255)                                                                                                                                                                    gAGGGGCGGAGCATCGCCATCGGCGGGGAGTGTACTGCTGTCACGGCATTAATGGAATTACGTTCAATGCACGGGATTCGTCCCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGCCAGCCGACTTCCGCGCCAGAGAt                                                                         <  2:157348/251‑1 (MQ=255)                                                                                                                                                                             cgcATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGtt                                                                <  1:28307/251‑1 (MQ=255)                                                                                                                                                                              gcATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTc                                                               >  1:290711/1‑251 (MQ=255)                                                                                                                                                                                            gcggcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAgcgc                                                 >  2:749/1‑251 (MQ=255)                                                                                                                                                                                               gcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGcc                                              <  1:58383/251‑1 (MQ=255)                                                                                                                                                                                               gcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGcc                                              >  1:336104/1‑251 (MQ=255)                                                                                                                                                                                                 gTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTACGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCACGCCa                                             >  1:144559/1‑250 (MQ=255)                                                                                                                                                                                                    tACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAAt                                           >  2:270258/1‑249 (MQ=255)                                                                                                                                                                                                          tttCACGGCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGAGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCcg                                   >  1:248085/1‑251 (MQ=255)                                                                                                                                                                                                                ggCATTAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCGCGGcc                              >  2:268659/1‑250 (MQ=255)                                                                                                                                                                                                                  cATTAATGGCATTACGTTCAATGCACGGGATTCGCATCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCCTTTTATACAGCGCGCCAATAGCCCCCGCGGCCAgc                           >  2:137463/1‑250 (MQ=255)                                                                                                                                                                                                                     tAATGGCATTACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCGCGGCCAGaaa                         >  1:170165/1‑250 (MQ=255)                                                                                                                                                                                                                             ttACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCCTTTTATACAGCGCGCCAATAGCCCCCGCGGCCAGAAAATAGCGGGc                >  1:198415/1‑251 (MQ=255)                                                                                                                                                                                                                             ttACGTTCAATGCACGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCGCGGCCAGAAAATAGCGGGc                >  2:120260/1‑251 (MQ=255)                                                                                                                                                                                                                                          aCGGGATTCGCACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCg   >  1:57616/1‑251 (MQ=255)                                                                                                                                                                                                                                           cGGGATTCGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGt  >  1:201935/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    GGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATACAGCGCGCCAATAGCCCCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGT  >  NC_000913/3258243‑3258727

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