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%	135.0 / NA	39	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 (22/17);  total (22/17)

ATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACCGG  >  NC_000913/4434890‑4435365                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  atatGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACtg                                                                                                                                                                                                                                   >  1:374249/1‑251 (MQ=255)                   ccGGAAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAAGCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGACGAGTGACTGTGGTTTCGTTCTGTTggg                                                                                                                                                                                                                 <  2:132026/251‑1 (MQ=255)                          aTTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTTAGGTTTCGTGCTGTTGGGATCAATc                                                                                                                                                                                                          >  1:344090/1‑251 (MQ=255)                                cgccgTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTa                                                                                                                                                                                                    <  2:374249/251‑1 (MQ=255)                                   cgTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAAc                                                                                                                                                                                                 >  2:58234/1‑251 (MQ=255)                                      aaGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAATGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAAGCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTg                                                                                                                                                                                               <  1:370608/250‑1 (MQ=255)                                            gATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGCTAATTCACACCATCAATAACATCAACGTTATAAGTGCGGAAACCATCCTAGTTGCTGAGTGACTTTGGTTTCGTGCt                                                                                                                                                                                                                       >  2:210579/1‑220 (MQ=255)                                            gATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCt                                                                                                                                                                                                                       <  1:210579/220‑1 (MQ=255)                                               aGTTATTGGTGGCACCCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAg                                                                                                                                                                                     <  2:26701/251‑1 (MQ=255)                                               aGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGAACTCGCCGTCATAACGGGCGGGCTGGGTAAAATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACAATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCgg                                                                                                                                                                                       >  1:220627/1‑249 (MQ=255)                                                       gtggCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTcc                                                                                                                                                                               >  1:292165/1‑249 (MQ=255)                                                            aaCCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCAct                                                                                                                                                                        >  2:75031/1‑251 (MQ=255)                                                             aCCAGGAACATGGCGTTCGGATTAATCAGCTTGCCCTTAAAGCTCAGGTTCTTACTCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACCTCAATCTGATAATTCACCCCATCAATCACATCCAAGTTATAAGTGCGGAAACCATCCTAGTTGATGCGTGACTGTGGTTTCGTGCTGTTGGTATCCATCTGGTTATACTGTCCCGCGGAGCCCTCCAGCCActc                                                                                                                                                                       >  2:370371/1‑251 (MQ=255)                                                                aGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCttt                                                                                                                                                                    >  1:44725/1‑251 (MQ=255)                                                                 ggAACATGGCGTTCGGATCAATCGGCTTGCCCTTAAAGATCAGGTTCTTCATCCTCTCCGCATTGGCATTACTCATCTGGCAACCGCCGTCATACAGGGCGGACTCAGTAACATCAATCTGCTAATTATCACCCTCAATAACTTCAAAGTTATAAATGCGGACACCATCCTCGTTGATGAGTGCCTGTGGTTTCGAGCTGTTGGGATCAATCTGGATAAACTGTCCCACGGAGCACTCGAGCCACTCTTTc                                                                                                                                                                   >  2:167540/1‑251 (MQ=255)                                                                  gAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCa                                                                                                                                                                  <  1:75031/251‑1 (MQ=255)                                                                    aCATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGACGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCAcc                                                                                                                                                                >  1:330860/1‑251 (MQ=255)                                                                           gTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCttt                                                                                                                                                           >  2:147670/1‑249 (MQ=255)                                                                           gTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCttt                                                                                                                                                           <  1:147670/249‑1 (MQ=255)                                                                               ggATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCt                                                                                                                                                     >  2:25230/1‑251 (MQ=255)                                                                                               cATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGc                                                                                                                                                                            >  1:65600/1‑212 (MQ=255)                                                                                               cATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGc                                                                                                                                                                            <  2:65600/212‑1 (MQ=255)                                                                                                              tCTTAAACATTTCCGCATTGGCATTAATCCTCTGGCACTCGCAGTCATAACGTGCGGGCTGGGTAACATCAATCTGATAATTCACACCCTCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTCGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCACGCGGAGCACTCCAGCAACTCTTTCACCTCTTTACCGCTGGATTTCACCACAATCAGCGAATTGGGAt                                                                                                                        >  2:193143/1‑249 (MQ=255)                                                                                                                ttAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATag                                                                                                                      >  1:133258/1‑249 (MQ=255)                                                                                                                          ccGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATcg                                                                                                          <  1:53730/251‑1 (MQ=255)                                                                                                                             cATTGGCATTAATCATCTGGCACTCGCTGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAATCCATCCTATTTGATGAGTGACTGTGGTTTCGGGCTGTTGGGATCAATCTGGTTAAACTGGCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATcggcg                                                                                                       <  2:404100/251‑1 (MQ=255)                                                                                                                                                        cGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCtttt                                                                            >  2:128943/1‑251 (MQ=255)                                                                                                                                                            aTAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCa                                                                        <  1:19656/251‑1 (MQ=255)                                                                                                                                                             tAACGGGCGGGCTGGGTAACATCACACTGATTATTCACACCATCAATAACATCAATGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTAGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGAACTCCAGCCACTCTTTCACCTCTTTACCGCTGGATTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGCAGGGAAAATGGCCTTTTTCCAc                                                                       >  2:338674/1‑251 (MQ=255)                                                                                                                                                                  ggcgggcTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCACGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGTGCTAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTcc                                                                   >  1:294342/1‑250 (MQ=255)                                                                                                                                                                                     aTCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCAtt                                               >  2:220029/1‑251 (MQ=255)                                                                                                                                                                                          aTAATTCACACCATCAAGAACATCAAAGTTATAAGTTCGGAACCCATCCTAGTTGATGAGTGACTGTTGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCTACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGc                                          <  1:340176/251‑1 (MQ=255)                                                                                                                                                                                              ttCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGacc                                      >  1:61150/1‑251 (MQ=255)                                                                                                                                                                                                     cATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCGTTTTCCACCTCCACATAGCTTGCCTGGGCATTTTTGCGACCACCGACt                               <  2:292165/251‑1 (MQ=255)                                                                                                                                                                                                      aTCAATAACATCAAAGTTATAAGTGAGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACtt                              <  2:238983/251‑1 (MQ=255)                                                                                                                                                                                                        cAATAACATCAAAGTTATAAGTGCGGAAACCAACCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATGTGGGTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGAGCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACttt                             <  1:392302/250‑1 (MQ=255)                                                                                                                                                                                                                              gctgAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCTCCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTTGCCTGTGTCCACCTCCACATAGCTGGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGc       <  1:237746/247‑1 (MQ=255)                                                                                                                                                                                                                                ggAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGCCTTTAAAAGGTGCGGCAGCTGAAAGCAcc    >  1:187523/1‑251 (MQ=255)                                                                                                                                                                                                                                  aaaCCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACCgg  <  1:25230/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  ATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACCGG  >  NC_000913/4434890‑4435365

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