Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	601,732	G→A	E1007K (GAG→AAG)	cusA →	copper/silver export system RND permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	601,732	0	G	A	100.0%	142.5 / NA	41	E1007K (GAG→AAG)	cusA	copper/silver export system RND permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/23);  total (18/23)

CTGGCTTTATCGCCCTCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAGAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTC  >  NC_000913/601496‑601972                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 cTGGCTTTATCGCCCTCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAAGGGGGGCTGGTTCAAAGGTGATGAGCCgg                                                                                                                                                                                                                                    >  2:67337/1‑251 (MQ=255)       ttATCGCCCTCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTg                                                                                                                                                                                                                                <  1:28124/249‑1 (MQ=255)        tATCGCCCTCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGTGGCTGGTTCAAAGGTGATGAGCCGGATTGcc                                                                                                                                                                                                                              <  1:147548/250‑1 (MQ=255)            gCCCTCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGCGGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGGGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGTTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGcc                                                                                                                                                                                                                         <  2:352810/251‑1 (MQ=255)                tCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGCCGTTGAATAATCCACAACCATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCgatg                                                                                                                                                                                                                     >  1:32878/1‑251 (MQ=255)                           gccgccGAATTTGGCGTGGTGATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCAt                                                                                                                                                                                                          >  1:321807/1‑251 (MQ=255)                                  aaTTTGGCGTGGGGATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGAGGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGTCGTCATGATCAcc                                                                                                                                                                                                   <  2:59110/251‑1 (MQ=255)                                      tGGCGTGGTGATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCa                                                                                                                                                                                                <  1:38110/250‑1 (MQ=255)                                        gcgttggtaTGCTGATGTATTTACGTCACGGCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGGGCAGAAGCTGGATGAGGCGTTATATCCCGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCGGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGGTGATCGGCGGCATGATGACCGCAccg                                                                                                                                                                                             <  2:193985/243‑2 (MQ=255)                                          gtggtgATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCttt                                                                                                                                                                                           <  1:106691/251‑1 (MQ=255)                                           tggtgATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTg                                                                                                                                                                                          >  1:279001/1‑251 (MQ=255)                                            ggtgATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAGTCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTg                                                                                                                                                                                          <  1:352396/250‑1 (MQ=255)                                            ggggATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGGGCAGAAGCTGGAGGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGTGTGTCTGCTGCCGATTCTGTGGGGGACGGGGTGGGGTTCAAAGGTTATGAGCCGGATTGCCGCGCCGATGACTGGCGGCATTAGCTCCGAACCTTTg                                                                                                                                                                                          <  2:235695/247‑1 (MQ=255)                                               gctGCTAATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATAGCCGGTCTGCTGCCCATTCTGTGGGGAACGGTGGCTGGTTCAAAGGTGATGAGCCGGATTTCCGCGCTGGTTAATGGCGGCATTGGCACCGCCCCTTTGCTgg                                                                                                                                                                                      >  2:316830/3‑250 (MQ=255)                                                      aTGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGTCATGATCACCGCACCTTTGCTTTCGCTGtt                                                                                                                                                                               <  1:157256/251‑1 (MQ=255)                                                        gTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGtttt                                                                                                                                                                             >  2:308183/1‑250 (MQ=255)                                                              aCGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTgg                                                                                                                                                                                                                <  2:92574/210‑1 (MQ=255)                                                              aCGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTgg                                                                                                                                                                                                                >  1:92574/1‑210 (MQ=255)                                                                     gCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCggcggc                                                                                                                                                                <  1:95853/251‑1 (MQ=255)                                                                           gagGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATaa                                                                                                                                                          <  1:141759/251‑1 (MQ=255)                                                                                           tGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCa                                                                                                                                           <  1:306093/250‑1 (MQ=255)                                                                                              ataatCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCg                                                                                                                                        >  1:1470/1‑250 (MQ=255)                                                                                                                   agaaGAAGCTGGATGAGGCGTTAAATCACGGAGAGGTCCTGCGCGTGCGCCCGAAAGCGATGAGGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGGTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAg                                                                                                                  <  2:328752/248‑1 (MQ=255)                                                                                                                           cTGGATGAGGCGTTATATCACGGCGCGGTCCCGCGAGTGCGCCCGAAAGCGATGACGGTGGCGGCGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCCTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCCCCGACATCGGGTACGGACATCAAAGCAGGATAc                                                                                                          >  2:245551/1‑251 (MQ=255)                                                                                                                             ggATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTTTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATAccc                                                                                                        <  1:287479/251‑1 (MQ=255)                                                                                                                              gATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTtattat                                                                                                                                                                         >  1:318379/1‑185 (MQ=255)                                                                                                                              gATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTtattat                                                                                                                                                                         <  2:318379/185‑1 (MQ=255)                                                                                                                               aTGAGGAGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCCTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGGACGGAAATAAAAGCAGGATACCCCg                                                                                                      >  2:303160/1‑251 (MQ=255)                                                                                                                                                           gcgcGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGTGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGcga                                                                           <  1:277672/250‑1 (MQ=255)                                                                                                                                                            cgcgTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGcgacg                                                                         <  2:364119/251‑1 (MQ=255)                                                                                                                                                                      tgAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAATCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGctat                                                                >  2:10708/2‑247 (MQ=255)                                                                                                                                                                        aagtggATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTaaa                                                             <  2:249799/246‑1 (MQ=255)                                                                                                                                                                                         ggTGATTATAGCCGGTCAGCTGCCGATTATGTGGGGCAAGGGGGCTGGTTCAAAGGCGATGAGCCGGATTGCCGAGCCGATGATTGGCGGCATGATCCCCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCAATAAATGTTAAAAGGTGccccc                                            >  2:184561/1‑250 (MQ=255)                                                                                                                                                                                          gTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTc                                           <  2:289651/251‑1 (MQ=255)                                                                                                                                                                                                   gCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACaaaaa                                   >  2:15875/1‑250 (MQ=255)                                                                                                                                                                                                           gctgcCGATTCTGTGGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGAcacaca                          <  1:15875/251‑1 (MQ=255)                                                                                                                                                                                                                       gtgGGGAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGt              >  1:20064/1‑251 (MQ=255)                                                                                                                                                                                                                         ggggAACGGGGGCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGa            <  1:145901/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggggCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAAcgcg    >  2:331827/1‑251 (MQ=255)                                                                                                                                                                                                                                  gggCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGt   >  2:194474/1‑251 (MQ=255)                                                                                                                                                                                                                                    gCTGGTTCAAAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTc  >  2:323354/1‑250 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 CTGGCTTTATCGCCCTCGCCGGGGTCGCCGCCGAATTTGGCGTGGTGATGCTGATGTATTTACGTCACGCCATAGAGGCCGTGCCGTCGTTGAATAATCCACAAACATTCAGCGAGCAGAAGCTGGATGAGGCGTTATATCACGGCGCGGTCCTGCGCGTGCGCCCGAAAGCGATGACGGTGGCGGTGATTATCGCCGGTCTGCTGCCGATTCTGTGGGGAACGGGGGCTGGTTCAGAGGTGATGAGCCGGATTGCCGCGCCGATGATTGGCGGCATGATCACCGCACCTTTGCTGTCGCTGTTTATTATCCCGGCGGCGTATAAGCTGATGTGGCTGCACCGACATCGGGTACGGAAATAAAAGCAGGATACCCCGTTTAACCGTGTGGATTGTGTCTTGCGACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTC  >  NC_000913/601496‑601972

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