Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,559,761	T→C	S50S (TCA→TCG)	ddpC ←	D,D‑dipeptide ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,559,761	0	T	C	100.0%	219.5 / NA	66	S50S (TCA→TCG)	ddpC	D,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (39/27);  total (39/27)

GGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCGCTTAGCAT  >  NC_000913/1559614‑1559907                                                                                                                                                    |                                                                                                                                                   ggCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaa                                                                                                                                                 <  2:257715/151‑1 (MQ=255) ggCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaa                                                                                                                                                 <  2:871621/151‑1 (MQ=255)  gCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaa                                                                                                                                                 >  1:251214/1‑150 (MQ=255)   cGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaa                                                                                                                                                 >  2:561966/1‑149 (MQ=255)    gagaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaaaa                                                                                                                                              <  2:39597/151‑1 (MQ=255)    gagaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaaaa                                                                                                                                              <  2:211232/151‑1 (MQ=255)    gagaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaaa                                                                                                                                               >  2:835563/1‑150 (MQ=255)     agaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCCGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCTaaaaa                                                                                                                                              >  1:477629/1‑150 (MQ=255)     agaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaaaa                                                                                                                                              >  2:375862/1‑150 (MQ=255)     agaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATAGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaaaaa                                                                                                                                             <  1:37673/151‑1 (MQ=255)       aaTTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAatc                                                                                                                                           <  1:697419/151‑1 (MQ=255)        attgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAatca                                                                                                                                          <  2:514420/151‑1 (MQ=255)         ctgaatgCTGACTGCCGACAAGTACGCGTCTAAACAGATCGCGTCCCACTTCATAGGTGCCAAACCGGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGGTAAATAAATTGCGTTGGGATCATGCGGCGCCAGACACGGCGaaaaaacac                                                                                                                                          >  2:474841/2‑147 (MQ=255)              tgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAatcatcatc                                                                                                                                     >  2:48386/1‑150 (MQ=255)                 tGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAAt                                                                                                                                  >  2:850629/1‑150 (MQ=255)                     tGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCAt                                                                                                                             >  2:286151/1‑151 (MQ=255)                       ccGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTa                                                                                                                           <  1:865013/151‑1 (MQ=255)                         gACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAAt                                                                                                                         >  2:489825/1‑151 (MQ=255)                         gACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAAt                                                                                                                         >  2:256958/1‑151 (MQ=255)                             aGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACa                                                                                                                      >  1:761318/1‑150 (MQ=255)                                      cTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTg                                                                                                            >  1:810335/1‑151 (MQ=255)                                        aaaCAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTgcg                                                                                                          >  2:413117/1‑151 (MQ=255)                                           cAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTg                                                                                                       >  2:531044/1‑151 (MQ=255)                                            aGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTgg                                                                                                      >  2:823968/1‑151 (MQ=255)                                             gATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGt                                                                                                     <  1:225247/151‑1 (MQ=255)                                              aTCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTc                                                                                                    >  2:267316/1‑151 (MQ=255)                                                 gcgTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACg                                                                                                 >  1:805293/1‑151 (MQ=255)                                                   gTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGt                                                                                               >  1:891596/1‑151 (MQ=255)                                                    tCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCCATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTg                                                                                              >  1:208624/1‑151 (MQ=255)                                                     cccACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTg                                                                                              <  2:141957/150‑1 (MQ=255)                                                            cATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGa                                                                                        <  2:209570/149‑1 (MQ=255)                                                            cATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACt                                                                                      <  1:827161/151‑1 (MQ=255)                                                            cATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACt                                                                                      >  1:629259/1‑151 (MQ=255)                                                                 gTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCtt                                                                                 >  1:628632/1‑151 (MQ=255)                                                                     cAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTc                                                                               >  2:296651/1‑149 (MQ=255)                                                                        aCCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCAt                                                                          <  2:312853/151‑1 (MQ=255)                                                                          cAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATc                                                                         <  1:526321/150‑1 (MQ=255)                                                                             tGCGCCGCAGGAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCAt                                                                      >  1:407374/1‑150 (MQ=255)                                                                                gccgcAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATAcc                                                                   >  1:363075/1‑150 (MQ=255)                                                                                     aGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGtttt                                                             <  1:540854/151‑1 (MQ=255)                                                                                      gAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTg                                                            <  2:580956/151‑1 (MQ=255)                                                                                       aaGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGc                                                           >  1:584273/1‑151 (MQ=255)                                                                                          gcggcAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCAcc                                                        <  1:17358/151‑1 (MQ=255)                                                                                           cggcAAAAGGCGGGCGGTTAAATCACTGGCGCTGGGCTCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAGTAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTAAGCATCCATACCAGTTTTGCACCg                                                       >  2:154254/1‑151 (MQ=255)                                                                                                 aaGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAAt                                                 >  2:334257/1‑151 (MQ=255)                                                                                                  aGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATc                                                >  2:361000/1‑151 (MQ=255)                                                                                                          ggTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCgtttgt                                        <  1:679855/151‑1 (MQ=255)                                                                                                             tAAATCAATGACGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGtttt                                     >  2:407416/1‑151 (MQ=255)                                                                                                              aaaTCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGttttt                                    <  1:593159/151‑1 (MQ=255)                                                                                                              aaaTCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGttttt                                    >  2:750190/1‑151 (MQ=255)                                                                                                              aaaTCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGtttt                                     <  2:230764/150‑1 (MQ=255)                                                                                                              aaaTCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGtttt                                     <  2:428738/150‑1 (MQ=255)                                                                                                                  cAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTgg                                <  1:587072/151‑1 (MQ=255)                                                                                                                      ggCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGt                            <  1:748786/151‑1 (MQ=255)                                                                                                                                  aTGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGtt                >  2:35506/1‑151 (MQ=255)                                                                                                                                     cGTCGCAAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTcc             >  2:495296/1‑151 (MQ=255)                                                                                                                                      gTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTcc             >  1:378812/1‑150 (MQ=255)                                                                                                                                        cGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTc           <  2:831462/150‑1 (MQ=255)                                                                                                                                        cGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCg          >  1:535311/1‑151 (MQ=255)                                                                                                                                          ccagccaCGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCGCt        >  1:906178/1‑151 (MQ=255)                                                                                                                                          ccagccaCGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCGCt        <  1:419548/151‑1 (MQ=255)                                                                                                                                              ccaCGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCGCTTAGc    >  1:831718/1‑151 (MQ=255)                                                                                                                                              ccaCGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCGCTTAGc    >  2:894377/1‑151 (MQ=255)                                                                                                                                               caCGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCGCTTAGca   <  2:821722/151‑1 (MQ=255)                                                                                                                                                aCGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCGCTTAGcat  <  2:540016/151‑1 (MQ=255)                                                                                                                                                 cGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCGCTTAGcat  >  2:703844/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCGCTTAGCAT  >  NC_000913/1559614‑1559907

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