Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	1,927,949	A→C	R297R (CGT→CGG)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	1,927,949	0	A	C	100.0%	306.0 / NA	87	R297R (CGT→CGG)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (41/46);  total (41/46)

TTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGC  >  NC_000913_3_bme_pgi/1927803‑1928091                                                                                                                                                   |                                                                                                                                               ttGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCgg                                                                                                                                             >  1:185796/1‑150 (MQ=255)  tGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGt                                                                                                                                            <  2:172121/150‑1 (MQ=255)    cGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGtat                                                                                                                                          >  3:519184/1‑150 (MQ=255)      gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATaa                                                                                                                                        >  2:575225/1‑150 (MQ=255)        tAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGc                                                                                                                                      <  2:749880/150‑1 (MQ=255)               gcgcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGtt                                                                                                                               <  4:257370/150‑1 (MQ=255)                cgcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGttt                                                                                                                              <  2:1097518/150‑1 (MQ=255)                 gcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGt                                                                                                                                            >  3:293881/1‑135 (MQ=255)                 gcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGt                                                                                                                                            <  4:293881/135‑1 (MQ=255)                 gcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGtttt                                                                                                                             >  2:213201/1‑150 (MQ=255)                     aCTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGcc                                                                                                                         <  3:345903/150‑1 (MQ=255)                             aCCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGt                                                                                                                 >  4:50731/1‑150 (MQ=255)                              cccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGtt                                                                                                                <  2:185796/150‑1 (MQ=255)                               ccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGTTGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                                                               <  1:227624/150‑1 (MQ=255)                               ccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                                                               >  3:550383/1‑150 (MQ=255)                                 gcgcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGa                                                                                                             >  3:678598/1‑150 (MQ=255)                                 gcgcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGa                                                                                                             <  3:266045/150‑1 (MQ=255)                                      gACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCa                                                                                                        <  4:519184/150‑1 (MQ=255)                                      cgaccactTCAACCACCAGCCAGTCGGTAAACAGCGAAAACCCTACCAGCATTATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGACCGGTATAAGCCAAATTTTTTGCCGTGGCGGTTGGAACGCa                                                                                                        <  4:294375/144‑1 (MQ=255)                                       aCGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAg                                                                                                       >  2:245584/1‑150 (MQ=255)                                          ctctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCg                                                                                                                                              <  1:503891/108‑1 (MQ=255)                                             ttCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGagatcg                                                                                                 >  3:294375/1‑145 (MQ=255)                                              tCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaaa                                                                                                >  2:137956/1‑150 (MQ=255)                                                 accaccAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGc                                                                                                                                      >  2:346750/1‑109 (MQ=255)                                                 accaccAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGc                                                                                                                                      <  1:346750/109‑1 (MQ=255)                                                 accaccAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGa                                                                                             >  2:1143867/1‑150 (MQ=255)                                                   caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATg                                                                                           <  2:901896/150‑1 (MQ=255)                                                   caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATg                                                                                           <  1:49798/150‑1 (MQ=255)                                                   caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGGCCGGTATAAGCCAAAGTTTGTGCCGTGGCGGTTGGAACGCAGATAAAAACGATg                                                                                           <  2:902362/150‑1 (MQ=255)                                                     ccagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCt                                                                                         >  2:46225/1‑150 (MQ=255)                                                      cagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTg                                                                                        <  1:709476/150‑1 (MQ=255)                                                       agccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTgg                                                                                       <  3:192246/150‑1 (MQ=255)                                                        gccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGt                                                                                      <  1:245584/150‑1 (MQ=255)                                                           agTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGt                                                                                   <  1:1047974/150‑1 (MQ=255)                                                              cGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGAt                                                                                <  4:75578/150‑1 (MQ=255)                                                              cGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGAt                                                                                >  2:417564/1‑150 (MQ=255)                                                               ggTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGttttt                                                                                                                            <  1:543386/105‑1 (MQ=255)                                                               ggTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGttttt                                                                                                                            >  2:543386/1‑105 (MQ=255)                                                                    aCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGc                                                                          >  4:605835/1‑150 (MQ=255)                                                                     cAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCt                                                                         >  2:532003/1‑150 (MQ=255)                                                                           aaaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTc                                                                   <  1:81639/150‑1 (MQ=255)                                                                            aaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAAc                                                                                               <  2:631583/121‑1 (MQ=255)                                                                            aaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAAc                                                                                               >  1:631583/1‑121 (MQ=255)                                                                             aCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGt                                                                 <  3:668842/150‑1 (MQ=255)                                                                             aCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGt                                                                 <  3:668642/150‑1 (MQ=255)                                                                                  tCCAGCATGATGTTTTCGCGTGGCGGTATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCt                                                            >  2:341774/1‑150 (MQ=255)                                                                                     aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTg                                                         >  2:994311/1‑150 (MQ=255)                                                                                      gCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgc                                                        <  1:341774/150‑1 (MQ=255)                                                                                      gCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgc                                                        <  1:575225/150‑1 (MQ=255)                                                                                        atgatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcgg                                                      <  3:694171/150‑1 (MQ=255)                                                                                            tgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcg                                                  <  4:164465/150‑1 (MQ=255)                                                                                              ttttCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGc                                                                                          >  1:518719/1‑108 (MQ=255)                                                                                              ttttCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGc                                                                                          <  2:518719/108‑1 (MQ=255)                                                                                              ttttCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcggc                                                >  2:358686/1‑150 (MQ=255)                                                                                               tttCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCa                                               <  3:50731/150‑1 (MQ=255)                                                                                                   gcgTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTgg                                                                                       >  1:39029/1‑106 (MQ=255)                                                                                                   gcgTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTgg                                                                                       <  2:39029/106‑1 (MQ=255)                                                                                                   gcgTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaaa                                           >  1:507455/1‑150 (MQ=255)                                                                                                    cgTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTgg                                                                                                                     >  4:621506/1‑75 (MQ=255)                                                                                                    cgTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTgg                                                                                                                     <  3:621506/75‑1 (MQ=255)                                                                                                     gTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaaaaa                                         <  1:137956/150‑1 (MQ=255)                                                                                                     gTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaaaaa                                         <  1:213201/150‑1 (MQ=255)                                                                                                             aTTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGc                                 <  1:455153/150‑1 (MQ=255)                                                                                                              ttAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCt                                >  1:174541/1‑150 (MQ=255)                                                                                                              ttAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCt                                >  1:355650/1‑150 (MQ=255)                                                                                                              ttAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCt                                >  2:446318/1‑150 (MQ=255)                                                                                                              ttAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTCGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGATGTATTCGTGTTCTTTGCGGCGCGGCAGAAAAACGAACGGCt                                >  2:466319/1‑150 (MQ=255)                                                                                                               tAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTc                               <  1:417564/150‑1 (MQ=255)                                                                                                                   tctTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGc                                                                                                         >  1:1135689/1‑72 (MQ=255)                                                                                                                   tctTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGc                                                                                                         <  2:1135689/72‑1 (MQ=255)                                                                                                                      tCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCg                        <  1:46225/150‑1 (MQ=255)                                                                                                                        ccACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTGGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGc                      >  1:641982/1‑150 (MQ=255)                                                                                                                        ccACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGc                      >  1:343390/1‑150 (MQ=255)                                                                                                                         cACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGcc                     <  2:355650/150‑1 (MQ=255)                                                                                                                         cACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGcc                     >  2:243196/1‑150 (MQ=255)                                                                                                                          aCTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCa                    >  3:573354/1‑150 (MQ=255)                                                                                                                              ttGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTc                >  1:1034109/1‑150 (MQ=255)                                                                                                                                 cTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaa                                                                                                 <  4:147313/66‑1 (MQ=255)                                                                                                                                 cTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaa                                                                                                 >  3:147313/1‑66 (MQ=255)                                                                                                                                 cTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgc             <  3:605835/150‑1 (MQ=255)                                                                                                                                 cTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgc             >  2:238719/1‑150 (MQ=255)                                                                                                                                  tcatcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGCCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgcg            <  1:532003/150‑1 (MQ=255)                                                                                                                                            tACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGc  >  1:591038/1‑150 (MQ=255)                                                                                                                                               gggTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCaa                                                                             >  3:363988/1‑72 (MQ=255)                                                                                                                                               gggTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCaa                                                                             <  4:363988/72‑1 (MQ=255)                                                                                                                                               gggTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTc          >  2:888233/1‑139 (MQ=255)                                                                                                                                               gggTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTc          <  1:888233/139‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               TTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGC  >  NC_000913_3_bme_pgi/1927803‑1928091

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