Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,042,663	A→G	R516R (CGT→CGC)	etk ←	tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,042,663	0	A	G	100.0%	226.2 / NA	64	R516R (CGT→CGC)	etk	tyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/36);  total (28/36)

TAAATAACACTTTTTGATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGC  >  NC_000913/1042515‑1042809                                                                                                                                                     |                                                                                                                                                   taaataaCACTTTTTGATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGcgcgcg                                                                                                                                                  >  2:691213/1‑151 (MQ=255)    ataaCACTTTTTGATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTAc                                                                                                                                               <  2:80754/151‑1 (MQ=255)        cacTTTTTGATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCt                                                                                                                                           >  2:884881/1‑151 (MQ=255)         acTTTTTGATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCtt                                                                                                                                          >  1:389029/1‑151 (MQ=255)                gATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCAcagc                                                                                                                                   >  1:81567/1‑151 (MQ=255)                  tCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCAcagc                                                                                                                                   >  2:282528/1‑149 (MQ=255)                    ggACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGa                                                                                                                                >  2:318293/1‑150 (MQ=255)                      aCTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATc                                                                                                                             <  1:469881/151‑1 (MQ=255)                            cGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCgg                                                                                                                       >  1:653273/1‑151 (MQ=255)                             gATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCggg                                                                                                                      >  2:726169/1‑151 (MQ=255)                             gATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCggg                                                                                                                      >  1:304430/1‑151 (MQ=255)                                  ccGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTAt                                                                                                                  <  1:857009/150‑1 (MQ=255)                                   cGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATcc                                                                                                                <  2:797534/151‑1 (MQ=255)                                         cAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGcc                                                                                                          >  1:895960/1‑151 (MQ=255)                                           gagTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCa                                                                                                         <  2:93386/150‑1 (MQ=255)                                           gagTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAg                                                                                                        <  1:336955/151‑1 (MQ=255)                                           gagTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAg                                                                                                        <  2:599015/151‑1 (MQ=255)                                             gTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAgg                                                                                                       >  2:533568/1‑150 (MQ=255)                                             gTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGa                                                                                                      <  1:689341/151‑1 (MQ=255)                                               tGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGaa                                                                                                     <  2:714102/150‑1 (MQ=255)                                                 aaCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAggg                                                                                                  >  2:782526/1‑151 (MQ=255)                                                     gACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGAATGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAAGCCGCCGGGTTATCCACCGCCAGGAAGGGGaga                                                                                              <  2:410095/151‑3 (MQ=255)                                                         aaCGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAtttt                                                                                          >  1:441450/1‑151 (MQ=255)                                                           cGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAttttt                                                                                         >  2:391322/1‑150 (MQ=255)                                                           cGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTa                                                                                        <  1:18776/151‑1 (MQ=255)                                                            gTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAg                                                                                       <  2:80699/151‑1 (MQ=255)                                                               ttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTa                                                                                     <  2:402713/150‑1 (MQ=255)                                                                       gTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGgc                                                                              <  2:465525/149‑1 (MQ=255)                                                                       gTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGgc                                                                              <  2:183705/149‑1 (MQ=255)                                                                         cTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgctg                                                                          >  2:441526/1‑151 (MQ=255)                                                                           ggCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgc                                                                         <  1:793470/150‑1 (MQ=255)                                                                            gCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGctgctg                                                                       <  2:13056/151‑1 (MQ=255)                                                                             cGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGa                                                                      <  1:745695/151‑1 (MQ=255)                                                                              gTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGAt                                                                     >  2:322782/1‑151 (MQ=255)                                                                                cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTa                                                                   <  2:816039/151‑1 (MQ=255)                                                                                     cGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaaa                                                               >  1:646760/1‑150 (MQ=255)                                                                                     cGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaaa                                                               >  2:667207/1‑150 (MQ=255)                                                                                         gATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataa                                                           <  2:698806/150‑1 (MQ=255)                                                                                          aTCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataaat                                                         >  1:253388/1‑151 (MQ=255)                                                                                          aTCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataaat                                                         <  2:683821/151‑1 (MQ=255)                                                                                          aTCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataaat                                                         <  2:669240/151‑1 (MQ=255)                                                                                            cATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                       >  2:262164/1‑151 (MQ=255)                                                                                             aTTAGAATGTTATGCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGTAGAATCCGCCGGGTTATCCTCCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                       <  2:228987/150‑1 (MQ=255)                                                                                               tAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttttt                                                     <  2:413285/150‑1 (MQ=255)                                                                                                aGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCt                                                   >  1:224610/1‑151 (MQ=255)                                                                                                  aaTGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                 <  1:832249/151‑1 (MQ=255)                                                                                                       tATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCa                                             <  1:214966/150‑1 (MQ=255)                                                                                                       tATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCa                                             <  1:260959/150‑1 (MQ=255)                                                                                                        aTTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAg                                            <  1:698226/150‑1 (MQ=255)                                                                                                             ccGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACgg                                        <  2:826776/149‑1 (MQ=255)                                                                                                             ccGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACg                                         <  2:645385/148‑1 (MQ=255)                                                                                                               gTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGc                                    >  2:576326/1‑151 (MQ=255)                                                                                                                   ccATCATAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGt                                  <  2:734733/149‑1 (MQ=255)                                                                                                                     atcatAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGa                                           <  1:107896/138‑1 (MQ=255)                                                                                                                     atcatAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGa                                           >  2:107874/1‑138 (MQ=255)                                                                                                                        atAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATcc                            >  1:404998/1‑150 (MQ=255)                                                                                                                        atAGCGAAATGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCa                           <  1:867653/151‑1 (MQ=255)                                                                                                                              aaaTGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACt                     <  1:235993/151‑1 (MQ=255)                                                                                                                               aaTGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTc                    >  1:445881/1‑151 (MQ=255)                                                                                                                                 tGCAGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCgg                  >  2:183055/1‑151 (MQ=255)                                                                                                                                    aGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGa                 >  2:45976/1‑149 (MQ=255)                                                                                                                                    aGACTGGTTCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGAc                >  2:360399/1‑150 (MQ=255)                                                                                                                                           ttCGTAGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGa         <  2:206847/150‑1 (MQ=255)                                                                                                                                                aGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTgg   >  2:265333/1‑151 (MQ=255)                                                                                                                                                aGCGCGCGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTgg   <  1:183170/151‑1 (MQ=255)                                                                                                                                                 gcgcgcGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGc  <  2:585725/151‑1 (MQ=255)                                                                                                                                                 gcgcgcGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAAAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTgg   <  2:751506/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TAAATAACACTTTTTGATCGGACTGGGCGATCACCGCTGCCAGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGC  >  NC_000913/1042515‑1042809

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