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%	222.2 / NA	65	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 (40/25);  total (40/25)

GCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACCGGCAGTT  >  NC_000913/4434881‑4435370                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       gCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGGTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCACTAACATCAAAGTTATAAGTGCGGACACCATCCTAGTtgat                                                                                                                                                                                                                                                 >  1:58322/1‑251 (MQ=255)  cAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAAAGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCACTAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTtgatg                                                                                                                                                                                                                                                >  2:364477/1‑251 (MQ=255)         gATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCCCACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACt                                                                                                                                                                                                                                         >  1:381770/1‑251 (MQ=255)              ggCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGATAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTgg                                                                                                                                                                                                                                     <  1:164673/250‑1 (MQ=255)                          tACCGGCAAATTTGCCGCCGTACGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTg                                                                                                                                                                                                                        >  1:471559/1‑251 (MQ=255)                           cccGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTgg                                                                                                                                                                                                                       <  1:281118/250‑1 (MQ=255)                           aCCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTg                                                                                                                                                                                                                        <  1:390448/250‑1 (MQ=255)                                 aaaTTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCa                                                                                                                                                                                                                  >  1:189043/1‑250 (MQ=255)                                    tttGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGCCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGAGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGCTCAATCt                                                                                                                                                                                                              >  2:162115/1‑251 (MQ=255)                                                       tAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCgg                                                                                                                                                                                            >  1:110227/1‑250 (MQ=255)                                                        aGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCgg                                                                                                                                                                                            >  2:63320/1‑249 (MQ=255)                                                        aGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCgg                                                                                                                                                                                            >  1:352233/1‑249 (MQ=255)                                                                       ccAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCActct                                                                                                                                                                           >  1:113644/1‑251 (MQ=255)                                                                       ccAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCActct                                                                                                                                                                           <  2:113644/251‑1 (MQ=255)                                                                          ggAACAGGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTAAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAAGAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTGGGGATCAATCTGGTTAAACTGTCCCGCGTAGCACTCCAGCCACTCTTTc                                                                                                                                                                        <  2:350156/251‑1 (MQ=255)                                                                             aCATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCAcc                                                                                                                                                                     >  1:236230/1‑251 (MQ=255)                                                                             aCATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTCTCCGCATTGGCATTAATCATCTGGAACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATACTTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGCTGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTcccc                                                                                                                                                                     >  1:571181/1‑248 (MQ=255)                                                                                     ttCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCATCTCTTTAcc                                                                                                                                                             >  2:441581/1‑251 (MQ=255)                                                                                      ccGGATCAATCGGCTTGCCATTAAAGGTCAGGTGCTTAATCCTTGCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGGGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCg                                                                                                                                                            <  2:236230/250‑1 (MQ=255)                                                                                          aTCAATCGGCTTGCCATTAAAGGTCCGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTgg                                                                                                                                                        >  2:486009/1‑251 (MQ=255)                                                                                          aTCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACGCCAGCCACTCTTTCACCTCTTTACCGCTgg                                                                                                                                                        <  1:454492/251‑1 (MQ=255)                                                                                             attcgacttgcaattcacggcaaGGTTGTTAATCCTTACCGCATTGACATTAATCATATGGCACTAGCCATACTACCGGGCGGGCTGACTAAAATCAATCTGATAATACCCACCATCAATAACATCAAAGTTGTAGGTGCGGACACCATCCTAGTTGCTGCGTGACTGGGGATTCGTGCTGTTGGGATCAATCTGGTTAAACTGGCCCGCGGAGCACTCCACCCACTCTATGACCTCTTTACCCCTGCCtt                                                                                                                                                     >  2:299476/23‑251 (MQ=255)                                                                                              aTCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCttt                                                                                                                                                    <  1:562292/251‑1 (MQ=255)                                                                                                     tGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCATCCACTCTTTCACCTCTTTACCGCTGGCTTTcacca                                                                                                                                               <  1:194793/249‑1 (MQ=255)                                                                                                      gCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACcaca                                                                                                                                             >  2:101239/1‑250 (MQ=255)                                                                                                         aTTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCTTTAATCATCTGGCACTCGCCGTAATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAACCCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACaacca                                                                                                                                         <  2:24973/251‑4 (MQ=255)                                                                                                         aTTAAAGGTCAGGTTCTTAATCCGTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAATTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACCCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCa                                                                                                                                         <  2:352233/251‑1 (MQ=255)                                                                                                          ttAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATc                                                                                                                                          >  2:328625/1‑249 (MQ=255)                                                                                                           tAAAGGGCAGGTTCTTAATCCTTTCCGCATGGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGTATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACAACAATCAGc                                                                                                                                       >  2:212069/1‑251 (MQ=255)                                                                                                                 tCAGGATCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTg                                                                                                                                 >  2:523609/1‑251 (MQ=255)                                                                                                                   aGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTggg                                                                                                                               >  1:408378/1‑251 (MQ=255)                                                                                                                   aGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTggg                                                                                                                               <  1:162132/251‑1 (MQ=255)                                                                                                                    ggTTCTTAATCCTTTCCGCCTTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAACCCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTggg                                                                                                                               >  1:118762/1‑250 (MQ=255)                                                                                                                      ttCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGATATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGAt                                                                                                                             >  2:129440/1‑250 (MQ=255)                                                                                                                          tAATCCTTTCCGCATTGGCATTCATCATCTGGAACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTCTTGGGATag                                                                                                                           >  2:379764/1‑248 (MQ=255)                                                                                                                          tAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGc                                                                                                                                                       <  1:22246/220‑1 (MQ=255)                                                                                                                          tAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGc                                                                                                                                                       >  2:22246/1‑220 (MQ=255)                                                                                                                          tAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATag                                                                                                                           <  1:307438/248‑1 (MQ=255)                                                                                                                                      cATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATcggcg                                                                                                            >  1:245484/1‑251 (MQ=255)                                                                                                                                             aTTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTAc                                                                                                     >  2:6895/1‑251 (MQ=255)                                                                                                                                                        ggCACTCGACGTAATAACGGGCGGGCTGGGTAACAACAATCTGATAATTCACACAATCAATCACATCAATGATATAAGTGCGGAAAACATCCTAGTTGATACGTGACTGTGGTTTCGAGCTGTTGGGATCAATCTGGTAAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCCGCGTATTGGGATAGAGATAAACTTCGGCGCCATTACGGAAGATCacc                                                                                          >  2:490791/1‑249 (MQ=255)                                                                                                                                                          cACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTg                                                                                        >  2:196948/1‑251 (MQ=255)                                                                                                                                                           aCTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAACATGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTgg                                                                                       >  1:237574/1‑251 (MQ=255)                                                                                                                                                                 cGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGc                                                                                                                                       >  2:184989/1‑197 (MQ=255)                                                                                                                                                                 cGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGc                                                                                                                                       <  1:184989/197‑1 (MQ=255)                                                                                                                                                                        aCGGTCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTAAAAGTACGGAAACCATCCTAGTGTATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCGGGTTAAACGGTCCCCCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAAGCAGCGTATTGGGATAGAGACAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCt                                                                          <  2:18342/251‑1 (MQ=255)                                                                                                                                                                          gggcgggcTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTcc                                                                        <  1:63320/251‑1 (MQ=255)                                                                                                                                                                           ggcgggcTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGTCTTTTTCCACCTCca                                                                       <  1:374665/251‑1 (MQ=255)                                                                                                                                                                              gggcTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACAt                                                                    <  2:166954/251‑1 (MQ=255)                                                                                                                                                                               ggcTGGGAAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACAt                                                                    >  1:294517/1‑250 (MQ=255)                                                                                                                                                                               ggcTGAGTAACCTCAATCTGATACTTCACACCATCAATAACATCAAAGCTATAAGCGCGGAAACCATCCTAGTTGACGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTCCGGCAGGTCAACTGGCCTTTTTCCACCTCCACAt                                                                    >  2:40852/1‑250 (MQ=255)                                                                                                                                                                                         caacaatcgtaaaaTTCACACCATCAATAACATCAAAGTTATAATTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTTCCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCggg                                                         <  2:467943/239‑1 (MQ=255)                                                                                                                                                                                            caatcggaggaggcggacgggcggtAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCACTCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACAGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGCTCGGCGGCATTACGGACGGTCCACTGGCCTTTTTCCACCTCCACCTAGCTTGCCGGGTc                                                       >  1:224175/25‑250 (MQ=255)                                                                                                                                                                                            cAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTAATGAGTGACTATGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAATATCGGAGGCATTAAGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTc                                                       >  2:177873/1‑250 (MQ=255)                                                                                                                                                                                                  gATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGGGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCATCTGGCCTTTTTCCACCTCCTCATAGCTTGCCGGGTCAttgttg                                                <  2:118471/251‑1 (MQ=255)                                                                                                                                                                                                             ccATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGCCCACCgcc                                     >  2:106679/1‑249 (MQ=255)                                                                                                                                                                                                                   aTAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTaaa                               <  1:364477/251‑1 (MQ=255)                                                                                                                                                                                                                    aaaCATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCGGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAAc                              <  1:523609/250‑1 (MQ=255)                                                                                                                                                                                                                     aaCATCAAAGTTATAAGTGAGGAACACATCCTAGTTGATGCGTGACTGTGGTTTCGTGCTGTTGGCATCCATCTGGTTAAACTGTCCCGCGGACCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGCTCGAGATAAAGATCGGCGGCATTACGGAAGGTCCACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACg                             >  2:577464/1‑251 (MQ=255)                                                                                                                                                                                                                              gTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAg                    <  1:101239/251‑1 (MQ=255)                                                                                                                                                                                                                               ttATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGTCTTTAAACGGTGCGGCAGc                   >  2:145886/1‑251 (MQ=255)                                                                                                                                                                                                                                  taAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTg                 >  1:78348/1‑250 (MQ=255)                                                                                                                                                                                                                                         ggAAACCATCCTAGTTGCTGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGCGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCAACGACTTTAAACGGTGCGGCAGCTGAAAGCAcc         >  2:275892/1‑251 (MQ=255)                                                                                                                                                                                                                                         ggAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCAcc         <  1:145886/251‑1 (MQ=255)                                                                                                                                                                                                                                                aTCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACCGGCAGtt  >  2:186070/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       GCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACCGGCAGTT  >  NC_000913/4434881‑4435370

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: ‑