Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,440,378	G→T	R218S (CGT→AGT)	rpoA ←	RNA polymerase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,440,378	0	G	T	100.0%	281.1 / NA	71	R218S (CGT→AGT)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (33/38);  total (33/38)

AGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACGGC  >  NC_000913/3440232‑3440524                                                                                                                                                   |                                                                                                                                                   agagCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATc                                                                                                                                                >  2:834398/1‑151 (MQ=255)  gagCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCg                                                                                                                                               <  1:18703/151‑1 (MQ=255)  gagCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCg                                                                                                                                               <  2:111682/151‑1 (MQ=255)             cAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAgg                                                                                                                                     >  1:98501/1‑150 (MQ=255)                ttCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCg                                                                                                                                 >  2:954515/1‑151 (MQ=255)                 tCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCg                                                                                                                                 <  2:31646/150‑1 (MQ=255)                     gATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATtgt                                                                                                                            <  1:323350/151‑1 (MQ=255)                     gATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATtgt                                                                                                                            >  1:615397/1‑151 (MQ=255)                        cttCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGcc                                                                                                                         >  1:1038966/3‑151 (MQ=255)                        cGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGcc                                                                                                                         >  1:439159/1‑151 (MQ=255)                         gTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCg                                                                                                                        >  2:177948/1‑151 (MQ=255)                         gTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCg                                                                                                                        <  1:942399/151‑1 (MQ=255)                          tCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGt                                                                                                                       >  1:59634/1‑151 (MQ=255)                            aaCAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTg                                                                                                                     >  1:803628/1‑151 (MQ=255)                                gggCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGttt                                                                                                                 >  2:1011149/1‑151 (MQ=255)                                 ggCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTc                                                                                                                >  2:914119/1‑151 (MQ=255)                                  gcgcAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTcc                                                                                                               <  1:746943/151‑1 (MQ=255)                                      agcagGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCAttt                                                                                                           <  1:1019688/151‑1 (MQ=255)                                      agcagGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCAttt                                                                                                           <  1:1063460/151‑1 (MQ=255)                                       gcagGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTc                                                                                                          <  1:421518/151‑1 (MQ=255)                                        cagGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTc                                                                                                          <  2:84439/150‑1 (MQ=255)                                         agGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCg                                                                                                         >  2:456264/1‑150 (MQ=255)                                          ggatcggatcGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGAt                                                                                                       <  2:60059/151‑1 (MQ=255)                                           gatcggatcgAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATg                                                                                                      >  2:317648/1‑151 (MQ=255)                                              cggatcgAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGAcc                                                                                                   <  2:334031/151‑1 (MQ=255)                                                  tcgATCTCTAATTTCTCTTCTTTCACTTCAGGCTGACGTACATCCCGTAAGTCTACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGTTCGATTGCGCCGTTGGTTTCCATTTCGATGAcaa                                                                                                  >  2:120767/1‑146 (MQ=255)                                                   cgAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCt                                                                                               <  2:27091/150‑1 (MQ=255)                                                     aaCTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGt                                                                                            <  1:302207/151‑1 (MQ=255)                                                      aCTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTc                                                                                           >  2:514579/1‑151 (MQ=255)                                                        tctGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTcc                                                                                          >  2:514660/1‑150 (MQ=255)                                                           ggTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCa                                                                                         <  2:482353/148‑1 (MQ=255)                                                           ggTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCa                                                                                         <  1:297628/148‑1 (MQ=255)                                                           ggTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGt                                                                                      >  2:468140/1‑151 (MQ=255)                                                               tctcTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCgg                                                                                   <  2:412584/150‑1 (MQ=255)                                                                      tttCGCTTCAGGCGGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCCTGTCCGGGTCGGTACg                                                                               <  2:245990/147‑1 (MQ=255)                                                                      tttCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGt                                                                           >  2:1027748/1‑151 (MQ=255)                                                                      tttCACTTCAGGCTGACGCACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTcc                                                                                          <  2:462559/136‑1 (MQ=255)                                                                      tttCACTTCAGGCTGACGCACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTcc                                                                                          >  1:462561/1‑136 (MQ=255)                                                                         cACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCt                                                                        >  1:799796/1‑151 (MQ=255)                                                                              cAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTacac                                                                    <  2:538479/150‑1 (MQ=255)                                                                              cAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACAcg                                                                   >  2:632883/1‑151 (MQ=255)                                                                               aGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACAcgc                                                                  >  2:230132/1‑151 (MQ=255)                                                                               aGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACAcgc                                                                  >  1:663560/1‑151 (MQ=255)                                                                                 gCTGACGTACATCGCGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACgcgc                                                                <  2:349818/151‑1 (MQ=255)                                                                                    gACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCgctgc                                                             >  2:699646/1‑151 (MQ=255)                                                                                    gACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCgctgc                                                             <  2:98501/151‑1 (MQ=255)                                                                                      cGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCtt                                                           <  2:939829/151‑1 (MQ=255)                                                                                          cATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGCTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAAc                                                       >  2:895718/1‑151 (MQ=255)                                                                                             cACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACAtt                                                    <  2:594800/151‑1 (MQ=255)                                                                                                   aGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGc                                              <  2:1038952/151‑1 (MQ=255)                                                                                                   aGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGc                                              <  2:596716/151‑1 (MQ=255)                                                                                                       aaCGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATa                                          <  1:450149/151‑1 (MQ=255)                                                                                                          gAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGCTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGc                                       >  2:1061959/1‑151 (MQ=255)                                                                                                           aaaGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCt                                      <  2:495079/151‑1 (MQ=255)                                                                                                             aGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTcc                                    <  1:935894/151‑1 (MQ=255)                                                                                                             aGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTcc                                    >  2:792088/1‑151 (MQ=255)                                                                                                               cTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCcac                                  >  2:279343/1‑151 (MQ=255)                                                                                                                   cAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAggg                              <  1:894803/151‑1 (MQ=255)                                                                                                                     gttgttCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCt                            <  2:56749/151‑1 (MQ=255)                                                                                                                     gttgttCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCt                            >  1:544204/1‑151 (MQ=255)                                                                                                                        gttCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGt                          <  1:587840/150‑1 (MQ=255)                                                                                                                          tCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAg                        <  1:632329/150‑1 (MQ=255)                                                                                                                          tCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGc                       >  1:522743/1‑151 (MQ=255)                                                                                                                            aGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCAt                     <  2:426432/151‑1 (MQ=255)                                                                                                                               cAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCg                  <  2:295252/151‑1 (MQ=255)                                                                                                                                   aTGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCg               >  1:400478/1‑150 (MQ=255)                                                                                                                                   aTGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGa              <  2:824311/151‑1 (MQ=255)                                                                                                                                     ggTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGAcc            >  1:116717/1‑151 (MQ=255)                                                                                                                                      gTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACca           <  1:343606/151‑1 (MQ=255)                                                                                                                                            gCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGAc     >  2:123890/1‑151 (MQ=255)                                                                                                                                             cACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATATGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACg    <  1:177950/151‑1 (MQ=255)                                                                                                                                               cgactaATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACgg   <  2:653731/150‑1 (MQ=255)                                                                                                                                               cgactaATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACGGc  >  2:600266/1‑151 (MQ=255)                                                                                                                                               cgactaATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACGGc  <  1:492170/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACGGC  >  NC_000913/3440232‑3440524

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