Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	471,711	A→G	N359S (AAC→AGC)	mdlB →	putative multidrug ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	471,711	0	A	G	100.0%	258.8 / NA	66	N359S (AAC→AGC)	mdlB	putative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/33);  total (33/33)

GTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTC  >  NC_000913/471569‑471855                                                                                                                                               |                                                                                                                                                 gttgttGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAAt                                                                                                                                           <  1:630215/150‑1 (MQ=255) gttgttGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAAt                                                                                                                                           <  1:539453/150‑1 (MQ=255)    gttgCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAAtctct                                                                                                                                       >  2:451968/1‑151 (MQ=255)    gttgCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAAtctct                                                                                                                                       <  1:764173/151‑1 (MQ=255)    gttgCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAAtctc                                                                                                                                        <  1:219259/150‑1 (MQ=255)      tgCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAAtctctc                                                                                                                                      <  2:203990/150‑1 (MQ=255)         tGGAGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACACAGTGGCACCATAGAAGTCGATAACGTTTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCtg                                                                                                                                    >  2:96705/1‑149 (MQ=255)          ggTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGc                                                                                                                                 <  2:39850/151‑1 (MQ=255)             gAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCtt                                                                                                                              <  1:460295/151‑1 (MQ=255)              aGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTc                                                                                                                             >  1:649072/1‑151 (MQ=255)                cgcgTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTcgc                                                                                                                           >  1:65000/1‑151 (MQ=255)                   gtgtTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCa                                                                                                                        <  2:181535/151‑1 (MQ=255)                      tttGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAAtt                                                                                                                     >  1:285738/1‑151 (MQ=255)                       ttGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAAtt                                                                                                                     <  1:389241/150‑1 (MQ=255)                          aaCTGATGGACGGACCTCGCCAGCAATCTGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTtgt                                                                                                                 >  1:355459/1‑151 (MQ=255)                              gATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGgc                                                                                                              <  2:434566/150‑1 (MQ=255)                                tGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCt                                                                                                           >  1:742337/1‑151 (MQ=255)                                    cggacCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTc                                                                                                       <  1:529822/151‑1 (MQ=255)                                      gacCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCgg                                                                                                     >  1:282677/1‑151 (MQ=255)                                        ccGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGc                                                                                                   <  2:413386/151‑1 (MQ=255)                                           cgcCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATa                                                                                                <  1:359512/151‑1 (MQ=255)                                            gcCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATAc                                                                                               >  1:465993/1‑151 (MQ=255)                                                 cAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGc                                                                                           >  2:657008/1‑150 (MQ=255)                                                 cAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCa                                                                                          <  2:26255/151‑1 (MQ=255)                                                     atGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTgg                                                                                      <  1:740344/151‑1 (MQ=255)                                                     atGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTgg                                                                                      <  2:443541/151‑1 (MQ=255)                                                      tGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGc                                                                                     >  1:57640/1‑151 (MQ=255)                                                         cAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCaa                                                                                   >  2:408294/1‑150 (MQ=255)                                                          aaTGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCaaa                                                                                  >  2:456985/1‑150 (MQ=255)                                                           atgatgATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAg                                                                                <  2:613696/151‑1 (MQ=255)                                                                gatCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCAccc                                                                           >  2:30195/1‑151 (MQ=255)                                                                 atCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCt                                                                          <  2:169980/151‑1 (MQ=255)                                                                  tCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTc                                                                         >  1:760979/1‑151 (MQ=255)                                                                   cccccGTTACAGAGTGGCACCATCGAAGTCGATAACTTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCCCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCAcccgc                                                                         <  2:222740/148‑3 (MQ=255)                                                                   cGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTc                                                                         <  1:555383/150‑1 (MQ=255)                                                                          tACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCg                                                                        >  2:550457/1‑144 (MQ=255)                                                                          tACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGtt                                                                  <  1:162189/150‑1 (MQ=255)                                                                            cAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTAt                                                               <  2:603205/151‑1 (MQ=255)                                                                             agagTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCACTTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAACCGCACCCTCGCCAGTTTAtt                                                              >  2:142542/1‑151 (MQ=255)                                                                               agTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTg                                                             >  1:632249/1‑150 (MQ=255)                                                                                gTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGAt                                                           <  1:139479/151‑1 (MQ=255)                                                                                   gCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATggg                                                        >  2:350472/1‑151 (MQ=255)                                                                                          cGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTa                                                  <  1:462128/150‑1 (MQ=255)                                                                                          cGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTAc                                                 >  1:644130/1‑151 (MQ=255)                                                                                          cGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTAc                                                 <  2:433701/151‑1 (MQ=255)                                                                                           gAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTAcc                                                >  1:595268/1‑151 (MQ=255)                                                                                            aaGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAACCTCTCTGTGTCTTCGAGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTAccc                                               <  2:368501/151‑1 (MQ=255)                                                                                             aGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCg                                              >  1:115232/1‑151 (MQ=255)                                                                                                cGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCt                                            <  1:312219/150‑1 (MQ=255)                                                                                                     aCGTGTAATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCTCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACGCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACg                                        <  2:496576/149‑1 (MQ=255)                                                                                                      cGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGaa                                     <  2:457846/151‑1 (MQ=255)                                                                                                       gtgtCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGa                                      >  1:642699/1‑149 (MQ=255)                                                                                                         gtCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAgg                                   <  2:450349/150‑1 (MQ=255)                                                                                                           cATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTga                                <  1:641696/151‑1 (MQ=255)                                                                                                            aTTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTg                                 >  2:617501/1‑149 (MQ=255)                                                                                                              ttGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTg                                 >  2:323429/1‑147 (MQ=255)                                                                                                                    aTCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCt                       >  2:686662/1‑151 (MQ=255)                                                                                                                     tCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCtt                      >  2:763734/1‑151 (MQ=255)                                                                                                                        cgATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGAt                   <  2:198724/151‑1 (MQ=255)                                                                                                                            gACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGgtc               >  2:306127/1‑151 (MQ=255)                                                                                                                                  cTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTcc           >  2:597448/1‑149 (MQ=255)                                                                                                                                   tGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCAtt        >  2:557218/1‑151 (MQ=255)                                                                                                                                     gTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTaa      <  2:54348/151‑1 (MQ=255)                                                                                                                                      tGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTa       >  1:65528/1‑149 (MQ=255)                                                                                                                                       gCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCAGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGt    >  1:362897/1‑151 (MQ=255)                                                                                                                                          aaaGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTc  >  2:672396/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                 GTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTC  >  NC_000913/471569‑471855

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