Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,183,795	A→G	T851A (ACC→GCC)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,183,795	0	A	G	100.0%	248.4 / NA	64	T851A (ACC→GCC)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/26);  total (38/26)

AGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAA  >  NC_000913/4183650‑4183942                                                                                                                                                  |                                                                                                                                                    aGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCt                                                                                                                                                <  1:80155/151‑1 (MQ=255)   cGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTg                                                                                                                                               >  2:13206/1‑150 (MQ=255)    gTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGAc                                                                                                                                             >  2:678510/1‑151 (MQ=255)     ttCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGAc                                                                                                                                             >  1:457731/1‑150 (MQ=255)       cATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACAt                                                                                                                                           <  2:133267/150‑1 (MQ=255)           ccGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCg                                                                                                                                       >  1:653617/1‑150 (MQ=255)            cGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGa                                                                                                                                      >  1:1016367/1‑150 (MQ=255)               ggAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACg                                                                                                                                   >  2:569790/1‑150 (MQ=255)                gAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGTTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGt                                                                                                                                  <  1:909484/150‑1 (MQ=255)                 aaTGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTgg                                                                                                                                <  1:428043/151‑1 (MQ=255)                    ggTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTg                                                                                                                             >  1:636199/1‑151 (MQ=255)                      ttACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGa                                                                                                                            >  1:279378/1‑150 (MQ=255)                         cAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCt                                                                                                                        >  2:13281/1‑151 (MQ=255)                          aaCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTg                                                                                                                       >  2:718670/1‑151 (MQ=255)                          aaCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTg                                                                                                                       >  1:222952/1‑151 (MQ=255)                          aaCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTg                                                                                                                       >  1:30780/1‑151 (MQ=255)                            cTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTgcg                                                                                                                     >  1:1168514/1‑151 (MQ=255)                               cGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGct                                                                                                                   <  1:475434/150‑1 (MQ=255)                                gAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGctct                                                                                                                 >  1:891518/1‑151 (MQ=255)                                gAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGctct                                                                                                                 >  2:625704/1‑151 (MQ=255)                                gAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGctc                                                                                                                  >  2:541693/1‑150 (MQ=255)                                     cTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGAGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCaa                                                                                                             >  1:577528/1‑150 (MQ=255)                                       ccATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAAc                                                                                                           <  2:208719/150‑1 (MQ=255)                                         aTCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTgg                                                                                                        <  1:18441/151‑1 (MQ=255)                                               gTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGGGAAGCTGCGCTCTCCAAACTGGATGAAt                                                                                                  >  2:1154575/1‑151 (MQ=255)                                                 aTCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATcc                                                                                                >  1:58011/1‑151 (MQ=255)                                                 aTCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATcc                                                                                                <  2:436002/151‑1 (MQ=255)                                                   ccGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCgg                                                                                              >  1:693544/1‑151 (MQ=255)                                                     gAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTa                                                                                            <  2:154032/151‑1 (MQ=255)                                                      aGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTAt                                                                                           <  2:895210/151‑1 (MQ=255)                                                      aGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTAt                                                                                           >  1:148199/1‑151 (MQ=255)                                                        cGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCg                                                                                         >  2:307174/1‑151 (MQ=255)                                                          tgttgtTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGtt                                                                                       >  1:1068448/1‑151 (MQ=255)                                                               ttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACa                                                                                   >  1:884648/1‑150 (MQ=255)                                                                     aaGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGc                                                                            <  1:484423/151‑1 (MQ=255)                                                                       gACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCgg                                                                          >  2:591606/1‑151 (MQ=255)                                                                         ccGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGaa                                                                        <  1:325333/151‑1 (MQ=255)                                                                         ccGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGa                                                                         >  2:1094977/1‑150 (MQ=255)                                                                          cGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGa                                                                         >  2:487176/1‑149 (MQ=255)                                                                           gTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGt                                                                      <  2:1135336/151‑1 (MQ=255)                                                                           gTTTCACCACCATCCACATTCAGGAACTGGCGTGGGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAg                                                                       >  2:621451/1‑150 (MQ=255)                                                                               caccacCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGAcc                                                                  <  1:608929/151‑1 (MQ=255)                                                                                 ccaccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCgg                                                                <  2:1088381/151‑1 (MQ=255)                                                                                     caTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGc                                                            >  1:434037/1‑151 (MQ=255)                                                                                        ccACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGa                                                          <  1:126704/150‑1 (MQ=255)                                                                                              tCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCt                                                    <  2:638540/150‑1 (MQ=255)                                                                                               cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg                                                  >  2:1192915/1‑151 (MQ=255)                                                                                               cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg                                                  <  2:867310/151‑1 (MQ=255)                                                                                                 ggAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGGCCGGTGGCGACATTCtggtt                                                <  2:513901/151‑1 (MQ=255)                                                                                                  gAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg                                               <  1:551116/151‑1 (MQ=255)                                                                                                     cTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttggt                                             >  1:980320/1‑150 (MQ=255)                                                                                                        gCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaagg                                         >  2:556150/1‑151 (MQ=255)                                                                                                         cGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt                                        >  1:769890/1‑151 (MQ=255)                                                                                                            gtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggtaa                                      <  2:80331/150‑1 (MQ=255)                                                                                                                gtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGcc                                  <  1:502491/150‑1 (MQ=255)                                                                                                                 tCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCg                                 >  2:1171544/1‑150 (MQ=255)                                                                                                                   ccGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGaa                               >  1:1069607/1‑150 (MQ=255)                                                                                                                     gTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAgg                            >  1:465496/1‑151 (MQ=255)                                                                                                                                  gggTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCt                <  1:938026/150‑1 (MQ=255)                                                                                                                                   ggTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGa              <  1:282682/151‑1 (MQ=255)                                                                                                                                   ggTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGa              <  2:185577/151‑1 (MQ=255)                                                                                                                                          aaGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAcccc          >  1:1065298/1‑148 (MQ=255)                                                                                                                                             agaTCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCaga       >  1:1103860/1‑148 (MQ=255)                                                                                                                                               aTCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGaaa  <  1:197911/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAA  >  NC_000913/4183650‑4183942

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: ‑