Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,926,991	G→A	L5L (CTG→TTG)	asnC ←	DNA‑binding transcriptional dual regulator AsnC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,926,991	0	G	A	100.0%	139.7 / NA	37	L5L (CTG→TTG)	asnC	DNA‑binding transcriptional dual regulator AsnC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/26);  total (11/26)

CAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGA  >  NC_000913/3926746‑3927228                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              cAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATCCGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAaataat                                                                                                                                                                                                                                          <  1:538598/251‑1 (MQ=255)                 tCTTTGGCGCTCTTTAATATAATGACGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGCCTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTc                                                                                                                                                                                                                          >  2:336753/1‑251 (MQ=255)                      ggCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATg                                                                                                                                                                                                                     >  1:35892/1‑251 (MQ=255)                            ctTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGAGGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGa                                                                                                                                                                                                                 <  2:488756/249‑1 (MQ=255)                                 atatAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATtgatg                                                                                                                                                                                                           >  1:169599/1‑250 (MQ=255)                                     aaggCCGAAAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATCCGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGattc                                                                                                                                                                                                       <  2:431371/247‑1 (MQ=255)                                     aaTGCCGATAAAGCACCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGattc                                                                                                                                                                                                       <  1:442356/250‑1 (MQ=255)                                                agcagcCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAg                                                                                                                                                                                           <  1:463595/251‑1 (MQ=255)                                                      cTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCcttct                                                                                                                                                                                     >  1:370400/1‑251 (MQ=255)                                                        gcGTCATATCCGAGCTGCTTCGGGCTGCCATCAATACGCGCCCCGGTAATGTTCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCcttct                                                                                                                                                                                     <  2:161209/248‑1 (MQ=255)                                                            cATACCCGAGCTGCTTCGTGCTGACATCAATACGCGCCCCGGTAATGTTCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCATATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCtttttt                                                                                                                                                                                <  2:514584/250‑1 (MQ=255)                                                                       tgctTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGCCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCaaa                                                                                                                                                                     >  2:380391/1‑250 (MQ=255)                                                                        gctTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCaaaa                                                                                                                                                                    <  1:150293/250‑1 (MQ=255)                                                                        gctTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTGTAGCCTTCTTTTTTAATGAATCaaaa                                                                                                                                                                    <  1:513000/250‑1 (MQ=255)                                                                                        aaTACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTATCTACTCGAACTTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGCGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTC‑ATAGGATTTTCTTATGCGTATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTAt                                                                                                                                                  >  1:315699/1‑251 (MQ=255)                                                                                           aCGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTAtt                                                                                                                                                 <  1:494790/250‑1 (MQ=255)                                                                                             ccgcCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGGCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTTTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGaa                                                                                                                                              <  2:283684/250‑1 (MQ=255)                                                                                               gcCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAAGCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAAGGAATCAAAAGTGAGTTAGGCTTTTTATTGAATg                                                                                                                                            <  1:29858/251‑1 (MQ=255)                                                                                                              cccGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCAtgtgtg                                                                                                                             <  2:113883/251‑1 (MQ=255)                                                                                                                                          cGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGTATTTTCTTATGCGGATTTATGATTCATTCGATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATTAGTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCaa                                                                                                 <  1:571151/251‑1 (MQ=255)                                                                                                                                                       aggACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTCTTTTAGCCTTCTTTTTTAATGAATCAAATGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATACGCAACAGGACGCAggcg                                                                                    >  2:324638/2‑249 (MQ=255)                                                                                                                                                                 ccAAGTTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAAt                                                                          <  1:147758/251‑1 (MQ=255)                                                                                                                                                                     aTTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGCCGCAGGAGTATAAAAAATGaa                                                                       >  2:302788/1‑250 (MQ=255)                                                                                                                                                                     aTTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGaaa                                                                      <  1:446168/251‑1 (MQ=255)                                                                                                                                                                       tGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGaaaa                                                                     <  1:242242/250‑1 (MQ=255)                                                                                                                                                                        gTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAAcc                                                                   <  2:183359/251‑1 (MQ=255)                                                                                                                                                                                    cGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGc                                                        <  1:530103/250‑1 (MQ=255)                                                                                                                                                                                          aaGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTAGTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAaaca                                                  <  2:263351/250‑1 (MQ=255)                                                                                                                                                                                           aGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTTCCAaacaac                                                <  1:424488/251‑1 (MQ=255)                                                                                                                                                                                                    cgcATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTa                                       <  1:473773/251‑1 (MQ=255)                                                                                                                                                                                                             cATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGa                              <  2:35892/251‑1 (MQ=255)                                                                                                                                                                                                              aTTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGa                              >  1:305482/1‑250 (MQ=255)                                                                                                                                                                                                                 aaTGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAAtc                          <  1:156505/251‑1 (MQ=255)                                                                                                                                                                                                                           ggATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTTAGTTAGGCTTTTTATTGAATAATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtct                <  2:398942/251‑1 (MQ=255)                                                                                                                                                                                                                            gATGCCACGGTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtctc               <  1:557233/251‑1 (MQ=255)                                                                                                                                                                                                                                    ggTCCAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTg       >  1:28395/1‑251 (MQ=255)                                                                                                                                                                                                                                        cAGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGTTTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGga     >  1:132213/1‑249 (MQ=255)                                                                                                                                                                                                                                         aGATTGTCGATCAAATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGgaaga  <  2:232403/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              CAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGA  >  NC_000913/3926746‑3927228

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