Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,316	A→G	V27A (GTT→GCT)	metJ ←	transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,316	0	A	G	100.0%	227.8 / NA	61	V27A (GTT→GCT)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/33);  total (28/33)

ATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTA  >  NC_000913/4128168‑4128463                                                                                                                                                     |                                                                                                                                                    aTCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGTACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGcc                                                                                                                                                   >  1:984267/1‑151 (MQ=255)   cGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGt                                                                                                                                                 >  1:585047/1‑151 (MQ=255)   cGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGt                                                                                                                                                 >  1:478348/1‑151 (MQ=255)     gCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTaa                                                                                                                                               <  2:745865/151‑1 (MQ=255)         catcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAGTGGAAGCCGTAATCtt                                                                                                                                           >  2:142890/1‑151 (MQ=255)         catcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCtt                                                                                                                                           >  1:848061/1‑151 (MQ=255)          atcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCttt                                                                                                                                          <  2:984288/151‑1 (MQ=255)            ccGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCttttt                                                                                                                                        <  1:839633/151‑1 (MQ=255)               gCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGAc                                                                                                                                     >  1:236968/1‑151 (MQ=255)                        gCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCAct                                                                                                                            <  1:757192/151‑1 (MQ=255)                            ggTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACtcttc                                                                                                                        <  2:848013/151‑1 (MQ=255)                            ggTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACtcttc                                                                                                                        >  1:383185/1‑151 (MQ=255)                              tAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTcttctt                                                                                                                      >  1:47485/1‑151 (MQ=255)                               aaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTg                                                                                                                     >  2:569727/1‑151 (MQ=255)                                  ggCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCg                                                                                                                  >  1:293786/1‑151 (MQ=255)                                   gCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGt                                                                                                                 >  1:154148/1‑151 (MQ=255)                                     aTGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGc                                                                                                               >  2:345815/1‑151 (MQ=255)                                         aGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAg                                                                                                           <  1:401636/151‑1 (MQ=255)                                           aaaCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCg                                                                                                         >  1:744592/1‑151 (MQ=255)                                            aaCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGt                                                                                                        <  2:778656/151‑1 (MQ=255)                                              cGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTAt                                                                                                      <  1:558831/151‑1 (MQ=255)                                               gtttCGTACAGCAGCTCGCTGTTAGAAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATAGGCGAGGATTTTTAACACCTTACGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCCTGCCGTGCTCAGCGt                                                                                                        >  2:868729/3‑148 (MQ=255)                                                cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATgg                                                                                                    >  2:3005/1‑151 (MQ=255)                                                  tCGCACATCAGCTCGCTGTTGGTAGCGTGACGCAGTTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATg                                                                                                     <  2:197727/148‑1 (MQ=255)                                                  tCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg                                                                                                   <  2:791840/150‑1 (MQ=255)                                                    gCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg                                                                                                <  1:977034/151‑1 (MQ=255)                                                    gCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg                                                                                                <  1:633182/151‑1 (MQ=255)                                                     cacaGCAGCTCGCTGTTGTTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCTTATGGGCTg                                                                                                <  2:804460/150‑1 (MQ=255)                                                          cagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGatatat                                                                                          <  2:740291/151‑1 (MQ=255)                                                          cagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGatatat                                                                                          >  1:4224/1‑151 (MQ=255)                                                           agcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATAtt                                                                                         <  2:910366/151‑1 (MQ=255)                                                           agcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATAtt                                                                                         >  1:846753/1‑151 (MQ=255)                                                                gCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTcgccg                                                                                    >  2:222327/1‑151 (MQ=255)                                                                    ttGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTcc                                                                                >  2:256296/1‑151 (MQ=255)                                                                      ggTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCa                                                                               <  2:883095/150‑1 (MQ=255)                                                                       gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAtt                                                                             <  1:556103/151‑1 (MQ=255)                                                                             tGACGCAGGTTGTTCTCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAAAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAg                                                                          <  2:1061199/148‑1 (MQ=255)                                                                                  cAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                                     >  1:447242/1‑148 (MQ=255)                                                                                   aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGAt                                                                 <  2:41080/151‑1 (MQ=255)                                                                                        gttCTCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTa                                                            <  2:894736/151‑1 (MQ=255)                                                                                           cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATc                                                         >  2:349702/1‑151 (MQ=255)                                                                                             ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                                       <  2:356754/151‑1 (MQ=255)                                                                                             ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                                       <  2:633884/151‑1 (MQ=255)                                                                                                  cgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCg                                                  <  2:900878/151‑1 (MQ=255)                                                                                                          gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATaa                                           >  1:62375/1‑150 (MQ=255)                                                                                                           cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAt                                         >  2:510637/1‑151 (MQ=255)                                                                                                            gtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                        >  1:638085/1‑151 (MQ=255)                                                                                                               cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgag                                     <  2:98830/151‑1 (MQ=255)                                                                                                                aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgaga                                    <  1:578721/151‑1 (MQ=255)                                                                                                                      gAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                              >  2:604705/1‑151 (MQ=255)                                                                                                                       aGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                             >  1:130103/1‑151 (MQ=255)                                                                                                                        ggATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAgaccagacca                            <  2:916028/151‑1 (MQ=255)                                                                                                                            ttttAACGCCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTTCCGTGCTCAGCGTATGTTCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGGCCAGACcacac                         <  1:543059/150‑1 (MQ=255)                                                                                                                            ttttAAAACCTTAAGCGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGTTATATTCGCCGCTCCATTCAGCCATGAGATACTTAATTTTCGTCGTCAATAAATTGAGACCAGACcac                           >  2:588333/1‑148 (MQ=255)                                                                                                                                 acacCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGa                   <  2:161095/151‑1 (MQ=255)                                                                                                                                     cTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgtg               >  2:250868/1‑151 (MQ=255)                                                                                                                                       tAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTgg              <  1:174983/150‑1 (MQ=255)                                                                                                                                        aaGAGGAATGTAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGCGACCAGACCACACAGTTGATGTGGGt            <  2:270891/151‑1 (MQ=255)                                                                                                                                           agGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt         <  1:934294/151‑1 (MQ=255)                                                                                                                                               aTGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAcc     <  2:305977/151‑1 (MQ=255)                                                                                                                                                  gAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTa  <  2:378458/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTA  >  NC_000913/4128168‑4128463

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