Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,124	C→T	W481* (TGG→TAG)	cpdB ←	2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,124	0	C	T	100.0%	241.0 / NA	68	W481* (TGG→TAG)	cpdB	2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/38);  total (30/38)

TGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACC  >  NC_000913/4434897‑4435363                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                tGTCGCCCGTAACGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATAGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCAATCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTc                                                                                                                                                                                                                          >  2:473779/1‑251 (MQ=255)    cccgcgttccGGCAACTTTGCCGCCGTAAGCGCGATAGTGTTTGGTGTCAACCAGGACCATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCTCATGTGCATTAATCATCTGGCACTCGCCGGCATCACGGTCGGGCTTGGTATCATCAATCTGACAACTCACACCATCAAGAACATCTGCGTGATAATTGAGGAAACCAGCCTAGTTGATGAGTGGCGGTGTTTTCGtt                                                                                                                                                                                                                       <  2:659576/243‑2 (MQ=255)     gCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGc                                                                                                                                                                                                                      <  2:631385/251‑1 (MQ=255)           aCCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTg                                                                                                                                                                                                                 <  1:731214/250‑1 (MQ=255)           aCCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGGTGATGATTGACTGTGGTTTCGTGCTGTTgg                                                                                                                                                                                                                <  2:631702/251‑1 (MQ=255)            ccGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTggg                                                                                                                                                                                                               >  2:616848/1‑251 (MQ=255)                   aTTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATc                                                                                                                                                                                                        >  1:328106/1‑251 (MQ=255)                        ccgccgTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCACTAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGtt                                                                                                                                                                                                   >  1:553350/1‑251 (MQ=255)                          gccgTACGCGCGAGAGTCATTGGTGGCAACAAGGAACCTGGAGTTCGGATCAATCGGCTTGCCATTCAAGGCCAGGTCCTTACTCCTTTACGCATTGGCATAACTCATCTGGCCCTCGCCGTCATAACGGTCGGGCTGGGTAACATCAATCTGATAATTCATACCATCAATAACATCAAAGTTATCAGTGCGGAACCCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATATGGTTa                                                                                                                                                                                                  >  2:657073/1‑250 (MQ=255)                           ccgTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGGCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTaaa                                                                                                                                                                                                <  1:407562/251‑1 (MQ=255)                             gTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACt                                                                                                                                                                                              <  1:159869/251‑1 (MQ=255)                               aaGCGCGATAGTTAGTGGTGGCAACCAGGAACATGGCGGTCGGATCAATCGGCTTGCCATTAAAGGGCAGGTTCTTAATCCTTTCCGCAGTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAAGTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTTTGGTTTCGTGCTGTTGGGATCAATCTGTTTAAAcggt                                                                                                                                                                                            <  2:694811/251‑4 (MQ=255)                                   gcgATAGTTAGTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTTGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCcg                                                                                                                                                                                        <  1:370055/251‑1 (MQ=255)                                          ttATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGc                                                                                                                                                                                  >  1:162729/1‑250 (MQ=255)                                            aTTGGTGGCAAACAGTAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAACCTGTCCCGCGGAGCACt                                                                                                                                                                               <  2:53746/251‑1 (MQ=255)                                               ggtggCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCAGCCTAGTTGATGAGTGACTGTGGTGTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTcc                                                                                                                                                                             <  2:55215/250‑1 (MQ=255)                                                   gCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCa                                                                                                                                                                        <  2:328106/251‑1 (MQ=255)                                                        cAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCtt                                                                                                                                                                   >  1:55135/1‑251 (MQ=255)                                                         aGGAACATGGCGTTCGGATCAATCGGCTTGCCATTACAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCttt                                                                                                                                                                  >  2:67288/1‑251 (MQ=255)                                                               agggCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCGTAATCCTTGCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTAGATGAGTGACTGGGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACctc                                                                                                                                                            <  2:608749/249‑1 (MQ=255)                                                                 ggCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACctct                                                                                                                                                           <  1:480902/250‑1 (MQ=255)                                                                   cGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTa                                                                                                                                                        <  2:127115/251‑1 (MQ=255)                                                                        ggATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCt                                                                                                                                                   >  2:265035/1‑251 (MQ=255)                                                                         gATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTg                                                                                                                                                  <  1:228270/251‑1 (MQ=255)                                                                            cAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCt                                                                                                                                               <  2:231634/251‑1 (MQ=255)                                                                            aaagcgacTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGTCGCTCGCCGTCATAACGGGCGGGCTGGGTAACAGCAATCAGATAATTCACACCAGCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGTCTGGGTTTTCCGGCTTTTGGGATCAATCTGTCTAAACTGCCCCGCGGAGCGCTGCAGCCTCTGTTTCACCTCTGTACCGCTGGCt                                                                                                                                               <  2:671899/244‑1 (MQ=255)                                                                              aTCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCttt                                                                                                                                             <  2:157689/251‑1 (MQ=255)                                                                                      gCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCAGAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGAGGAGTGACTGTGGTTTCGTGCTGTTGGGAGCAATCTGGTTAAACTGTCCCGCGGAGCGCCCCAGCCACTCTTGCACCTCTTGACCGCTGGCTTTCACcaca                                                                                                                                      <  2:379845/250‑1 (MQ=255)                                                                                       ccATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGAGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGAGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAAt                                                                                                                                    >  1:352777/1‑251 (MQ=255)                                                                                         aTTAAAGGTCATGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAAGTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGGTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCa                                                                                                                                  <  2:33942/251‑1 (MQ=255)                                                                                             aaGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTTTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGt                                                                                                                              <  1:439681/251‑1 (MQ=255)                                                                                                gTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTAtt                                                                                                                           >  2:353558/1‑251 (MQ=255)                                                                                                       tCTTAAGCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGTAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATa                                                                                                                     <  2:671741/250‑1 (MQ=255)                                                                                                        cTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGAt                                                                                                                      >  1:339648/1‑248 (MQ=255)                                                                                                          tAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATag                                                                                                                    <  1:234003/248‑1 (MQ=255)                                                                                                                  cccgtaTTGACATTAATCATCTGGCATTCGCCGTCATAACTGGCGGGCTGGGAAACATCAATCTGATATGTCACACCATCAATAACATCAAAGTGATAAGTGCGGAAACCAGCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTGCACCACAATCAGCGTATTGGTATAGAGATAAAGATc                                                                                                         <  2:266922/246‑1 (MQ=255)                                                                                                                      cATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATcggcg                                                                                                     >  2:30046/1‑251 (MQ=255)                                                                                                                      cATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATcggcg                                                                                                     <  1:473779/251‑1 (MQ=255)                                                                                                                       aTTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGc                                                                                                                                >  2:608856/1‑223 (MQ=255)                                                                                                                       aTTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGc                                                                                                                                <  1:608856/223‑1 (MQ=255)                                                                                                                        ttGGCATCAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCCATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCggcggc                                                                                                   >  2:653506/1‑251 (MQ=255)                                                                                                                           gCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCAtt                                                                                                <  1:67288/251‑1 (MQ=255)                                                                                                                                          cACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTCTTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTg                                                                                 >  1:26793/1‑251 (MQ=255)                                                                                                                                           aCTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTATGGAAGGTCAACTgg                                                                                >  1:556830/1‑251 (MQ=255)                                                                                                                                              cGCCTTCATAACGGGCTGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCt                                                                             <  1:442608/251‑1 (MQ=255)                                                                                                                                                ccGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGCGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCttt                                                                           >  1:112715/1‑251 (MQ=255)                                                                                                                                                  gTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACAATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTAACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGACttttt                                                                         <  2:30256/251‑1 (MQ=255)                                                                                                                                                   tCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTc                                                                        >  1:646410/1‑251 (MQ=255)                                                                                                                                                               ggcTGGGTAACATCAACCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCAGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCGTTCACCACAATCACCGTATTGTGATAGAGATAAAGATCGGCGGCATTACGGAAGGACAACTGGCCTTTTTCCACCTCCCCATa                                                            <  2:248645/251‑1 (MQ=255)                                                                                                                                                                   gggTAACATCAATCTGATAATTCACACCAACAATAACATCAAAGTTAGAAGTGCGGAAACCATCCTAGGTGATGAGTGACTGTGGGTTCGTGCTGTTGGGATCAATCTGGTTAAACGGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCGCCACATAGCtt                                                        <  2:112715/251‑1 (MQ=255)                                                                                                                                                                         cATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCggg                                                  <  1:616848/251‑1 (MQ=255)                                                                                                                                                                                  gATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCAtttt                                           <  1:353558/249‑1 (MQ=255)                                                                                                                                                                                    tAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCg                                       >  2:673156/1‑251 (MQ=255)                                                                                                                                                                                            acCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGGCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTGTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGa                               <  1:653506/251‑1 (MQ=255)                                                                                                                                                                                                tCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACttt                           >  1:126059/1‑251 (MQ=255)                                                                                                                                                                                                 cAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTTCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACttt                           <  2:120971/250‑1 (MQ=255)                                                                                                                                                                                                      aCATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACgg                     >  2:352900/1‑251 (MQ=255)                                                                                                                                                                                                        aTCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCAttttt                                          >  2:183226/1‑228 (MQ=255)                                                                                                                                                                                                        aTCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCAttttt                                          <  1:183226/228‑1 (MQ=255)                                                                                                                                                                                                        aTCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTAGTGCTGTTGGGATCAATCTGCTTAAACTGTCCCACGGAGCACTCCAGCAACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCCGCGTATTGGGATAGAGATAAAGATCGGCGGCCTCACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTAGCCGGGTCATTTTTGCGACCACCGACTTTCAACGGTg                   >  1:737327/1‑251 (MQ=255)                                                                                                                                                                                                          cAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCg                 >  1:359382/1‑251 (MQ=255)                                                                                                                                                                                                                  taAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTg          >  2:14578/1‑250 (MQ=255)                                                                                                                                                                                                                     gTGCGGACACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCACTCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATagag                                                                                                                  >  2:97635/1‑143 (MQ=255)                                                                                                                                                                                                                     gTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATagag                                                                                                                  <  1:97635/143‑1 (MQ=255)                                                                                                                                                                                                                      ggCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAAGCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGc     <  1:265035/250‑1 (MQ=255)                                                                                                                                                                                                                        cGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCACTCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCAc   >  1:303101/1‑251 (MQ=255)                                                                                                                                                                                                                        cGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGCGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCAACCTCCACATAGCTTGCCGGGTCATTTTTGCGCCCACCGCCTTTAAACGGTGCGGCAGCTGAAAGCAc   >  2:625091/1‑251 (MQ=255)                                                                                                                                                                                                                         ggAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGACGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCAcc  >  1:563137/1‑251 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                TGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACC  >  NC_000913/4434897‑4435363

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