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%	237.2 / NA	61	N359S (AAC→AGC)	mdlB	putative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/37);  total (24/37)

CAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGT  >  NC_000913/471563‑471853                                                                                                                                                     |                                                                                                                                               cAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCa                                                                                                                                              <  1:291250/151‑1 (MQ=255)   ggCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCAtt                                                                                                                                            <  2:899710/151‑1 (MQ=255)   ggCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCAtt                                                                                                                                            <  1:452174/151‑1 (MQ=255)      tgttgtTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAAt                                                                                                                                         >  2:830484/1‑151 (MQ=255)       gttgttGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTaa                                                                                                                                          <  2:757920/149‑1 (MQ=255)       gttgttGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAAtc                                                                                                                                        <  2:102664/151‑1 (MQ=255)       gttgttGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAAt                                                                                                                                         <  1:489965/150‑1 (MQ=255)             gCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATctctct                                                                                                                                   <  1:298074/150‑1 (MQ=255)                ggTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGc                                                                                                                               <  1:854618/151‑1 (MQ=255)                   gAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCt                                                                                                                             <  1:746523/150‑1 (MQ=255)                     gcgcGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTcg                                                                                                                          >  1:350362/1‑151 (MQ=255)                          tgtTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCa                                                                                                                      >  1:59826/1‑150 (MQ=255)                          tgtTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCa                                                                                                                      >  2:657310/1‑150 (MQ=255)                           gtTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAAt                                                                                                                    <  2:497939/151‑1 (MQ=255)                           gtTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAAt                                                                                                                    <  1:561035/151‑1 (MQ=255)                            tttGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAAtt                                                                                                                   >  2:948944/1‑151 (MQ=255)                             ttGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAAttt                                                                                                                  >  1:965370/1‑151 (MQ=255)                                aaCTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTtgt                                                                                                               <  1:190086/151‑1 (MQ=255)                                 aCTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTtgtg                                                                                                              <  2:1032308/151‑1 (MQ=255)                                    gATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTgg                                                                                                             <  1:597180/149‑1 (MQ=255)                                    gATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTgg                                                                                                             <  2:510965/149‑1 (MQ=255)                                    gATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGgcg                                                                                                           <  1:46118/151‑1 (MQ=255)                                     aTGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGgcgc                                                                                                          <  2:257797/151‑1 (MQ=255)                                       ggacggacCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGGGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTg                                                                                                        <  1:485403/151‑1 (MQ=255)                                         acggacCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGt                                                                                                      >  2:428713/1‑151 (MQ=255)                                                        aaTATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAg                                                                                       <  2:828003/151‑1 (MQ=255)                                                         atatGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCGTACCGGCAGt                                                                                      <  2:200780/151‑1 (MQ=255)                                                                       atCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCAccc                                                                         >  1:431630/1‑150 (MQ=255)                                                                         cGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTc                                                                       <  1:393536/150‑1 (MQ=255)                                                                           cccGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGcc                                                                    <  1:1037680/151‑1 (MQ=255)                                                                            ccGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATTCCGGCAGTGGCAAAAGCACCCTCGCCa                                                                   <  2:863660/151‑1 (MQ=255)                                                                                                cGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTAc                                               >  1:249977/1‑151 (MQ=255)                                                                                                cGAAGTCGATAAAGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTAc                                               >  2:1024896/1‑151 (MQ=255)                                                                                                  aaGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTAccc                                             >  2:58121/1‑151 (MQ=255)                                                                                                  aaGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTAccc                                             <  1:947850/151‑1 (MQ=255)                                                                                                       gATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTa                                         >  2:48914/1‑150 (MQ=255)                                                                                                         tAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACg                                      >  2:362180/1‑151 (MQ=255)                                                                                                          aaCGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACgg                                     <  1:259477/151‑1 (MQ=255)                                                                                                             gtgtCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAg                                  <  1:568353/151‑1 (MQ=255)                                                                                                             gtgtCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAg                                  <  1:731517/151‑1 (MQ=255)                                                                                                             gtgtCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAg                                  <  1:850665/151‑1 (MQ=255)                                                                                                             gtgtCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAg                                  <  2:393353/151‑1 (MQ=255)                                                                                                              tgtCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAgg                                 >  1:67580/1‑151 (MQ=255)                                                                                                               gtCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAgg                                 >  2:236567/1‑150 (MQ=255)                                                                                                                tCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTg                               >  1:1042807/1‑151 (MQ=255)                                                                                                                 cATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTga                              <  1:271295/151‑1 (MQ=255)                                                                                                                         tATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGcc                      >  1:959365/1‑151 (MQ=255)                                                                                                                         tATCGCGATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGc                       <  2:709232/150‑1 (MQ=255)                                                                                                                            cgcgATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCtt                    <  1:1033656/150‑1 (MQ=255)                                                                                                                             gcgATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTg                   >  1:16962/1‑150 (MQ=255)                                                                                                                             gcgATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGa                  <  1:201844/151‑1 (MQ=255)                                                                                                                              cgATGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGAt                 >  2:432012/1‑151 (MQ=255)                                                                                                                                aTGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATgg               <  2:244492/151‑1 (MQ=255)                                                                                                                                aTGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATgg               >  2:142000/1‑151 (MQ=255)                                                                                                                                 tGACAATCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGgt              >  2:271547/1‑151 (MQ=255)                                                                                                                                      aTCTGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTcc         >  2:394436/1‑151 (MQ=255)                                                                                                                                         tGGTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCAtt      >  1:922648/1‑151 (MQ=255)                                                                                                                                          ggTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTcc         <  2:129162/147‑1 (MQ=255)                                                                                                                                          ggTGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCa        <  2:89135/148‑1 (MQ=255)                                                                                                                                            tGCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTa     >  2:1010410/1‑149 (MQ=255)                                                                                                                                             gCTAAAGAGCATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGt  <  2:124274/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               CAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACCGCGCCAGCAATATGGCAATGATGATCGCCCGTTACAGAGTGGCACCATCGAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGT  >  NC_000913/471563‑471853

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