Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	450,073	C→T	W190* (TGG→TAG)	cyoB ←	cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	450,073	0	C	T	100.0%	259.4 / NA	65	W190* (TGG→TAG)	cyoB	cytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/30);  total (35/30)

GTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAGTTGTTGA  >  NC_000913/449930‑450220                                                                                                                                                |                                                                                                                                                    gTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAAt                                                                                                                                               >  2:142367/1‑150 (MQ=255)      ttCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGacac                                                                                                                                         <  1:671922/151‑1 (MQ=255)         gcgcACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGT‑ACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGa                                                                                                                                      >  1:327094/1‑150 (MQ=255)          cgcACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGa                                                                                                                                      >  1:191658/1‑150 (MQ=255)            cacaGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACt                                                                                                                                    <  2:670337/150‑1 (MQ=255)                gTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCgg                                                                                                                               >  1:252208/1‑151 (MQ=255)                gTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCgg                                                                                                                               >  1:252220/1‑151 (MQ=255)                   aTGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGAc                                                                                                                             <  2:478577/150‑1 (MQ=255)                   aTGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGAc                                                                                                                             <  2:359917/150‑1 (MQ=255)                   aTGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGAc                                                                                                                             <  1:124634/150‑1 (MQ=255)                   aTGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGAc                                                                                                                             <  1:331344/150‑1 (MQ=255)                    tGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTg                                                                                                                           <  2:590152/151‑1 (MQ=255)                     gCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGa                                                                                                                              >  1:363210/1‑147 (MQ=255)                       ccAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAc                                                                                                                        <  2:643142/151‑1 (MQ=255)                       ccAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAc                                                                                                                        <  2:673009/151‑1 (MQ=255)                       ccAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAc                                                                                                                        <  2:80012/151‑1 (MQ=255)                          ggTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCGGGACTGTAct                                                                                                                       >  1:471859/1‑149 (MQ=255)                          ggTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAct                                                                                                                       >  1:108286/1‑149 (MQ=255)                           gTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTActct                                                                                                                     <  2:511888/150‑1 (MQ=255)                           gTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTActct                                                                                                                     <  1:732188/150‑1 (MQ=255)                           gTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTActct                                                                                                                     <  1:312646/150‑1 (MQ=255)                           gTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTActc                                                                                                                      <  1:297029/149‑1 (MQ=255)                            tAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTAt                                                                                                                   <  2:607444/151‑1 (MQ=255)                            tAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTAt                                                                                                                   <  2:275208/151‑1 (MQ=255)                              aaTACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTc                                                                                                                 <  2:526423/151‑1 (MQ=255)                               aTACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTTGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTc                                                                                                                 >  1:135816/1‑150 (MQ=255)                               aTACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTcc                                                                                                                <  2:164522/151‑1 (MQ=255)                               aTACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTcc                                                                                                                <  1:91209/151‑1 (MQ=255)                               aTACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTc                                                                                                                 >  1:684629/1‑150 (MQ=255)                                tACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTcc                                                                                                                >  2:350230/1‑150 (MQ=255)                                 aCTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCg                                                                                                              <  2:250841/151‑1 (MQ=255)                                 aCTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCg                                                                                                              <  2:227043/151‑1 (MQ=255)                                    ggCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATa                                                                                                           >  2:427843/1‑151 (MQ=255)                                        tCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCgg                                                                                                       <  2:140468/151‑1 (MQ=255)                                          ttGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTg                                                                                                     <  2:719115/151‑1 (MQ=255)                                            gAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGa                                                                                                   >  2:40487/1‑151 (MQ=255)                                             aaCATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGAt                                                                                                  <  2:212306/151‑1 (MQ=255)                                                  ggTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgc                                                                                              >  1:553365/1‑150 (MQ=255)                                                   gTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgcca                                                                                            >  1:577726/1‑151 (MQ=255)                                                   gTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgcc                                                                                             <  1:126505/150‑1 (MQ=255)                                                   gTCATGCCCGGAGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgcca                                                                                            >  1:321722/1‑151 (MQ=255)                                                     cATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGACAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgccag                                                                                           >  1:190303/1‑150 (MQ=255)                                                                         tCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTTGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAAtt                                                                       >  1:614007/1‑150 (MQ=255)                                                                            gAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGc                                                                    >  1:576851/1‑150 (MQ=255)                                                                            gAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGc                                                                    >  2:425776/1‑150 (MQ=255)                                                                             aaTGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGc                                                                    <  2:519579/149‑1 (MQ=255)                                                                                       aagaagTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCgagag                                                        <  1:402880/151‑1 (MQ=255)                                                                                           agTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGaaa                                                     >  2:394104/1‑150 (MQ=255)                                                                                            gTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACg                                                   >  1:409171/1‑151 (MQ=255)                                                                                                    cGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAg                                           >  2:145464/1‑151 (MQ=255)                                                                                                       tAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAAt                                        <  2:619945/151‑1 (MQ=255)                                                                                                            gtcgtACTTATACTGGATAGCTGGAGACTCAATATCTAGTAATAGACACCGACTCCCGGACTGTACTCTATTCCCGAAAGCTGTGGTAATGCCAGCCTGCCGGTCAGCGCAAAATCGCCCACACCGAGAGAAACGTTAACCAGAATcaca                                    >  1:325846/1‑150 (MQ=255)                                                                                                             tcgtACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCacac                                   >  1:227440/1‑150 (MQ=255)                                                                                                             tcgtACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCacac                                   <  1:121659/150‑1 (MQ=255)                                                                                                                tACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACcaa                                >  1:167913/1‑150 (MQ=255)                                                                                                                tACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACcaa                                >  2:296342/1‑150 (MQ=255)                                                                                                                tACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACcaa                                >  1:195784/1‑150 (MQ=255)                                                                                                                tACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACcaa                                >  1:258723/1‑150 (MQ=255)                                                                                                                tACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACcaa                                >  1:288267/1‑150 (MQ=255)                                                                                                                       aCCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTa                        >  1:230683/1‑151 (MQ=255)                                                                                                                          ggATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTaaa                      >  2:536965/1‑150 (MQ=255)                                                                                                                            aTAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAAcc                    <  2:297902/150‑1 (MQ=255)                                                                                                                             tAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAg                  >  1:514259/1‑151 (MQ=255)                                                                                                                                       cTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAgt        >  1:180582/1‑151 (MQ=255)                                                                                                                                             atCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAGTTGTTGa  >  2:666806/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                    GTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAGTTGTTGA  >  NC_000913/449930‑450220

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