Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,512,774	A→G	Y374C (TAC→TGC)	ydcS →	putative ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,512,774	0	A	G	100.0%	270.2 / NA	68	Y374C (TAC→TGC)	ydcS	putative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (48/20);  total (48/20)

AAGGGTGTAAAGCCAGTCCGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTACATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGTTCTTCTCTATG  >  NC_000913/1512633‑1512921                                                                                                                                              |                                                                                                                                                    aaGGGTGTAAAGCCAGTCCGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCAt                                                                                                                                            <  1:268283/151‑1 (MQ=255)   gggTGTAAAGCCAGTCCGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTa                                                                                                                                          <  1:526779/151‑1 (MQ=255)     gtgtAAAGCCAGTCCGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATg                                                                                                                                        >  2:460727/1‑151 (MQ=255)        tAAAGCCAGTCCGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGc                                                                                                                                     >  1:477476/1‑151 (MQ=255)          aaGCCAGTCCGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCgg                                                                                                                                   >  2:785310/1‑151 (MQ=255)           aGCCAGTCCGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGt                                                                                                                                  <  2:751411/151‑1 (MQ=255)                 tCCGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTa                                                                                                                             >  1:303592/1‑150 (MQ=255)                  ccGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCGGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAAc                                                                                                                           >  1:667217/1‑151 (MQ=255)                    gTTATTAGGCGAAAAAGGTTGTGAAACCAACGGATTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGTTAACtt                                                                                                                         >  2:82038/1‑151 (MQ=255)                     ttattaGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTc                                                                                                                        >  2:413491/1‑151 (MQ=255)                      tattaGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTc                                                                                                                        <  1:432577/150‑1 (MQ=255)                          aGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTgg                                                                                                                   >  1:517015/1‑151 (MQ=255)                             cGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGt                                                                                                                <  1:661978/151‑1 (MQ=255)                               aaaaaGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGc                                                                                                              >  2:600323/1‑151 (MQ=255)                                aaaaGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCt                                                                                                             >  2:1152334/1‑151 (MQ=255)                                 aaaGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCtt                                                                                                            >  2:804104/1‑151 (MQ=255)                                       tgtgAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGAc                                                                                                      >  1:479919/1‑151 (MQ=255)                                        gtgAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACg                                                                                                     >  1:1131965/1‑151 (MQ=255)                                          gAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGt                                                                                                    >  1:227278/1‑150 (MQ=255)                                          gAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGt                                                                                                    >  1:683123/1‑150 (MQ=255)                                          gAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGt                                                                                                    >  1:1123596/1‑150 (MQ=255)                                          gAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTa                                                                                                   <  1:460724/151‑1 (MQ=255)                                            aaCCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACg                                                                                                 >  2:438782/1‑151 (MQ=255)                                            aaCCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGGGCTTTATGACGTACg                                                                                                 >  1:74760/1‑151 (MQ=255)                                                 aCGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTg                                                                                            >  2:994760/1‑151 (MQ=255)                                                      tttAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGtt                                                                                       >  2:555776/1‑151 (MQ=255)                                                             aTTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAAc                                                                                >  2:691838/1‑151 (MQ=255)                                                              tttCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACg                                                                               >  2:304941/1‑151 (MQ=255)                                                                 cGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGtct                                                                            >  1:577768/1‑151 (MQ=255)                                                                 cGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGtc                                                                             <  2:869607/150‑1 (MQ=255)                                                                  gACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGtctc                                                                           <  1:447911/151‑1 (MQ=255)                                                                  gACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGtctc                                                                           >  2:872380/1‑151 (MQ=255)                                                                    cAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTcgc                                                                         >  1:879060/1‑151 (MQ=255)                                                                     aaaaTCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTcgcg                                                                        >  2:144296/1‑151 (MQ=255)                                                                      aaaTCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCgg                                                                       >  2:911047/1‑151 (MQ=255)                                                                      aaaTCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCgg                                                                       >  2:442316/1‑151 (MQ=255)                                                                       aaTCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGt                                                                      >  1:648784/1‑151 (MQ=255)                                                                        aTCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGtt                                                                     >  1:200889/1‑151 (MQ=255)                                                                         tCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTg                                                                    >  2:821688/1‑151 (MQ=255)                                                                              ttCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACgg                                                               >  2:784828/1‑151 (MQ=255)                                                                              ttCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACgg                                                               >  1:699627/1‑151 (MQ=255)                                                                              ttCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACgg                                                               >  2:493927/1‑151 (MQ=255)                                                                                   gAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACg                                                          >  2:609370/1‑151 (MQ=255)                                                                                    aaaaCGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACg                                                          >  1:13740/1‑150 (MQ=255)                                                                                     aaaCGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTg                                                        >  2:475909/1‑151 (MQ=255)                                                                                        cGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTgcgc                                                     <  2:355327/151‑1 (MQ=255)                                                                                                 cAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGt                                                >  2:1116321/1‑147 (MQ=255)                                                                                                 cAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGt                                                >  2:258054/1‑147 (MQ=255)                                                                                                  aGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATg                                           >  1:506076/1‑151 (MQ=255)                                                                                                  aGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATg                                           <  2:1031429/151‑1 (MQ=255)                                                                                                           cAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTGTGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGt                                   >  1:168757/1‑150 (MQ=255)                                                                                                           cAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTa                                  <  1:609361/151‑1 (MQ=255)                                                                                                                    tCCCTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATa                          <  2:1199429/150‑1 (MQ=255)                                                                                                                     cccTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATaa                         <  2:945509/150‑1 (MQ=255)                                                                                                                     cccTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATaa                         <  1:876853/150‑1 (MQ=255)                                                                                                                     cccTACAGTCGACGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGGGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGATGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATaa                         >  1:503532/1‑150 (MQ=255)                                                                                                                       cTACAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATaaa                        >  2:354548/1‑149 (MQ=255)                                                                                                                         aCAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAg                      <  2:886559/149‑1 (MQ=255)                                                                                                                          cAGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATg                   >  2:1014805/1‑151 (MQ=255)                                                                                                                           aGTCGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATgg                  >  1:1033080/1‑151 (MQ=255)                                                                                                                              cGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGa               <  1:144294/151‑1 (MQ=255)                                                                                                                              cGCTGGACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGa               <  1:459894/151‑1 (MQ=255)                                                                                                                                  ggACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGttc           <  1:3447/151‑1 (MQ=255)                                                                                                                                  ggACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGttc           >  2:989208/1‑151 (MQ=255)                                                                                                                                   gACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGttct          <  2:899160/151‑1 (MQ=255)                                                                                                                                   gACGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGttct          >  1:795932/1‑151 (MQ=255)                                                                                                                                     cGCAGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGttcttc        >  1:628736/1‑151 (MQ=255)                                                                                                                                        aGGATTGCATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGTTCTTctct     <  2:1057346/151‑1 (MQ=255)                                                                                                                                           attgcattgcCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGTTCTTCTCTATg  <  1:795934/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    AAGGGTGTAAAGCCAGTCCGTTATTAGGCGAAAAAGGTTGTGAAACCAACGGTTTTAACTATTTCGACAAAATCGCCTTCTGGAAAACGCCTATAGCAGAAGGGGGCAAGTTTGTTCCCTACAGTCGCTGGACGCAGGATTACATTGCCATTATGGGCGGTCGCTAACTTCGCTGGGGTGCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGTTCTTCTCTATG  >  NC_000913/1512633‑1512921

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