Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	568,660	C→T	intergenic (+13/‑242)	emrE → / → ybcK	DLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	568,660	0	C	T	100.0%	248.1 / NA	64	intergenic (+13/‑242)	emrE/ybcK	DLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/34);  total (30/34)

GTCCTGATTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCA  >  NC_000913/568519‑568804                                                                                                                                              |                                                                                                                                                 gTCCTGATTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGAc                                                                                                                                          <  2:545239/150‑1 (MQ=255)   ccTGATTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTaa                                                                                                                                       <  2:733876/151‑1 (MQ=255)   ccTGATTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTa                                                                                                                                        >  1:29599/1‑150 (MQ=255)   ccTGATTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTa                                                                                                                                        <  1:773854/150‑1 (MQ=255)         tAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAatat                                                                                                                                  <  2:98470/150‑1 (MQ=255)         tAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAatat                                                                                                                                  <  1:291634/150‑1 (MQ=255)          aGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATg                                                                                                                                 <  2:238225/150‑1 (MQ=255)          aGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATg                                                                                                                                 <  2:920830/150‑1 (MQ=255)            cTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATgg                                                                                                                                >  2:634131/1‑149 (MQ=255)                  gTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGctct                                                                                                                         <  2:316224/150‑1 (MQ=255)                   tCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCtt                                                                                                                        >  1:931121/1‑150 (MQ=255)                     aTGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCttat                                                                                                                      <  1:216487/150‑1 (MQ=255)                           aTTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCttatatttata                                                                                                               <  2:916209/151‑1 (MQ=255)                            tttttCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCttatatttatat                                                                                                              >  2:432069/1‑151 (MQ=255)                             ttttCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCttatatttatat                                                                                                              <  1:196624/150‑1 (MQ=255)                               ttCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATg                                                                                                             >  2:868135/1‑149 (MQ=255)                                 cGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAg                                                                                                           <  2:89217/149‑1 (MQ=255)                                   gCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTc                                                                                                       >  2:90689/1‑151 (MQ=255)                                      aaCGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGtt                                                                                                    <  1:513995/151‑1 (MQ=255)                                      aaCGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGtt                                                                                                    >  2:162895/1‑151 (MQ=255)                                       aCGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGttt                                                                                                   <  1:873819/151‑1 (MQ=255)                                         ggCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTa                                                                                                 <  1:600438/151‑1 (MQ=255)                                          gCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTAt                                                                                                <  2:456518/151‑1 (MQ=255)                                           cTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATg                                                                                               >  2:206708/1‑151 (MQ=255)                                                   gccagccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGt                                                                                       >  1:739553/1‑151 (MQ=255)                                                      agccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTaa                                                                                    <  1:206843/151‑1 (MQ=255)                                                      agccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTa                                                                                     >  2:55039/1‑150 (MQ=255)                                                        ccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATg                                                                                  >  2:12342/1‑151 (MQ=255)                                                         caTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGt                                                                                 >  2:946371/1‑151 (MQ=255)                                                               aGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAt                                                                           >  2:945526/1‑151 (MQ=255)                                                                    tgatgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtt                                                                          <  2:21380/147‑1 (MQ=255)                                                                      atgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                                    <  2:171159/151‑1 (MQ=255)                                                                      atgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                                    <  2:350917/151‑1 (MQ=255)                                                                           gATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTAtt                                                               <  2:592359/151‑1 (MQ=255)                                                                           gATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTAtt                                                               <  1:815273/151‑1 (MQ=255)                                                                            aTTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTc                                                              >  2:591817/1‑151 (MQ=255)                                                                            aTTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTc                                                              >  2:316385/1‑151 (MQ=255)                                                                             tttGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCt                                                             >  2:793448/1‑151 (MQ=255)                                                                                 tgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtg                                                         >  2:227051/1‑151 (MQ=255)                                                                                  gCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtg                                                         <  1:90903/150‑1 (MQ=255)                                                                                      gtgtgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCtt                                                    <  1:860753/151‑1 (MQ=255)                                                                                              attattAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGAtt                                            >  2:61485/1‑151 (MQ=255)                                                                                              attattAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGAtt                                            <  1:690314/151‑1 (MQ=255)                                                                                                 attaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGCCTAAAATATGGATGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTTTGCTTTCAAGAttc                                           <  1:411945/149‑2 (MQ=255)                                                                                                   taATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATgg                                        <  1:717025/150‑1 (MQ=255)                                                                                                   taATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGc                                       >  2:316188/1‑151 (MQ=255)                                                                                                     aTTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGt                                     <  2:646411/151‑1 (MQ=255)                                                                                                       ttATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTa                                    >  2:809933/1‑150 (MQ=255)                                                                                                        tATTGTCACGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTaa                                   <  2:365741/150‑1 (MQ=255)                                                                                                               aCGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTg                           >  2:306337/1‑151 (MQ=255)                                                                                                                cGAAGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGc                          <  1:132911/151‑1 (MQ=255)                                                                                                                   aGCACACCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAAt                       >  1:510427/1‑151 (MQ=255)                                                                                                                      acacCACATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTACTTTCAAGATTATGGCGTAAGAAAATTGCAATACg                    >  2:480940/1‑151 (MQ=255)                                                                                                                        accacATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGat                  <  2:648784/151‑1 (MQ=255)                                                                                                                         ccacATTAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGatt                 >  1:728292/1‑151 (MQ=255)                                                                                                                              ttAAAATAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgt            >  1:129971/1‑151 (MQ=255)                                                                                                                              ttAAAATAATTAGTTTTTAAAGGACTAAAATAAGGCGCCTCTTATATTTATACGAGCCCCGTTTTATGCATTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGCGCTTACAAGATTATGGAGTAAGCAAATTGCAATACGATTAttgt            >  1:487489/1‑151 (MQ=255)                                                                                                                                aaaaTAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgttg          >  2:615089/1‑151 (MQ=255)                                                                                                                                aaaaTAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgttg          >  2:178292/1‑151 (MQ=255)                                                                                                                                 aaaTAATTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATtgttgt         >  1:799700/1‑151 (MQ=255)                                                                                                                                   ataatTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGtat       <  2:12387/151‑1 (MQ=255)                                                                                                                                   ataatTTGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGtat       <  1:53739/151‑1 (MQ=255)                                                                                                                                         tGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCa  <  1:40384/150‑1 (MQ=255)                                                                                                                                         tGTTTTTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCa  >  1:368998/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                 GTCCTGATTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCA  >  NC_000913/568519‑568804

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