Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,606,883	G→A	E145K (GAG→AAG)	zntA →	zinc, cobalt and lead efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,606,883	0	G	A	100.0%	287.1 / NA	73	E145K (GAG→AAG)	zntA	zinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/36);  total (37/36)

TTCGTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGGAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTAT  >  NC_000913/3606737‑3607025                                                                                                                                                   |                                                                                                                                               ttCGTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGc                                                                                                                                             >  1:763005/1‑150 (MQ=255)  tCGTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGc                                                                                                                                             >  1:683075/1‑149 (MQ=255)   cGTGCACAAGTTGAATCTGCGCTGCAAAAAGCATGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGt                                                                                                                                          >  2:1001785/1‑151 (MQ=255)   cGTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGt                                                                                                                                          <  1:441391/151‑1 (MQ=255)    gTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGtt                                                                                                                                         >  1:939877/1‑151 (MQ=255)    gTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGt                                                                                                                                          >  2:82456/1‑150 (MQ=255)        acaAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAat                                                                                                                                     >  2:535221/1‑151 (MQ=255)          aaGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAatca                                                                                                                                   <  2:144830/151‑1 (MQ=255)              tGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATcc                                                                                                                                <  2:411005/150‑1 (MQ=255)              tGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATcc                                                                                                                                <  1:475154/150‑1 (MQ=255)              tGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCg                                                                                                                               >  1:702143/1‑151 (MQ=255)               gAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGt                                                                                                                              <  2:1089030/151‑1 (MQ=255)                actCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACCGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACACTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGtt                                                                                                                             <  2:668422/149‑1 (MQ=255)                aaTCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGtt                                                                                                                             <  2:983950/151‑1 (MQ=255)                aaTCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATAATCCGtt                                                                                                                             >  2:230506/1‑151 (MQ=255)                 aTCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTc                                                                                                                            >  2:42105/1‑151 (MQ=255)                    tGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCggg                                                                                                                         >  1:1045777/1‑151 (MQ=255)                     gcgcTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCgg                                                                                                                          <  1:519456/149‑1 (MQ=255)                        cTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAAc                                                                                                                     >  1:721026/1‑151 (MQ=255)                          gCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACt                                                                                                                    <  1:444249/150‑1 (MQ=255)                           cAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGACGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTg                                                                                                                   <  2:578257/150‑1 (MQ=255)                            aaaaaGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTgg                                                                                                                  <  2:1037870/150‑1 (MQ=255)                             aaaaGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCg                                                                                                                >  1:1234287/1‑151 (MQ=255)                             aaaaGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCg                                                                                                                >  1:967175/1‑151 (MQ=255)                             aaaaGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCg                                                                                                                >  1:1169458/1‑151 (MQ=255)                              aaaGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGt                                                                                                               >  1:1225889/1‑151 (MQ=255)                                  cAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTAt                                                                                                           <  1:777937/151‑1 (MQ=255)                                   aGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATc                                                                                                          >  2:1053006/1‑151 (MQ=255)                                   aGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATc                                                                                                          >  2:440679/1‑151 (MQ=255)                                   aGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATc                                                                                                          >  1:859753/1‑151 (MQ=255)                                       tATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGAAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATcgc                                                                                                        >  2:1168884/1‑149 (MQ=255)                                              tGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGc                                                                                                >  1:78520/1‑150 (MQ=255)                                                  cgATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGtt                                                                                           >  1:709132/1‑151 (MQ=255)                                                      gAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGc                                                                                       >  2:936624/1‑151 (MQ=255)                                                              cgccgAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCg                                                                               <  1:118196/151‑1 (MQ=255)                                                              cgccgAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCg                                                                               <  1:1228955/151‑1 (MQ=255)                                                                ccgAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGAt                                                                             <  1:380534/151‑1 (MQ=255)                                                                  gaagaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTg                                                                           >  1:1125602/1‑151 (MQ=255)                                                                    agaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCt                                                                         <  2:1045782/151‑1 (MQ=255)                                                                      aaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTc                                                                        <  2:518297/150‑1 (MQ=255)                                                                                aTCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCAtt                                                              <  1:1136277/150‑1 (MQ=255)                                                                                  cACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACg                                                           <  2:544085/151‑1 (MQ=255)                                                                                    cGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGt                                                         <  1:1090003/151‑1 (MQ=255)                                                                                     gCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGtt                                                        >  1:1217062/1‑151 (MQ=255)                                                                                       cTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTg                                                       >  1:711399/1‑150 (MQ=255)                                                                                         gAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATc                                                    >  1:321979/1‑151 (MQ=255)                                                                                          aaaGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATc                                                    >  1:1232557/1‑150 (MQ=255)                                                                                          aaaGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATc                                                    >  2:6316/1‑150 (MQ=255)                                                                                          aaaGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCa                                                   <  1:1185990/151‑1 (MQ=255)                                                                                              agaATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATc                                               <  1:1239697/151‑1 (MQ=255)                                                                                              agaATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATc                                               <  1:502991/151‑1 (MQ=255)                                                                                               gaATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATcc                                              <  1:893177/151‑1 (MQ=255)                                                                                               gaATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATcc                                              <  2:882310/151‑1 (MQ=255)                                                                                                    tGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGc                                           >  1:857886/1‑149 (MQ=255)                                                                                                      ccGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCt                                       <  1:914087/151‑1 (MQ=255)                                                                                                        gCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTAc                                     <  1:583509/151‑1 (MQ=255)                                                                                                         cTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACt                                    >  1:727320/1‑151 (MQ=255)                                                                                                          tGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACtt                                   >  2:1218054/1‑151 (MQ=255)                                                                                                          tGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACtt                                   >  1:1141800/1‑151 (MQ=255)                                                                                                           gATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTc                                  >  2:478859/1‑151 (MQ=255)                                                                                                            aTTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCg                                 <  1:606778/151‑1 (MQ=255)                                                                                                                 gCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCAtt                            <  1:1001786/151‑1 (MQ=255)                                                                                                                 gCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCAtt                            <  2:967173/151‑1 (MQ=255)                                                                                                                  cTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTg                           >  1:1004427/1‑151 (MQ=255)                                                                                                                     aTCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGaaa                        <  1:160400/151‑1 (MQ=255)                                                                                                                      tCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAAc                       >  1:1092974/1‑151 (MQ=255)                                                                                                                         tgatgaTGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCtt                    >  1:226763/1‑151 (MQ=255)                                                                                                                              atgGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATg                >  2:1075225/1‑150 (MQ=255)                                                                                                                                 gCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCg            <  2:730081/151‑1 (MQ=255)                                                                                                                                   aaTCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTa          <  1:163884/151‑1 (MQ=255)                                                                                                                                    aTCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAg         <  1:32555/151‑1 (MQ=255)                                                                                                                                           ggggTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCt    <  2:1098688/149‑1 (MQ=255)                                                                                                                                           ggggTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTAt  <  1:144184/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               TTCGTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGGAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTAT  >  NC_000913/3606737‑3607025

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