Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	598,133	A→G	T219A (ACC→GCC)	cusB →	copper/silver export system membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	598,133	0	A	G	100.0%	216.6 / NA	58	T219A (ACC→GCC)	cusB	copper/silver export system membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/40);  total (18/40)

CGGCGACCCAGACTGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCACCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTT  >  NC_000913/597999‑598269                                                                                                                                       |                                                                                                                                         cGGCGACCCAGACTGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGGTATTCGCCCCCCGATCGCCCCGCAAAAAAACCAGACTCGCCTTACGCCCAAAGCGCCCCTTGATGGCGTGGTCGCcccg                                                                                                                                     >  2:1329524/1‑140 (MQ=255) cGGCGACCCAGACTGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCcgcg                                                                                                                                     <  2:1458483/140‑1 (MQ=255)      aCCCAGACTGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATCCCGGCGGCGGGATTTCGTCCCCAGAACGCCACCCAAAAAATCCCGGCCCCCCTTTCGCCCAAAGCGCCCATTTACGGCGGGATCTCCGCCGTTGa                                                                                                                                >  2:1431279/1‑140 (MQ=255)         cAGACTGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAACGCCCATTGATGGCGTGATCGCCGCGTTTGATCt                                                                                                                             >  2:256874/1‑140 (MQ=255)            aCTGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTgcg                                                                                                                          <  2:859795/140‑1 (MQ=255)            aCTGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTgcg                                                                                                                          <  1:382464/140‑1 (MQ=255)              tGAAGTCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTgcgcg                                                                                                                        <  2:1172054/140‑1 (MQ=255)              tGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTgcgcg                                                                                                                        <  1:475847/140‑1 (MQ=255)              tGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTgcgcg                                                                                                                        <  2:920487/140‑1 (MQ=255)                      ttCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGgaatg                                                                                                                >  1:1031725/1‑140 (MQ=255)                      ttCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGgaatg                                                                                                                >  1:1092368/1‑140 (MQ=255)                          tGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAata                                                                                                            <  2:1301549/140‑1 (MQ=255)                          tGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAata                                                                                                            <  1:784656/140‑1 (MQ=255)                          tGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAACGCCCCTTGATGGCGTGGTCGCCCCGTTTGATCTGCGCGCGGGAAAGAata                                                                                                            >  2:956787/1‑140 (MQ=255)                              cGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGc                                                                                                        <  2:405615/140‑1 (MQ=255)                                gCTGCGACTGGCGGGAATGCCGGAGGCGGATATTTGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGGATGAATATCGccc                                                                                                      >  2:88933/1‑139 (MQ=255)                                gCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCa                                                                                                      <  1:973099/140‑1 (MQ=255)                                gCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCa                                                                                                      >  1:708048/1‑140 (MQ=255)                                 cTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCaa                                                                                                     <  1:1244979/140‑1 (MQ=255)                                          gCGGGAATGCCGGAGGCGGATATTCGCCCCCCGATCGCCACGCAAAAAAACCAGACTCGCTTTACGCTCCAAGCGCCCCTTGATGGCGTGATCGCCCCGGTTGATCTGCGCGCGGGAATGAAAATCGCCAAAGAAAACgg                                                                                            >  2:1224673/1‑139 (MQ=255)                                              gAATGCCGGAGGCGGATATTCGCCGCCCGATCGCCACGCAAAAAAACCCGACTCGCTTTACGCCCAAAGCGCCCCTTGATGGCGTGATCGCCCCGGTTGATCTGCGCCCGGGGAGGAAAATCGCCCAAGAAAACGTgggg                                                                                        >  2:791319/1‑138 (MQ=255)                                                        ggCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTc                                                                              <  1:515092/140‑1 (MQ=255)                                                        ggCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTc                                                                              <  1:414722/140‑1 (MQ=255)                                                         gCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCa                                                                             <  1:1331129/140‑1 (MQ=255)                                                         gCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCa                                                                             <  1:412685/140‑1 (MQ=255)                                                              tatTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTa                                                                        <  2:1057932/140‑1 (MQ=255)                                                              tatTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTa                                                                        <  1:85250/140‑1 (MQ=255)                                                                  cgccgcTTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGa                                                                    <  2:1103859/140‑1 (MQ=255)                                                                   gccgccTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGAc                                                                   <  1:588695/140‑1 (MQ=255)                                                                   gccgccTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGAc                                                                   <  1:1329524/140‑1 (MQ=255)                                                                   gccgccTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGAc                                                                   >  1:1322934/1‑140 (MQ=255)                                                                    ccgccTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGAcc                                                                  <  2:258046/140‑1 (MQ=255)                                                                    ccgccTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGAcc                                                                  <  2:351017/140‑1 (MQ=255)                                                                    ccgccTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGAcc                                                                  <  1:1431279/140‑1 (MQ=255)                                                                    ccgccTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGAcc                                                                  <  2:651180/140‑1 (MQ=255)                                                                    ccgccTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGAcc                                                                  <  1:1282808/140‑1 (MQ=255)                                                                       ccTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCgg                                                               >  2:61088/1‑140 (MQ=255)                                                                         tGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGgtg                                                             <  2:222026/140‑1 (MQ=255)                                                                         tGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGgtg                                                             <  1:1247991/140‑1 (MQ=255)                                                                           aTCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGtgtg                                                           >  1:1034138/1‑140 (MQ=255)                                                                                  cGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCAct                                                    >  1:598468/1‑140 (MQ=255)                                                                                   gCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCActg                                                   >  1:589558/1‑140 (MQ=255)                                                                                       aaaaTCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACtgctgc                                               <  2:120214/140‑1 (MQ=255)                                                                                             cAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATccc                                         >  2:1330934/1‑140 (MQ=255)                                                                                               gACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCgg                                       <  2:863724/140‑1 (MQ=255)                                                                                               gACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCgg                                       <  2:1031725/140‑1 (MQ=255)                                                                                               gACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCgg                                       <  1:479267/140‑1 (MQ=255)                                                                                                aCTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGa                                      <  1:713751/140‑1 (MQ=255)                                                                                                    gCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCt                                  <  2:1048484/140‑1 (MQ=255)                                                                                                         aCGCTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGc                             <  2:708048/140‑1 (MQ=255)                                                                                                            cTCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCctg                          >  1:953600/1‑140 (MQ=255)                                                                                                             tCAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCctgg                         <  1:1060580/140‑1 (MQ=255)                                                                                                              cAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCctggc                        >  2:290289/1‑140 (MQ=255)                                                                                                              cAAAGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCctggc                        >  1:1044688/1‑140 (MQ=255)                                                                                                                aaGCGCCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCctggctg                      >  2:105620/1‑140 (MQ=255)                                                                                                                  gcgcCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGt                    <  2:589558/140‑1 (MQ=255)                                                                                                                    gcCCATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGa                  >  2:252630/1‑140 (MQ=255)                                                                                                                     cccATTGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCCGGCTGGTGaa                 >  1:83653/1‑140 (MQ=255)                                                                                                                         ttGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGAt             >  1:606640/1‑140 (MQ=255)                                                                                                                         ttGATGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGAt             >  1:502544/1‑140 (MQ=255)                                                                                                                             tGGCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCt         <  2:1092368/140‑1 (MQ=255)                                                                                                                              ggCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTc        <  2:158744/140‑1 (MQ=255)                                                                                                                              ggCGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTc        <  2:1322934/140‑1 (MQ=255)                                                                                                                                cGTGATCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGc      <  1:722644/140‑1 (MQ=255)                                                                                                                                    aTCGCCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACGCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGGTGCCTCGCAGtt  >  1:221563/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                         CGGCGACCCAGACTGAAGGCATTCTTGAGCGGCTGCGACTGGCGGGAATGCCGGAGGCGGATATTCGCCGCCTGATCGCCACGCAAAAAATCCAGACTCGCTTTACGCTCAAAGCGCCCATTGATGGCGTGATCACCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTT  >  NC_000913/597999‑598269

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