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%	290.7 / NA	84	R297R (CGT→CGG)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (43/41);  total (43/41)

CTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGCGAA  >  NC_000913_3_bme_pgi/1927802‑1928094                                                                                                                                                    |                                                                                                                                                  cTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCgg                                                                                                                                                <  1:766026/151‑1 (MQ=255)  ttGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGt                                                                                                                                               >  1:138099/1‑151 (MQ=255)  ttGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGt                                                                                                                                               >  1:718892/1‑151 (MQ=255)    gCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGtat                                                                                                                                             >  2:778858/1‑151 (MQ=255)     cGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGtata                                                                                                                                            >  1:272151/1‑151 (MQ=255)      ggTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATaa                                                                                                                                           >  1:905214/1‑151 (MQ=255)       gtgaaTTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCAGCAGCCAGGCGGTATACAGCGTAAATCCTTCCAGCATGATGTCTTGTCTTGGCGCAATTAACTCGTCCCACTGTTGCTCATCACGCATACGGGCCCGGTATaa                                                                                                                                           >  4:40639/4‑150 (MQ=255)       gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATaa                                                                                                                                           >  4:280458/1‑150 (MQ=255)       gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAg                                                                                                                                          <  1:879734/151‑1 (MQ=255)           aTTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCaa                                                                                                                                      >  1:886013/1‑151 (MQ=255)             ttGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTTCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAg                                                                                                                                    <  2:711206/151‑1 (MQ=255)               ggCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGt                                                                                                                                   <  4:472219/150‑1 (MQ=255)                gcgcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGtt                                                                                                                                  >  4:438379/1‑150 (MQ=255)                 cgcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGtttt                                                                                                                                <  1:555864/151‑1 (MQ=255)                  gcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCATGCATACGGGTCCGGTATAAGCCAAAGtttt                                                                                                                                <  3:480959/150‑1 (MQ=255)                  gcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGttttt                                                                                                                               <  1:933446/151‑1 (MQ=255)                   cAAACTGGTTAACCCGCGCTCACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTCGCGGAATTAACTCTTCCCACCGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTg                                                                                                                              <  2:503458/151‑1 (MQ=255)                       cTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGt                                                                                                                          >  2:772870/1‑151 (MQ=255)                        tGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGt                                                                                                                          <  4:388744/150‑1 (MQ=255)                         ggTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTgg                                                                                                                        <  2:886013/151‑1 (MQ=255)                          gTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGc                                                                                                                       <  2:169113/151‑1 (MQ=255)                             aaCCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGta                                                                                                                                              <  1:773784/125‑1 (MQ=255)                             aaCCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGta                                                                                                                                              >  2:773784/1‑125 (MQ=255)                               cccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                                                                  >  1:776705/1‑151 (MQ=255)                               cccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                                                                  <  2:756302/151‑1 (MQ=255)                               cccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                                                                  >  1:659424/1‑151 (MQ=255)                               cccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                                                                  <  2:428728/151‑1 (MQ=255)                                 cgcgCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTgg                                                                                                                 >  3:471838/1‑150 (MQ=255)                                  gcgcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGaa                                                                                                               <  1:119020/151‑1 (MQ=255)                                      tGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCa                                                                                                           <  2:776705/151‑1 (MQ=255)                                        aCGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAg                                                                                                          <  3:453689/150‑1 (MQ=255)                                        aCGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGa                                                                                                         >  2:358432/1‑151 (MQ=255)                                             ctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaa                                                                                                    <  1:158917/151‑1 (MQ=255)                                              ttCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaa                                                                                                    <  3:186752/150‑1 (MQ=255)                                               tCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaaa                                                                                                   <  4:311468/150‑1 (MQ=255)                                               tCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAAc                                                                                                  >  2:224747/1‑151 (MQ=255)                                                    caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGaa                                                                                                               >  2:394254/1‑133 (MQ=255)                                                    caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGaa                                                                                                               <  1:394254/133‑1 (MQ=255)                                                         gccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTa                                                                                        >  2:358786/1‑151 (MQ=255)                                                          ccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTa                                                                                        >  4:55373/1‑150 (MQ=255)                                                                  tAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAg                                                                               <  2:256340/151‑1 (MQ=255)                                                                    aaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAgg                                                                              >  4:451415/1‑150 (MQ=255)                                                                    aaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAgg                                                                              <  3:581208/150‑1 (MQ=255)                                                                             aaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGt                                                                    <  2:138099/151‑1 (MQ=255)                                                                             aaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGt                                                                    <  2:718892/151‑1 (MQ=255)                                                                                ccTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGa                                                                  >  4:197236/1‑150 (MQ=255)                                                                                      aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTg                                                            >  3:578101/1‑150 (MQ=255)                                                                                      aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTg                                                            >  3:69085/1‑150 (MQ=255)                                                                                      aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTg                                                            <  3:476679/150‑1 (MQ=255)                                                                                      aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTCCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTg                                                            >  3:249369/1‑150 (MQ=255)                                                                                       gCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcg                                                          <  1:778858/151‑1 (MQ=255)                                                                                          tgatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcg                                                       >  2:182691/1‑151 (MQ=255)                                                                                               ttttCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCa                                                  >  1:765665/1‑151 (MQ=255)                                                                                                 ttCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAg                                                 <  3:545082/150‑1 (MQ=255)                                                                                                      gTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaac                                           <  2:272151/151‑1 (MQ=255)                                                                                                       tGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaac                                           <  3:212018/150‑1 (MQ=255)                                                                                                         gCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaacaaa                                        >  1:440923/1‑151 (MQ=255)                                                                                                          cGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaacaaac                                       <  1:772870/151‑1 (MQ=255)                                                                                                          cGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaacaaac                                       >  1:463915/1‑151 (MQ=255)                                                                                                          cGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaacaaa                                        >  4:61800/1‑150 (MQ=255)                                                                                                          cGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaacaaa                                        >  3:220133/1‑150 (MQ=255)                                                                                                                   ctctTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTc                                                                      <  3:465311/111‑1 (MQ=255)                                                                                                                   ctctTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTc                                                                      >  4:465311/1‑111 (MQ=255)                                                                                                                    tctTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGa                                                                                                <  3:354519/84‑1 (MQ=255)                                                                                                                    tctTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGa                                                                                                >  4:354519/1‑84 (MQ=255)                                                                                                                    tctTCCCACTGTTGCTCATCACGCATACGGGTCAGAGATATGTCAATGTGTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCAt                             >  2:285198/1‑151 (MQ=255)                                                                                                                         ccACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCCCGGCATCGGcc                        >  1:348510/1‑151 (MQ=255)                                                                                                                          cACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCa                       >  2:29183/1‑151 (MQ=255)                                                                                                                                tGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcg                 <  2:905214/151‑1 (MQ=255)                                                                                                                                  cTCATCACTCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGTCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgcg               <  1:732208/151‑1 (MQ=255)                                                                                                                                   tcatcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgcg               >  3:365807/1‑150 (MQ=255)                                                                                                                                    catcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTc             <  2:659424/151‑1 (MQ=255)                                                                                                                                     atcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgc                                                           >  2:83487/1‑104 (MQ=255)                                                                                                                                     atcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgc                                                           <  1:83487/104‑1 (MQ=255)                                                                                                                                     atcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTcc            >  2:634261/1‑151 (MQ=255)                                                                                                                                      tcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTcc            >  3:241297/1‑150 (MQ=255)                                                                                                                                         cGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCcagca        >  1:105005/1‑151 (MQ=255)                                                                                                                                           cATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                                                             <  2:666398/96‑1 (MQ=255)                                                                                                                                           cATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                                                             >  1:666398/1‑96 (MQ=255)                                                                                                                                            aTACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAgg      <  3:451415/150‑1 (MQ=255)                                                                                                                                            aTACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAgg      >  3:327247/1‑150 (MQ=255)                                                                                                                                               cGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGa                                                                                                         <  2:220037/48‑1 (MQ=255)                                                                                                                                               cGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGa                                                                                                         >  1:220037/1‑48 (MQ=255)                                                                                                                                                gggTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGCGaa  <  4:220133/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGCGAA  >  NC_000913_3_bme_pgi/1927802‑1928094

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: ‑