Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,805,674	T→G	L287* (TTA→TGA)	proV →	glycine betaine/proline ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,805,674	0	T	G	100.0%	224.8 / NA	61	L287* (TTA→TGA)	proV	glycine betaine/proline ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (31/30);  total (31/30)

AAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTCCATCGATTCGCTTAAAACCGCGTTAACGCA  >  NC_000913/2805531‑2805818                                                                                                                                                |                                                                                                                                                 aaaTGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGt                                                                                                                                           >  1:184485/1‑151 (MQ=255)      gTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAAc                                                                                                                                      <  2:894605/151‑1 (MQ=255)       tGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAAcc                                                                                                                                     >  1:772507/1‑151 (MQ=255)        gAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAAccc                                                                                                                                    <  2:101304/151‑1 (MQ=255)        gAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAAccc                                                                                                                                    <  2:41333/151‑1 (MQ=255)        gAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAAcc                                                                                                                                     <  2:325227/150‑1 (MQ=255)         aaGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAGcccc                                                                                                                                   <  2:291180/151‑1 (MQ=255)          aGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCt                                                                                                                                  <  1:633513/151‑1 (MQ=255)           gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTg                                                                                                                                 >  2:388054/1‑151 (MQ=255)             ggTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGc                                                                                                                               >  2:299235/1‑151 (MQ=255)               tACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCtt                                                                                                                             >  1:624246/1‑151 (MQ=255)                  aGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCgg                                                                                                                          >  1:665518/1‑151 (MQ=255)                    gTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGcc                                                                                                                        <  2:639797/151‑1 (MQ=255)                      cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGccc                                                                                                                       >  2:531995/1‑150 (MQ=255)                       ggCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCAc                                                                                                                     <  2:693480/151‑1 (MQ=255)                          acacCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGtt                                                                                                                  <  2:890914/151‑1 (MQ=255)                           cacCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTc                                                                                                                 <  1:235111/151‑1 (MQ=255)                                   gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGa                                                                                                         <  1:545893/151‑1 (MQ=255)                                      aTTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAat                                                                                                      >  1:797292/1‑151 (MQ=255)                                        tctcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAatt                                                                                                     >  2:377423/1‑150 (MQ=255)                                         ctcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAattat                                                                                                   >  2:658430/1‑151 (MQ=255)                                             ataatCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGc                                                                                                >  2:364838/1‑150 (MQ=255)                                             ataatCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCa                                                                                               <  2:547386/151‑1 (MQ=255)                                                      cGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGa                                                                                      <  1:34204/151‑1 (MQ=255)                                                       gAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGAt                                                                                     <  1:725175/151‑1 (MQ=255)                                                        aaTGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATg                                                                                          >  1:858074/1‑145 (MQ=255)                                                           gATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCTCGTTCGGCACTGAAATTATTGCAGGATGAAGATcgcg                                                                                 >  2:582142/1‑151 (MQ=255)                                                           gATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCGTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATcgc                                                                                  <  1:751112/150‑1 (MQ=255)                                                            aTTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATcgcg                                                                                 >  1:391137/1‑150 (MQ=255)                                                            aTTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATcgcg                                                                                 <  2:242749/150‑1 (MQ=255)                                                             ttATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGa                                                                                >  2:923382/1‑150 (MQ=255)                                                                 gTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATg                                                                           <  2:830558/151‑1 (MQ=255)                                                                    cGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCt                                                                        <  1:438897/151‑1 (MQ=255)                                                                      tACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTa                                                                       >  1:580318/1‑150 (MQ=255)                                                                       aCCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACg                                                                     >  2:71365/1‑151 (MQ=255)                                                                           tcttcCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAACTTATTGCAGGATGAAGATCGCGAATATGGCTACGTTAt                                                                 >  1:40757/1‑151 (MQ=255)                                                                           tcttcCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGtt                                                                   <  2:292361/149‑1 (MQ=255)                                                                           tcttcCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTa                                                                  <  2:519220/150‑1 (MQ=255)                                                                                cGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAAc                                                            >  2:881298/1‑151 (MQ=255)                                                                                cGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAAc                                                            <  1:922722/151‑1 (MQ=255)                                                                                  tGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAAcgc                                                          >  2:511967/1‑151 (MQ=255)                                                                                    gCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCgg                                                        >  2:182903/1‑151 (MQ=255)                                                                                      gTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGta                                                      <  1:839574/151‑1 (MQ=255)                                                                                              tAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGtt                                               <  2:476119/150‑1 (MQ=255)                                                                                               aGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTg                                             <  2:44180/151‑1 (MQ=255)                                                                                                  cAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCg                                          <  1:826202/151‑1 (MQ=255)                                                                                                   aGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCg                                          <  1:653494/150‑1 (MQ=255)                                                                                                     gTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGgcg                                       >  1:747334/1‑151 (MQ=255)                                                                                                     gTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGgcg                                       >  2:740030/1‑151 (MQ=255)                                                                                                      tATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGgcgc                                      >  1:674950/1‑151 (MQ=255)                                                                                                      tATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGgcgc                                      >  2:321175/1‑151 (MQ=255)                                                                                                          cAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGtc                                  >  2:660640/1‑151 (MQ=255)                                                                                                             tGCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTcc                               >  1:580185/1‑151 (MQ=255)                                                                                                              gCGAAAGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTcc                               >  1:776925/1‑150 (MQ=255)                                                                                                                  aaGATATTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTCCATCGa                          <  1:21426/151‑1 (MQ=255)                                                                                                                     atatTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTCCATCGATTc                       <  2:217217/151‑1 (MQ=255)                                                                                                                       atTGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGGTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGAAATAAGTTTGTCGGCGCAGTCTCCATCgtt                         >  1:504183/1‑145 (MQ=255)                                                                                                                         tGCCCGCCGGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTCCATCGATTCGCt                    <  2:687618/150‑1 (MQ=255)                                                                                                                                cgGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTCCATCGATTCGCTTAAAAcc             >  1:780738/1‑150 (MQ=255)                                                                                                                                cgGACACCGAATGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTCCATCGATTCGCTTAAAACcg            >  2:13042/1‑151 (MQ=255)                                                                                                                                          aTGGCTGAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTCCATCGATTCGCTTAAAACCGCGTTAACgca  <  2:608146/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 AAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAAATTATTGCAGGATGAAGATCGCGAATATGGCTACGTTATCGAACGCGGTAATAAGTTTGTCGGCGCAGTCTCCATCGATTCGCTTAAAACCGCGTTAACGCA  >  NC_000913/2805531‑2805818

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