Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,558,907	T→A	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,558,907	0	T	A	94.6%	94.3 / NA	37	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (16/19);  minor base C (1/1);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTT  >  NC_000913/4558666‑4559153                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       gTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                                                                                                                                                                                                               <  1:57854/251‑1 (MQ=255)    gCCGTCGGAGCTGAGGGTGACGCGTTCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACGCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGc                                                                                                                                                                                                                                            <  1:235555/251‑1 (MQ=255)    gCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGc                                                                                                                                                                                                                                            >  1:259250/1‑251 (MQ=255)         gggggcggAGGGTGCCGCGGGCCAGCGGAATCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                                                                                                                                                                                                       <  2:225246/244‑1 (MQ=255)          ggagctgggtgttccgctggCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                                                                                                                                                                                                      <  2:239003/232‑1 (MQ=255)                  gggTGACGCGTGCCAGCGGAATGCCCGCCTGAACGGCGCGGGCAATGCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACg                                                                                                                                                                                                                              <  1:5647/251‑1 (MQ=255)                                    gcaTACCCGCCTGAACGGCGCGGGCAATACCTTCTGCAGGTGCGACCGGCTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTTCTCAAACACCGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGATTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCTCCACTGTGGAACACGGTGACGCCAGGTTTAccc                                                                                                                                                                                                            >  2:84259/3‑250 (MQ=255)                                                gAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc                                                                                                                                                                                                 <  1:276857/250‑1 (MQ=255)                                                            cAATTCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCGCATCGCAGTTTGCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGAt                                                                                                                                                                                    <  1:132086/251‑1 (MQ=255)                                                                aCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTg                                                                                                                                                                                >  1:96485/1‑251 (MQ=255)                                                                    tcgtcAGGGGCGACCGTTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GTTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGGTTACCGCCGAGCAAACCGCCAACGCGGGATTCTgccg                                                                                                                                                                             <  1:204089/246‑1 (MQ=255)                                                                         aGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCatat                                                                                                                                                                       <  1:222107/251‑1 (MQ=255)                                                                            ggCGACCGGCTAGTCAAAGCTGCTGGTGATATCGATGGTGCAGCCTTTGCGCGCGAACTCCAGCGCCTGCTAAAACACCGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTATACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCCACGCGGGATTCTGCCGCCATATTgg                                                                                                                                                                    >  1:249943/1‑251 (MQ=255)                                                                                             ggctgctGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGTCCAGGTGATAAACGTcc                                                                                                                                                   <  1:219102/250‑1 (MQ=255)                                                                                                    ggggATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTCCGTTGCGGTTATCGTG‑GGTCGGCAGCATCGTGCTGATCGGCACATCCCATTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACCCCAGGTTTACCGCCGAGCAAACCTCCAACGCGGGAGTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCgg                                                                                                                                            <  1:139428/248‑1 (MQ=255)                                                                                                     gtgttttcaaTGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACTTG‑GGTCGGCAGCAGCTTGCTTATCGGCACATTGCAGTTTTCCAGCAGGTCATAAATAGGCGGTAACGCCATTTTTCTGTCGCCCATGTGGACCCCGGTGACGCCAGGTTTACCGCCGAGCAAACCCCCAACGCGGGATTCTGCCGCCATATTGGCCAGTTGATAAACGTCCGGTGCGGc                                                                                                                                           <  1:210346/242‑1 (MQ=255)                                                                                                          atCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAAc                                                                                                                                      <  1:94430/251‑1 (MQ=255)                                                                                                              aTGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTg                                                                                                                                  >  1:208563/1‑251 (MQ=255)                                                                                                                ggTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGAt                                                                                                                                <  1:250329/251‑1 (MQ=255)                                                                                                                         tttGCGCGCGAACTCCCGCGCCTGCTCCAACAACGGTCCGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGAAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCCTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCACACCGCCAACGCGGGCTTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATgg                                                                                                                       >  2:178045/1‑251 (MQ=255)                                                                                                                               cgcgAACCCCAGCGCCTGATCAAACAGCGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCGGATCGGCACAGCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGGCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGGCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCAtt                                                                                                                 <  2:294955/251‑1 (MQ=255)                                                                                                                                                              ccGTTGCGGTTAAC‑TGTGGTCGGCAACAGTTTTCTGACCGGCATATCTCATTTTTCCATCAGTTCATAAAGAGGCTGAAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACCCCAGTTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCAc                                                                                  <  1:162676/250‑1 (MQ=255)                                                                                                                                                                gTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATcagcagcacgttaccgcgtgtt                                                                                 >  1:232093/1‑230 (MQ=255)                                                                                                                                                                                 gTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAcc                                                                >  2:4924/1‑251 (MQ=255)                                                                                                                                                                                    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGt                                                             >  1:114222/1‑251 (MQ=255)                                                                                                                                                                                              ggCTGATCGGAAAACCGCAGTTTCCCAGCAGGTCAGAAAAAGTCTGTACCGCCATTTTCCTGTCGCCCAGGGGGAACCCGGTGACGCCAGTTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAAGCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCg                                                    <  2:252373/249‑1 (MQ=255)                                                                                                                                                                                                   aTCGGCACATCGCCGTTTTCAAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGCTTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAgg                                              >  2:299608/1‑251 (MQ=255)                                                                                                                                                                                                       gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACAGTGACGCCAGGCTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGGGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCCTTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGCACCCGTAATGGTGGTGGAAggggc                                           >  2:233967/1‑248 (MQ=255)                                                                                                                                                                                                            tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGCTTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCACTCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATgct                                      >  2:90950/1‑248 (MQ=255)                                                                                                                                                                                                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGATCATCTGGAACACGGTCACGCCAGGTTTACCCACGAGCAACCCGCCAACGCGGGATTCTGCCCCCATATTGGCCAGGTGATAAACg                                                                                                                                                      >  1:200076/1‑130 (MQ=255)                                                                                                                                                                                                                      ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGt                                                                                >  2:96678/1‑198 (MQ=255)                                                                                                                                                                                                                      ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGt                                                                                <  1:96678/198‑1 (MQ=255)                                                                                                                                                                                                                       cagcggttCATAAAAAGGCGTTACCGCCATTTTCCTGTCGCCCATGTGGAACCGGGTGACGCCAGGTTTCCCGCCGACCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCGTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCa                          <  2:259250/244‑1 (MQ=255)                                                                                                                                                                                                                                aaaaaTAGGCTGAACCGCACTTTTACTGTCGCCCATGTGGACCACGGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAgg                 <  1:233967/249‑1 (MQ=255)                                                                                                                                                                                                                                    aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCt             >  2:218432/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggCTGTAACGCAATTTTACTGTCGCGCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGGGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTgct          >  2:129438/1‑249 (MQ=255)                                                                                                                                                                                                                                               cGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGCGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttctt  >  1:177481/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       GTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTG‑GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTT  >  NC_000913/4558666‑4559153

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