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	90.0%	112.2 / ‑3.2	41	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base A (19/17);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCTTTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCA  >  NC_000913/4558660‑4559155                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          gctACCGTTGCCGTCGTAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCattt                                                                                                                                                                                                                                                        <  1:277563/251‑1 (MQ=255)    cccGTTGCCGTCGGGGCTGAGGGGGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTAc                                                                                                                                                                                                                                                     <  2:80016/250‑1 (MQ=255)        gggCCGTCGGAGCGGAGGGTGACGCGTGCCAGCGGAATCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGG‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTc                                                                                                                                                                                                                                                 <  2:260371/249‑1 (MQ=255)         tGCCGTCGGGGCTGAGGGGGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTc                                                                                                                                                                                                                                                 <  2:205060/250‑1 (MQ=255)             gTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGccc                                                                                                                                                                                                                                             >  1:230826/1‑250 (MQ=255)                        gggTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑g                                                                                                                                                                                                                                <  1:108033/251‑1 (MQ=255)                                  tGCCAGCAGAAAACCGGCCTCAACGGCGCGGGCAATACCTTCGGAAGGGGCGACCGGGTCGTCAATGCCGCTGGTGATATCGATGGTGCCGCCTTTGAGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTAAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCA‑CACGGTCCTAAATAGGCTGTAA‑CGCCCTTTTACTGTCGCCCCTGTGGAAGAC‑TGTGACGcc                                                                                                                                                                                                                        >  2:196227/1‑249 (MQ=255)                                                    cTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTAGTCAATGCTGCTGGTGATATCGACGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGAACGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAAc                                                                                                                                                                                                    >  1:296930/1‑251 (MQ=255)                                                      gAACGGCTCGGGCAAGTCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTCCGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAAcc                                                                                                                                                                                                   <  2:78046/250‑1 (MQ=255)                                                                      aCCTTCGGCAGGGGCGACCGGTTCGTCAAGGCTGCTGGGGATACCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTg                                                                                                                                                                                  <  2:258311/251‑1 (MQ=255)                                                                        cTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCCAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTgcc                                                                                                                                                                                >  2:210584/1‑251 (MQ=255)                                                                                         ggTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTg                                                                                                                                                               >  2:297465/1‑251 (MQ=255)                                                                                          gTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAACCAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTg                                                                                                                                                               >  2:92718/1‑250 (MQ=255)                                                                                            gcGTCAATGCTGCTGGTGATATCGAGGGGGCCGCCTTGGCGCGCGAACTCCAGCGCCTGCTCAAACAGCGGTACGTTGCGGTGAACGGGGGT‑CGGCAGCAGCTGGCTGATCGGCATAGCGCAGTTT‑TCCAGCG‑GGTAATAAATAGGCTGAAA‑CGCCATTTTACTGTCGCCCATGTGGAACCC‑GGGGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATa                                                                                                                                                            <  2:203867/250‑1 (MQ=255)                                                                                                        ctgGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGc                                                                                                                                                >  1:75236/1‑251 (MQ=255)                                                                                                                  cGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACgg                                                                                                                                      >  1:113169/1‑251 (MQ=255)                                                                                                                   gATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGt                                                                                                                                     >  1:276784/1‑251 (MQ=255)                                                                                                                                ttGCGTGCGAACTACCGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGgc                                                                                                                        >  1:241959/1‑251 (MQ=255)                                                                                                                                ttGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATg                                                                                                                          <  1:229312/249‑1 (MQ=255)                                                                                                                                            cccAGCCCCTGCCCAAACAACGGTACTTTGCGGTTAACGTGGGT‑CGCCAGCGGCTTGCTGATCGGCACATCGCAGTTT‑ACCAGCC‑GCTCATAAATGGGCTGTAA‑CTCCTTTTTACTGCCGCCCTGGTGGAACAC‑GGTGACCCCGGGTTTCCCGCCGAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGcc                                                                                                            <  2:238369/250‑1 (MQ=255)                                                                                                                                              cAGCGCCTGCTCAAACAACGGGACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCAGTTCACGCCaa                                                                                                          <  1:92718/251‑1 (MQ=255)                                                                                                                                                  gcCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATcac                                                                                                      <  1:286775/251‑1 (MQ=255)                                                                                                                                                             caacGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCag                                                                                                                               >  2:228711/1‑215 (MQ=255)                                                                                                                                                             caacGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCag                                                                                                                               <  1:228711/215‑1 (MQ=255)                                                                                                                                                              aacGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAAt                                                                                          <  1:232470/251‑1 (MQ=255)                                                                                                                                                                  gTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGcc                                                                                      >  2:226490/1‑251 (MQ=255)                                                                                                                                                                     cGTTGCGGTTAACGTGGGT‑CGGCACCAGCTTGCTGAGCGGCACATCGCAGTTT‑CCCAGCA‑GGTCATAAATAGCCTGTAC‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATGCTGCCGCCATATTGGCCAGTTGAGAAACGTCCGGTGCGGCAGAAGGGTGATCAGAGAGGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACg                                                                                   <  1:143793/251‑1 (MQ=255)                                                                                                                                                                     cGTTGCGGTTAACG‑GGGTGCGGCAGCGGCTGGCTGAGCGGAACGTCGCAGTTT‑CCCGGCA‑GGTCAAAAAGAGGTTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGGGCGGCCAGGTTTACCGCCGAGCAACCCGCCAACGCGGGATTCGGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACg                                                                                   <  2:137843/251‑1 (MQ=255)                                                                                                                                                                         gCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCtt                                                                               >  1:12924/1‑251 (MQ=255)                                                                                                                                                                           ggTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACGGTCGCCCATGTGGAACCC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCtttt                                                                             <  1:271240/251‑1 (MQ=255)                                                                                                                                                                               aaCGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTcc                                                                          >  1:268980/1‑250 (MQ=255)                                                                                                                                                                                agcgggggccggccacagcTTGCTGACCGGCACAGCGC‑GTTTGTCCAGCA‑GGTAAAAAATGGGC‑GTAAACCCAATTTTCCGGCCGCCCAGGTGGAACAC‑GGTGCCGCCGGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCa                                                                         <  2:243222/236‑1 (MQ=255)                                                                                                                                                                                 gcgggggtcggcgactgctggttGATCGGCACTTCTCAGTTT‑TCCAGCA‑GTTAATAAAGAGGCTGTAA‑CGCCTTTTTACTTTCGCCCATGGGGAACAC‑GGGGACGCCGGGTTTGCCGCCGAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTGTTCCAc                                                                        <  2:230826/229‑1 (MQ=255)                                                                                                                                                                                    gggT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGa                                                                     >  2:253165/1‑250 (MQ=255)                                                                                                                                                                                     ggT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCC‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGAGCATTTCACGCCAATCACACGCTCGATAATCGCCACGTCTTTTTCCACGGAAc                                                                   >  1:270737/1‑251 (MQ=255)                                                                                                                                                                                            caggtgcttgcgtgtgggCACAGCGCAGTTT‑TCCAGCA‑GGTCATAAAAAGGCTGTAC‑CGCCATTTTACTGTCCCCCAGGGGGAACAC‑GGTGACGCCGGTTTTCCCGCCGAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGAGCGATAATCGCCACGTCTTTTTCCACGGAACCCGTaa                                                             <  2:8163/235‑1 (MQ=255)                                                                                                                                                                                             agcagcTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAAt                                                            >  1:170530/1‑251 (MQ=255)                                                                                                                                                                                                                caTCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGCGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATg                                          >  1:212713/1‑250 (MQ=255)                                                                                                                                                                                                                                        aTAAATAGGCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGCCATGATAAGCGCCGGTCAGCATCCAgg                   >  2:250006/1‑251 (MQ=255)                                                                                                                                                                                                                                             aagggtggaaaccccTTTTTACTGTCGCCCATG‑GGAACCCCGGTGCCGCCATGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTGTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCt               <  2:24149/238‑1 (MQ=255)                                                                                                                                                                                                                                               tgCTGTAA‑CGCCATTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGAt            >  2:161699/2‑251 (MQ=255)                                                                                                                                                                                                                                                          ccTTTTTACGGTCGCCCATGTGGAACGC‑GGTGACGCCAGTTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCa  <  2:272242/251‑1 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          GCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGT‑CGGCAGCAGCTTGCTGATCGGCACATCGCAGTTT‑TCCAGCA‑GGTCATAAATAGGCTGTAA‑CGCCTTTTTACTGTCGCCCATGTGGAACAC‑GGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCA  >  NC_000913/4558660‑4559155

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