Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	600,653	A→C	E647A (GAG→GCG)	cusA →	copper/silver efflux system, membrane component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	600,653	0	A	C	100.0%	231.7 / NA	63	E647A (GAG→GCG)	cusA	copper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (28/35);  total (28/35)

GCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACTGTGCTGGCGGATATCGACGC  >  NC_000913/600506‑600797                                                                                                                                                    |                                                                                                                                                 gCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGa                                                                                                                                               <  2:249950/151‑1 (MQ=255)      ccTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTgg                                                                                                                                          >  2:231157/1‑151 (MQ=255)         gAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGAt                                                                                                                                        >  2:613219/1‑150 (MQ=255)            gTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAac                                                                                                                                     >  2:636641/1‑150 (MQ=255)                 gcgGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTg                                                                                                                               >  1:846168/1‑151 (MQ=255)                 gcgGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTg                                                                                                                               >  2:647480/1‑151 (MQ=255)                      tATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCt                                                                                                                          >  1:639035/1‑151 (MQ=255)                       aTTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCt                                                                                                                          <  2:2522/150‑1 (MQ=255)                       aTTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTg                                                                                                                         >  2:372111/1‑151 (MQ=255)                        tttGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGc                                                                                                                        >  2:13812/1‑151 (MQ=255)                          tGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCg                                                                                                                      >  1:489199/1‑151 (MQ=255)                           ggCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCgg                                                                                                                     <  1:296797/151‑1 (MQ=255)                            gCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCggg                                                                                                                    >  1:907356/1‑151 (MQ=255)                               aaaCCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCt                                                                                                                 >  2:819916/1‑151 (MQ=255)                                  ccGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTgg                                                                                                               >  2:110489/1‑150 (MQ=255)                                   cGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCg                                                                                                             >  1:239331/1‑151 (MQ=255)                                    gggAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGa                                                                                                            <  1:385878/151‑1 (MQ=255)                                        aaGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCt                                                                                                        <  2:63167/151‑1 (MQ=255)                                           cGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCtgtg                                                                                                     >  1:340579/1‑151 (MQ=255)                                               aaCCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTg                                                                                                 >  2:104327/1‑151 (MQ=255)                                                aCCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTgc                                                                                                >  2:701522/1‑151 (MQ=255)                                                        cGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAAtt                                                                                        <  1:555927/151‑1 (MQ=255)                                                               gCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAAc                                                                                  <  1:631205/150‑1 (MQ=255)                                                                  ccGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGt                                                                               >  2:893515/1‑150 (MQ=255)                                                                      tGGAGATGGTAGAACCGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCGGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGCATCg                                                                           >  1:592869/1‑150 (MQ=255)                                                                       ggAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCg                                                                           <  2:917382/149‑1 (MQ=255)                                                                       ggAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCg                                                                           <  1:833949/149‑1 (MQ=255)                                                                       ggAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGa                                                                          <  1:207332/150‑1 (MQ=255)                                                                           aTGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGc                                                                     <  1:867446/151‑1 (MQ=255)                                                                           aTGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGc                                                                     >  1:318807/1‑151 (MQ=255)                                                                               tAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCtctc                                                                 >  2:864823/1‑151 (MQ=255)                                                                                aGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGCGGGTGCCGCCAATTCGTAACCGTATCGATATGctct                                                                  >  1:167257/1‑149 (MQ=255)                                                                                 gAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCaa                                                               <  1:104346/151‑1 (MQ=255)                                                                                 gAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCaa                                                               <  2:853934/151‑1 (MQ=255)                                                                                   aaCGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAAcc                                                             <  1:251106/151‑1 (MQ=255)                                                                                    aCGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCg                                                            <  1:58185/151‑1 (MQ=255)                                                                                    aCGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCg                                                            <  2:253097/151‑1 (MQ=255)                                                                                    aCGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCg                                                            >  1:854175/1‑151 (MQ=255)                                                                                      gACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGc                                                          <  2:887742/151‑1 (MQ=255)                                                                                       aCCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCa                                                         <  2:485864/151‑1 (MQ=255)                                                                                        ccATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCAt                                                        <  2:803058/151‑1 (MQ=255)                                                                                          aTCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCAtt                                                       >  1:317695/1‑150 (MQ=255)                                                                                            ccAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTaa                                                     <  2:814999/150‑1 (MQ=255)                                                                                            ccAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTa                                                      >  2:216784/1‑149 (MQ=255)                                                                                              cgCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGGCAAAATCATTGCGGAACTGGACAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCCATTCGTAACCGGATCGATATGCTCTCAACCGGCATTaaa                                                    >  1:56811/2‑149 (MQ=255)                                                                                                   aaGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCa                                             <  2:540760/151‑1 (MQ=255)                                                                                                    aGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCAt                                            <  2:861007/151‑1 (MQ=255)                                                                                                     gCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATc                                           <  1:258393/151‑1 (MQ=255)                                                                                                      ccGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCg                                          <  2:139615/151‑1 (MQ=255)                                                                                                      ccGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCg                                          <  2:309026/151‑1 (MQ=255)                                                                                                       cGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCgg                                         <  1:471444/151‑1 (MQ=255)                                                                                                          aGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCAt                                      <  2:456766/151‑1 (MQ=255)                                                                                                                 gTGGCGGCCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGtt                               <  2:516714/151‑1 (MQ=255)                                                                                                                      ggcCAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCgg                          >  2:702405/1‑151 (MQ=255)                                                                                                                        ccAGGCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCa                        <  1:403106/151‑1 (MQ=255)                                                                                                                            gCATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACtgt                    <  1:342108/151‑1 (MQ=255)                                                                                                                             cATGACGATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACtgtg                   <  1:150893/151‑1 (MQ=255)                                                                                                                                   gATGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACTGTGCTGGCg             <  1:25567/151‑1 (MQ=255)                                                                                                                                    aTGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACTGTGCTGGCgg            <  1:573186/151‑1 (MQ=255)                                                                                                                                    aTGGACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACTGTGCTGGCgg            >  1:761407/1‑151 (MQ=255)                                                                                                                                       gACAAAATCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACTGTGCTGGCGGat          >  1:497450/1‑150 (MQ=255)                                                                                                                                            aaTCATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACTGTGCTGGCGGATATCGAc    <  1:417118/151‑1 (MQ=255)                                                                                                                                               cATTGCGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACTGTGCTGGCGGATATCGAcgc  <  2:879991/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGGCACTGTGCTGGCGGATATCGACGC  >  NC_000913/600506‑600797

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