Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,394,081	Δ2 bp	coding (855‑856/1614 nt)	mppA →	murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,394,079	0	T	.	100.0%	291.6 / NA	60	coding (853/1614 nt)	mppA	murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (30/30);  total (30/30)
*	NC_000913	1,394,080	0	A	.	100.0%	291.6 / NA	60	coding (854/1614 nt)	mppA	murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (30/30);  total (30/30)

AAAACGGTACTGCAAAAAGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGA  >  NC_000913/1393932‑1394218                                                                                                                                                    ||                                                                                                                                           aaaaCGGTACTGCAAAAAGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTtacg                                                                                                                                          <  2:380046/151‑3 (MQ=255)  aaaCGGTACTGCAAAAAGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTtacgc                                                                                                                                         >  1:587467/1‑148 (MQ=255)         aCTGCAAAAAGTGTCCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCGCCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCGGGGGCAGGTT‑‑TACGCCGCCGCa                                                                                                                                <  1:41839/151‑1 (MQ=255)              aaaaaGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCg                                                                                                                           <  1:253632/151‑1 (MQ=255)              aaaaaGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCg                                                                                                                           >  1:450450/1‑151 (MQ=255)              aaaaaGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCg                                                                                                                           <  2:753388/151‑1 (MQ=255)               aaaaGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCgg                                                                                                                          <  1:218493/151‑1 (MQ=255)                 aaGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCggg                                                                                                                         >  2:186911/1‑150 (MQ=255)                 aaGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGa                                                                                                                        <  1:231745/151‑1 (MQ=255)                      aCCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCtat                                                                                                                   >  2:443682/1‑151 (MQ=255)                          tCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCtattat                                                                                                                <  2:967298/150‑1 (MQ=255)                                      aTCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAaca                                                                                                   <  1:695078/151‑1 (MQ=255)                                        cAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACg                                                                                                 <  1:290589/151‑1 (MQ=255)                                        cAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACg                                                                                                 <  2:528423/151‑1 (MQ=255)                                         aGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGc                                                                                                <  2:631083/151‑1 (MQ=255)                                           gAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCa                                                                                               >  2:464572/1‑150 (MQ=255)                                            aaTCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCaa                                                                                              >  2:825361/1‑150 (MQ=255)                                                   aGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGcc                                                                                      >  2:636146/1‑151 (MQ=255)                                                   aGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGcc                                                                                      <  1:870360/151‑1 (MQ=255)                                                    gCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCg                                                                                     <  2:315756/151‑1 (MQ=255)                                                     ccACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGa                                                                                    <  1:770446/151‑1 (MQ=255)                                                          aaGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCa                                                                               <  1:911833/151‑1 (MQ=255)                                                            gCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGa                                                                             <  1:644023/151‑1 (MQ=255)                                                              gTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATc                                                                           <  1:867857/151‑1 (MQ=255)                                                                tACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATc                                                                           >  2:670416/1‑149 (MQ=255)                                                                tACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCa                                                                          <  1:832557/150‑1 (MQ=255)                                                                 aCCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAgc                                                                        >  1:228247/1‑151 (MQ=255)                                                                    tCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGcgcg                                                                     >  2:962396/1‑151 (MQ=255)                                                                     cgcgGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAgcgc                                                                      <  1:487373/149‑1 (MQ=255)                                                                     cgcgGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAgcgc                                                                      <  1:369257/149‑1 (MQ=255)                                                                     cgcgGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGcgcg                                                                     >  2:937632/1‑150 (MQ=255)                                                                        gggggATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCg                                                                 <  2:98460/151‑1 (MQ=255)                                                                        gggggATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCg                                                                 >  1:890812/1‑151 (MQ=255)                                                                           ggATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCt                                                              <  2:391377/151‑1 (MQ=255)                                                                             atatTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTgg                                                            >  1:78780/1‑151 (MQ=255)                                                                                ttGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCAt                                                         >  2:898130/1‑151 (MQ=255)                                                                                 tGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCAtt                                                        >  1:854156/1‑151 (MQ=255)                                                                                       cACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTAt                                                  >  2:305633/1‑151 (MQ=255)                                                                                          cGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGAc                                               <  2:524827/151‑1 (MQ=255)                                                                                          cGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGAc                                               >  1:740601/1‑151 (MQ=255)                                                                                            aaTCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACg                                              <  2:587528/150‑1 (MQ=255)                                                                                            aaTCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACg                                              >  1:710084/1‑150 (MQ=255)                                                                                             aTCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGa                                             >  1:203125/1‑150 (MQ=255)                                                                                                 ttCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGAt                                            >  1:176117/1‑147 (MQ=255)                                                                                                   cccGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGAt                                            >  2:13768/1‑145 (MQ=255)                                                                                                   cccGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGAt                                            >  2:260517/1‑145 (MQ=255)                                                                                                   cccGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATc                                       <  1:660991/150‑1 (MQ=255)                                                                                                    ccGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATcgc                                     >  2:441004/1‑151 (MQ=255)                                                                                                             atGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCtgatg                            >  2:602802/1‑151 (MQ=255)                                                                                                              tGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCtgatg                            >  1:222649/1‑150 (MQ=255)                                                                                                                   cAGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGa                       >  1:313203/1‑150 (MQ=255)                                                                                                                    aGAAGCTGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGaaa                     <  2:450500/151‑1 (MQ=255)                                                                                                                        gatGTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTAGGGACCTATTATTATGCGTATAACACGCAAAAACGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGATGATAGATCGCCGCCTGATGACCGAAAAAg                  >  1:255694/3‑150 (MQ=255)                                                                                                                           gTTGAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTAtt              >  2:431734/1‑151 (MQ=255)                                                                                                                              gAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAgg           <  2:456026/151‑1 (MQ=255)                                                                                                                              gAAGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAg            <  2:417632/150‑1 (MQ=255)                                                                                                                                aGGATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAggg          >  2:109432/1‑150 (MQ=255)                                                                                                                                  gATATTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACg       >  1:433731/1‑151 (MQ=255)                                                                                                                                     atTCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGc    >  1:681063/1‑151 (MQ=255)                                                                                                                                       tCCGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCg   <  1:532141/150‑1 (MQ=255)                                                                                                                                        ccGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGa  <  1:883161/150‑1 (MQ=255)                                                                                                                                        ccGGGGCAGGTT‑‑TACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGa  <  2:348585/150‑1 (MQ=255)                                                                                                                                                    ||                                                                                                                                           AAAACGGTACTGCAAAAAGTGACCTTCCTGCCAATTAATCAGGAATCCGCAGCCACTAAGCGTTACCTCGCGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGA  >  NC_000913/1393932‑1394218

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