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%	239.4 / NA	59	R218S (CGT→AGT)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (23/36);  total (23/36)

GCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACGGCC  >  NC_000913/3440235‑3440525                                                                                                                                                |                                                                                                                                                    gCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGcc                                                                                                                                              >  1:531096/1‑151 (MQ=255)    gACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCctct                                                                                                                                           <  2:692900/151‑1 (MQ=255)                 aGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATtg                                                                                                                              <  2:538306/151‑1 (MQ=255)                        cAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGtt                                                                                                                       >  1:527872/1‑151 (MQ=255)                         aaCAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTg                                                                                                                      <  2:362272/151‑1 (MQ=255)                              ggCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTc                                                                                                                 <  1:97415/151‑1 (MQ=255)                              ggCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTc                                                                                                                 <  2:250840/151‑1 (MQ=255)                                   agccgGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGCACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCTCGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCAttt                                                                                                            >  2:883358/1‑151 (MQ=255)                                   agcagGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCAttt                                                                                                            >  1:377714/1‑151 (MQ=255)                                   agcagGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCAttt                                                                                                            <  2:115868/151‑1 (MQ=255)                                      agGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGa                                                                                                         <  2:783281/151‑1 (MQ=255)                                            ggatcgAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGAcc                                                                                                    >  1:222865/1‑150 (MQ=255)                                             gatcgAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAg                                                                                                  <  2:601251/151‑1 (MQ=255)                                                 gAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCtt                                                                                               >  2:258002/1‑150 (MQ=255)                                                  aaCTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGt                                                                                             <  1:42082/151‑1 (MQ=255)                                                     tctGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTcc                                                                                           >  2:53560/1‑150 (MQ=255)                                                        ggTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCa                                                                                          <  1:577092/148‑1 (MQ=255)                                                        ggTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGt                                                                                       <  2:176050/151‑1 (MQ=255)                                                           ttCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCg                                                                                     <  2:921120/150‑1 (MQ=255)                                                               cttcttTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGt                                                                                   <  2:980638/148‑1 (MQ=255)                                                                   tttCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGt                                                                            >  2:401271/1‑151 (MQ=255)                                                                          tCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTacac                                                                     >  1:544632/1‑151 (MQ=255)                                                                           cAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTacac                                                                     <  2:437844/150‑1 (MQ=255)                                                                            aGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACAcgc                                                                   >  1:874500/1‑151 (MQ=255)                                                                             ggCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACAcgcg                                                                  <  1:483657/151‑1 (MQ=255)                                                                                 gACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCgct                                                                <  2:627941/149‑1 (MQ=255)                                                                                      aCATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCaa                                                         <  1:571635/151‑1 (MQ=255)                                                                                      aCATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCa                                                          <  1:974623/150‑1 (MQ=255)                                                                                           aCGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACAtt                                                     <  2:240618/150‑1 (MQ=255)                                                                                            cGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTg                                                    <  2:941179/150‑1 (MQ=255)                                                                                               aaGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAgg                                                <  2:175636/151‑1 (MQ=255)                                                                                                   cAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGTCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAAt                                            <  2:859212/151‑1 (MQ=255)                                                                                                   cAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCaa                                             <  2:583550/150‑1 (MQ=255)                                                                                                   cAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAAt                                            >  1:839465/1‑151 (MQ=255)                                                                                                    aaCGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAAt                                            >  2:793243/1‑150 (MQ=255)                                                                                                    aaCGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAAt                                            >  2:894398/1‑150 (MQ=255)                                                                                                     aCGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATAc                                          <  2:126347/151‑1 (MQ=255)                                                                                                       gAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGc                                        >  2:17089/1‑151 (MQ=255)                                                                                                         aaGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTc                                      >  1:814492/1‑151 (MQ=255)                                                                                                           gCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCca                                    <  1:859720/151‑1 (MQ=255)                                                                                                              tCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCcac                                   <  2:527870/149‑1 (MQ=255)                                                                                                              tCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCcac                                   >  1:926138/1‑149 (MQ=255)                                                                                                                  gttgttCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCTGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCt                             >  1:669729/1‑151 (MQ=255)                                                                                                                  gttgttCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCt                             <  2:296950/151‑1 (MQ=255)                                                                                                                     gttCAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGt                           <  2:380944/150‑1 (MQ=255)                                                                                                                        cAGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCACTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGc                        <  1:951350/150‑1 (MQ=255)                                                                                                                         aGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCAt                      >  1:592333/1‑151 (MQ=255)                                                                                                                         aGCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCAt                      <  2:44145/151‑1 (MQ=255)                                                                                                                          gCCAGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATg                     <  2:761135/151‑1 (MQ=255)                                                                                                                             aGAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTCCGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCAcgcg                   >  1:392851/1‑150 (MQ=255)                                                                                                                              gAATGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTc                 >  2:201247/1‑151 (MQ=255)                                                                                                                               aaTGGTTGCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCg                <  1:599156/151‑1 (MQ=255)                                                                                                                                      gCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACcagca         <  1:411943/151‑1 (MQ=255)                                                                                                                                      gCCGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACcagc          <  1:279064/150‑1 (MQ=255)                                                                                                                                       ccGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACcagt          >  2:581058/1‑148 (MQ=255)                                                                                                                                       ccGCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACcagcag        <  2:850245/151‑1 (MQ=255)                                                                                                                                         gCACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGAc      >  1:220362/1‑151 (MQ=255)                                                                                                                                          cACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACg     >  1:440738/1‑151 (MQ=255)                                                                                                                                          cACGACTAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACg     <  2:286788/151‑1 (MQ=255)                                                                                                                                             gactaATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACGGcc  <  1:780881/151‑1 (MQ=255)                                                                                                                                             gactaATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACGGcc  >  1:642786/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                    GCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCCTCTTCAGGATCGATTGTGCCGTTGGTTTCCATTTCGATGACCAGCTTGTCCAGGTCGGTACGCTGTTCTACACGCGCTGCTTCAACATTGTAGGCAATACGCTCCACAGGGCTGTAGCATGCGTCGACCAGCAGACGGCC  >  NC_000913/3440235‑3440525

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