Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,181,852	A→C	K203T (AAA→ACA)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,181,852	0	A	C	100.0%	264.8 / NA	66	K203T (AAA→ACA)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/39);  total (27/39)

GGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGCATC  >  NC_000913/4181704‑4181999                                                                                                                                                     |                                                                                                                                                    ggCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTacac                                                                                                                                                   <  1:778141/151‑1 (MQ=255)   cGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACActg                                                                                                                                                 <  1:976451/151‑1 (MQ=255)    gTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACActgc                                                                                                                                                >  1:16563/1‑151 (MQ=255)        tcttTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACActgcctg                                                                                                                                             <  2:59864/150‑1 (MQ=255)        tcttTGACTCCAACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACActgcct                                                                                                                                              <  2:48477/149‑1 (MQ=255)         cttTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCg                                                                                                                                           >  2:433084/1‑151 (MQ=255)           ttGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCg                                                                                                                                           <  2:901884/149‑1 (MQ=255)           ttGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGa                                                                                                                                          <  1:764387/150‑1 (MQ=255)            tGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGAc                                                                                                                                         <  2:752667/150‑1 (MQ=255)               cTCCGACAAAGGTAAAACCCACTCTTCGGGGAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACcatc                                                                                                                                     >  1:698523/1‑151 (MQ=255)                tCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACcatc                                                                                                                                     >  2:176938/1‑150 (MQ=255)                 ccGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACcatcat                                                                                                                                   <  1:476004/151‑1 (MQ=255)                    aCAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTc                                                                                                                                 >  1:957338/1‑150 (MQ=255)                    aCAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCt                                                                                                                                <  2:518402/151‑1 (MQ=255)                      aaaGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTgc                                                                                                                              >  2:805296/1‑151 (MQ=255)                      aaaGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTgc                                                                                                                              >  1:911171/1‑151 (MQ=255)                        aGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTgcgc                                                                                                                            <  2:144695/151‑1 (MQ=255)                           tAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTGTCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTgcgcgc                                                                                                                          <  1:115194/150‑1 (MQ=255)                            aaaaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGccc                                                                                                                        <  2:998264/151‑1 (MQ=255)                             aaaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCt                                                                                                                       <  1:221037/151‑1 (MQ=255)                              aaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTg                                                                                                                      <  1:1024211/151‑1 (MQ=255)                                              aaaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACagag                                                                                                       <  1:306146/150‑1 (MQ=255)                                               aaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGc                                                                                                      >  1:1029330/1‑150 (MQ=255)                                                  tGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGc                                                                                                      >  1:409516/1‑147 (MQ=255)                                                   gCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGc                                                                                                      >  2:915023/1‑146 (MQ=255)                                                   gCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGAt                                                                                                  <  1:298495/150‑1 (MQ=255)                                                        ataACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGTCCTGAACTACACCACAGAGCAGATCCTc                                                                                              >  1:322652/1‑149 (MQ=255)                                                                   atcatcCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTg                                                                                 >  1:302454/1‑151 (MQ=255)                                                                    tcatcCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTg                                                                                 >  2:401006/1‑150 (MQ=255)                                                                     catcCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGaa                                                                               <  1:813672/151‑1 (MQ=255)                                                                     catcCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGaa                                                                               <  1:88195/151‑1 (MQ=255)                                                                     catcCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGa                                                                                <  2:348441/150‑1 (MQ=255)                                                                        cccTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGaaaa                                                                             <  2:11396/150‑1 (MQ=255)                                                                            tACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGtt                                                                         <  1:672783/150‑1 (MQ=255)                                                                            tACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGt                                                                          <  2:662698/149‑1 (MQ=255)                                                                             aCCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTAt                                                                       <  2:371799/151‑1 (MQ=255)                                                                             aCCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTAt                                                                       <  2:716024/151‑1 (MQ=255)                                                                              ccGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATc                                                                      <  2:854161/151‑1 (MQ=255)                                                                               cGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCt                                                                     <  2:201087/151‑1 (MQ=255)                                                                               cGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCt                                                                     <  1:836123/151‑1 (MQ=255)                                                                                  ggTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTg                                                                  <  2:783232/151‑1 (MQ=255)                                                                                    ttCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGa                                                                 >  2:655251/1‑150 (MQ=255)                                                                                         ggctggACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCg                                                           <  2:568669/151‑1 (MQ=255)                                                                                         ggctggACTGCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGGCCTGAGCTGCACCACAGAGCAGATCCTCGACCTGTGCTTTGAAAAAGTTATCTTTGAAATc                                                             <  1:225576/149‑1 (MQ=255)                                                                                              gACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGAt                                                       >  1:166212/1‑150 (MQ=255)                                                                                               aCTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCATCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATa                                                      >  2:533675/1‑150 (MQ=255)                                                                                                 ttCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAAc                                                    >  2:995076/1‑150 (MQ=255)                                                                                                  tCGAATTCGATCCGAATGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACa                                                   <  1:532958/150‑1 (MQ=255)                                                                                                   cGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACaa                                                  <  2:883426/150‑1 (MQ=255)                                                                                                    gAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGc                                                <  1:269604/151‑1 (MQ=255)                                                                                                     aaTTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCt                                               <  2:919644/151‑1 (MQ=255)                                                                                                     aaTTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCt                                               >  2:50430/1‑151 (MQ=255)                                                                                                              cGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCttt                                                                   >  1:254105/1‑122 (MQ=255)                                                                                                                 aGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGGGATAAGAAGAGGAAGAAGGaaat                                   >  1:148219/1‑149 (MQ=255)                                                                                                                  ggACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTg                                  <  1:603090/151‑1 (MQ=255)                                                                                                                      aaCCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGc                              >  1:882200/1‑151 (MQ=255)                                                                                                                        ccTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCg                            <  2:386953/151‑1 (MQ=255)                                                                                                                         cTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCgg                           >  2:361795/1‑151 (MQ=255)                                                                                                                           gTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGa                          >  2:818802/1‑150 (MQ=255)                                                                                                                              cgtacgtaTCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGc                      <  1:974967/151‑1 (MQ=255)                                                                                                                              cgtacgtaTCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGc                      <  2:655622/151‑1 (MQ=255)                                                                                                                                     aTCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCgt                <  1:805114/150‑1 (MQ=255)                                                                                                                                        gACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCgtggtg            >  2:927521/1‑151 (MQ=255)                                                                                                                                            gTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAAc        >  1:441831/1‑151 (MQ=255)                                                                                                                                               gCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGc     <  1:678038/151‑1 (MQ=255)                                                                                                                                                 cGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGCAt   >  2:477946/1‑151 (MQ=255)                                                                                                                                                  gTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGCATc  >  1:845192/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGCATC  >  NC_000913/4181704‑4181999

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