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%	245.9 / NA	61	K203T (AAA→ACA)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/34);  total (27/34)

GTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGCAT  >  NC_000913/4181707‑4181998                                                                                                                                                  |                                                                                                                                                   gTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACActgc                                                                                                                                               <  1:1035912/151‑1 (MQ=255) gTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACActgc                                                                                                                                               >  2:329890/1‑151 (MQ=255)     tcttTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACActgcct                                                                                                                                             <  1:9898/149‑1 (MQ=255)      cttTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCg                                                                                                                                          >  1:737292/1‑151 (MQ=255)        ttGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACActgcctgc                                                                                                                                           <  2:71624/148‑1 (MQ=255)              ccGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACcatcat                                                                                                                                  <  1:561843/151‑1 (MQ=255)                gACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTc                                                                                                                                <  2:387664/151‑1 (MQ=255)                   aaaGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACTGTCGCCGTACACTGCCTGCGACCATCATTCTgc                                                                                                                             >  1:39321/1‑151 (MQ=255)                    aaGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTgcg                                                                                                                            <  2:66059/151‑1 (MQ=255)                      ggTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTgcgcg                                                                                                                          <  2:954558/151‑1 (MQ=255)                         aaaaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGccc                                                                                                                       >  2:729440/1‑151 (MQ=255)                         aaaaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGccc                                                                                                                       <  1:531966/151‑1 (MQ=255)                           aaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTg                                                                                                                     <  2:581190/151‑1 (MQ=255)                              ccACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAAc                                                                                                                  >  1:855053/1‑151 (MQ=255)                               cACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAGCCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACt                                                                                                                 <  2:335837/151‑1 (MQ=255)                               cACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAAcg                                                                                                                 <  1:778821/151‑2 (MQ=255)                                       gggTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACcaca                                                                                                         <  1:53477/151‑1 (MQ=255)                                        ggTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACag                                                                                                        <  2:604356/151‑1 (MQ=255)                                         gTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTAcaccac                                                                                                          >  2:412271/1‑148 (MQ=255)                                          tAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGCGGCTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGCACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACag                                                                                                        >  1:701874/1‑149 (MQ=255)                                           aaaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGc                                                                                                     >  1:517538/1‑151 (MQ=255)                                            aaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCa                                                                                                    <  1:366878/151‑1 (MQ=255)                                             aGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGc                                                                                                     >  1:442246/1‑149 (MQ=255)                                              gTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGa                                                                                                  <  1:367990/151‑1 (MQ=255)                                                  tGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATcc                                                                                               >  1:822416/1‑150 (MQ=255)                                                   gTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATcc                                                                                               <  1:233590/149‑1 (MQ=255)                                                     ataACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCg                                                                                            >  1:1013361/1‑150 (MQ=255)                                                         cgcgcgTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTg                                                                                       <  2:856564/151‑1 (MQ=255)                                                         cgcgcgTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTg                                                                                       <  2:52778/151‑1 (MQ=255)                                                              gTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCtt                                                                                  >  1:630403/1‑151 (MQ=255)                                                                atcatcCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCttt                                                                                 <  2:676272/150‑1 (MQ=255)                                                                   atcCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGaa                                                                              >  2:378597/1‑150 (MQ=255)                                                                     cccTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGaaaaa                                                                           <  1:5746/151‑1 (MQ=255)                                                                         tACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGtt                                                                        <  2:798341/150‑1 (MQ=255)                                                                         tACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGt                                                                         <  1:93102/149‑1 (MQ=255)                                                                          aCCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTa                                                                       <  2:456828/150‑1 (MQ=255)                                                                            cGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACACCCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCt                                                                    >  1:707057/1‑151 (MQ=255)                                                                                 ttCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGa                                                                <  1:784320/150‑1 (MQ=255)                                                                                 ttCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGa                                                                <  2:629097/150‑1 (MQ=255)                                                                                  tCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGaa                                                               >  1:921618/1‑150 (MQ=255)                                                                                  tCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGaa                                                               <  1:412073/150‑1 (MQ=255)                                                                                      ggctggACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCg                                                          <  2:921974/151‑1 (MQ=255)                                                                                       gctggACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGt                                                         >  2:351881/1‑151 (MQ=255)                                                                                       gctggACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGt                                                         <  2:212715/151‑1 (MQ=255)                                                                                          ggACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGAt                                                      <  2:533488/151‑1 (MQ=255)                                                                                           gACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCTGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGAt                                                      >  2:296850/1‑150 (MQ=255)                                                                                                        aTCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGc                                            >  1:56803/1‑147 (MQ=255)                                                                                                          ccGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATgg                                      >  2:514834/1‑151 (MQ=255)                                                                                                          ccGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATgg                                      >  1:450841/1‑151 (MQ=255)                                                                                                           cGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATgg                                      >  1:794427/1‑150 (MQ=255)                                                                                                           cGAAAGACAACCTGTTCGTTCCTATCGACCGGCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGAt                                        <  2:899486/148‑1 (MQ=255)                                                                                                            gAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGaa                                    <  1:204568/151‑1 (MQ=255)                                                                                                               ggACAACCTGTTCGTACGTATCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTg                                 >  1:941036/1‑151 (MQ=255)                                                                                                                           cgtacgtaTCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGc                     >  2:386921/1‑151 (MQ=255)                                                                                                                              acgtaTCGACCGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTg                  >  2:1032572/1‑151 (MQ=255)                                                                                                                                        cGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGaa         >  1:723960/1‑150 (MQ=255)                                                                                                                                        cGTCGCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAAAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGaa         >  2:1034517/1‑150 (MQ=255)                                                                                                                                            gCCGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGc    <  1:774739/151‑1 (MQ=255)                                                                                                                                             ccGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGc    >  2:231545/1‑150 (MQ=255)                                                                                                                                              cGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGCAt  <  2:1013736/151‑1 (MQ=255)                                                                                                                                              cGTACACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGCAt  <  2:396515/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCCTGAACTACACCACAGAGCAGATCCTCGACCTGTTCTTTGAAAAAGTTATCTTTGAAATCCGTGATAACAAGCTGCAGATGGAACTGGTGCCGGAACGCCTGCGTGGTGAAACCGCAT  >  NC_000913/4181707‑4181998

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