Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,802	T→C	G297G (GGT→GGC)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,802	0	T	C	100.0%	178.0 / NA	54	G297G (GGT→GGC)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/31);  total (23/31)

TGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAG‑C‑TCAGTGATGGTGGTATG‑GGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTG  >  NC_000913/4626556‑4627045                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        tGTCATGGTTGGACACGTAACCAACGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTGTCTTAATTCGCGGAGATGAAGTGACCTCCGGCAG‑c                                                                                                                                                                                                                                                     <  1:483537/250‑1 (MQ=255)       ggttggacACGTAACCAAAGATGGCTCGCTGGCGGGCCCTCAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGGGATGAAGTGCCCTCCGGCAG‑C‑TCCGTg                                                                                                                                                                                                                                              >  2:473478/8‑250 (MQ=255)       ggttggacACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTAGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGCCTTCGCGATGACCGAACCGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCCTAAGTCGCGGAAGTGAAGTGACCTCCGGCAG‑C‑Tccgtt                                                                                                                                                                                                                                              >  2:404505/8‑246 (MQ=255)           ggacACGTAACCAAAGATGGCTCGGTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTGAATTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGTtggt                                                                                                                                                                                                                                         <  1:130828/248‑1 (MQ=255)                           aTGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACAATGTATCGCCTGTTCGGAGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCCTAAAAACCGCTTCGACGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGGTGGAGGTGTG‑GGAGGGAAcg                                                                                                                                                                                                                        >  1:425064/1‑251 (MQ=255)                             ggctcggtggcgggcCCTAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGATTCACGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGGCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAAcgcg                                                                                                                                                                                                                      <  2:302943/244‑1 (MQ=255)                                         ggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCtggt                                                                                                                                                                                                          >  1:91866/1‑251 (MQ=255)                                         ggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTGTTCGCACTTTACGCAGCCATAAAACCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCtggt                                                                                                                                                                                                          <  1:120539/251‑1 (MQ=255)                                           ccctaaaGTGCTGGAACCCTGTATCGACTGTTCGGTGCTTTTGGATGTCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑T‑TCCGTGATGGTGGTGTG‑TGAAGGAACGCGTccgtttctggtgg                                                                                                                                                                                                        <  1:171858/247‑12 (MQ=255)                                                aaGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGAtt                                                                                                                                                                                                   >  1:343875/1‑251 (MQ=255)                                                       ggAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGgc                                                                                                                                                                                              <  1:404747/249‑1 (MQ=255)                                                        gAACACTGTATCGACTGGTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCt                                                                                                                                                                                           <  1:161784/251‑1 (MQ=255)                                                                     aCTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACTGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCg                                                                                                                                                                              <  1:39303/251‑1 (MQ=255)                                                                           cGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGAtggtgg                                                                                                                                                                                                                                       <  1:29904/189‑1 (MQ=255)                                                                           cGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGAtggtgg                                                                                                                                                                                                                                       >  2:29904/1‑189 (MQ=255)                                                                              tGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAACAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCCG‑C‑TCCGTGCTGGTGGTGTG‑GGAAGGACCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCCCTCGGTGCTGGCg                                                                                                                                                                     >  2:212386/1‑251 (MQ=255)                                                                                      gATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGaatccgcg                                                                                                                                                             >  1:446654/1‑245 (MQ=255)                                                                                         gttgttGCCGACTCCCGTTTTGTCACGTTACGCGTCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATTACCGAACAGGGGCGGCGTGAAGTTAGCAACCCTTCGGCAATTTTCTTAAGGCGCGGAGATGAAGTTACCGCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑AGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATggttaatccgcgcc                                                                                                                                                           <  2:392566/245‑13 (MQ=255)                                                                                                ccGACTCCCGTTTTCGCACTTTACGCAGCCATAAAACCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGCTCACTCGATGATGGCGAATCCGCGCCGCGTGGcc                                                                                                                                                   >  2:280454/1‑250 (MQ=255)                                                                                                ccGACCCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑GCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCa                                                                                                                                                   <  2:76635/251‑1 (MQ=255)                                                                                                  gcctccctgtttCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGGGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGt                                                                                                                                                 <  2:406707/242‑1 (MQ=255)                                                                                                        cGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCt                                                                                                                                           <  1:9160/251‑1 (MQ=255)                                                                                                           tctccaccgttcctcagcaaaaaaaaCCGCTTCGGCGCGGTGAAGGAGCTGGGTGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCACCCGGTCGGCAAGTTTCGTAAGGCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGGCCGTTGCTGGTGGAGATTCAGGCGCTGGGAGATCACTTGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTgg                                                                                                                                         <  2:345242/229‑1 (MQ=255)                                                                                                                cactttacgcggcaatagaaaCCGCTTCGGCGCGGTGCATGAGCTGGTCGTCATCGCGATGACCGAACAGGGGCTGCGTGAAGTAAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGAGATGATGG‑AGGTGGAAAGCACGCGTCCGTTGCTGGAGGAGCTTCAGGCGCTGGACGATCACTCTATAATGGCGCATCAGCGCCGCATGGCAGTGGGGCTGGGACac                                                                                                                                    >  2:170939/19‑249 (MQ=255)                                                                                                                       cgcAGCCATAAAAACCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCt                                                                                                                            <  2:129521/251‑1 (MQ=255)                                                                                                                          aGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGGTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGGGGCAGTGGGTCTGGCACAAAACCGTATGGc                                                                                                                         >  2:764/1‑251 (MQ=255)                                                                                                                                   aaCCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCctgct                                                                                                                <  1:357774/251‑1 (MQ=255)                                                                                                                                    aCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCctgctg                                                                                                               >  1:59988/1‑251 (MQ=255)                                                                                                                                           cGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCtgtc                                                                                                         <  1:292824/250‑2 (MQ=255)                                                                                                                                                 ggTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCAcc                                                                                                  <  1:359403/251‑1 (MQ=255)                                                                                                                                                      aTGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGCTTCAGGCGCTGGTCGATCACTCGATGATGGCGTATCCGCGCCGCGTGGCAGTGGGGCTGGAACACAACCGTATGGCAATCCTGCTGGCTGCCTTGCACCGtacc                                                                                             >  2:413595/1‑248 (MQ=255)                                                                                                                                                             gggCGTCGGCGCGATGAACGAACAGGGGATGCGTGAAGTCAGCAACCCTTCGGCAATTTTCGTAAGTCGCGGAGATGAAGTGACCGCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACggtggt                                                                                       <  2:51192/250‑1 (MQ=255)                                                                                                                                                              ggCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCTCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTTCACCGTACCGGTGGTCt                                                                                     >  2:214503/1‑251 (MQ=255)                                                                                                                                                               gggccTTCGCGATGACCGAACAGGGGCTTCGGGAAGTCAGCAACCCTGCGGCAATTTTCGTAAGTCGCGGATATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTg                                                                                    <  2:423422/247‑1 (MQ=255)                                                                                                                                                                         gATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCg                                                                          >  1:72453/1‑251 (MQ=255)                                                                                                                                                                               cGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAgg                                                                    <  1:96014/251‑1 (MQ=255)                                                                                                                                                                                 aaCAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGAt                                                                  >  1:233169/1‑251 (MQ=255)                                                                                                                                                                                 aaCAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTAGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGGGGAGATTCAGGCACTGGTCGGTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGCCTGGCAATCCTGCTGGCTGTCTTTCACCGGCAGGGTGGTCTGCAAATGGCCGATCAggct                                                                  >  1:9517/1‑249 (MQ=255)                                                                                                                                                                                      gggCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGa                                                                         >  1:213920/1‑239 (MQ=255)                                                                                                                                                                                      gggCTGCGTGAAGTCAGCAACCCTGCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGa                                                                         <  2:213920/239‑1 (MQ=255)                                                                                                                                                                                           gCGTGAAGTCAGCATCCCTGCGGTAATTTTCTTAATGCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGGGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTtgt                                                          <  2:347776/249‑1 (MQ=255)                                                                                                                                                                                              ggAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGGTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACg                                                     <  2:413023/250‑1 (MQ=255)                                                                                                                                                                                                  gTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGAGGAGATTCAGGAGCTGGTCGATCACTCGATGATGGCGAATCCGCGACGAGTGGCAGTGGGGCTGGAAAAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTATGCAAATGGCCGTTAAGGATGTGTTGTTGAACG‑TGGt                                                >  1:334216/1‑251 (MQ=255)                                                                                                                                                                                                    cAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGaa                                                                                                                                       <  1:26416/163‑1 (MQ=255)                                                                                                                                                                                                     aGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTcgg                                             >  1:92991/1‑251 (MQ=255)                                                                                                                                                                                                               cGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGCTCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGa                                     >  1:49711/1‑249 (MQ=255)                                                                                                                                                                                                                 gCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGCTGATGGCGACTCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACCGTGGTCTGCCAATGGCCGATCCGGCTGTGTTTATGAACGCTGG‑CGGCGGCGTGACGGt                                 >  2:17055/1‑251 (MQ=255)                                                                                                                                                                                                                  cAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGt                                 >  1:341481/1‑250 (MQ=255)                                                                                                                                                                                                                             aGTCGGGGAGATGAAGTGACGTCCGGCAGGC‑GC‑GTGATGGTGGTGTG‑GGAAGGAACGCGTTCGTTGCTGGTGGAGAGTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAg                     <  2:137365/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggagATGAAGTGACCTCCGGCAG‑CCTC‑GTGATGGGGGGGGG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAAGCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGa               <  2:204339/248‑1 (MQ=255)                                                                                                                                                                                                                                      gATGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACtt            <  1:45882/251‑1 (MQ=255)                                                                                                                                                                                                                                       aTGAAGTGACCTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTa           <  1:389398/251‑1 (MQ=255)                                                                                                                                                                                                                                           aGTGACCTCCGGCAG‑C‑GCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGt       <  1:396544/251‑1 (MQ=255)                                                                                                                                                                                                                                                ccTCCGGCAG‑C‑TCCGTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTActg  <  1:414251/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        TGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAG‑C‑TCAGTGATGGTGGTATG‑GGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTG  >  NC_000913/4626556‑4627045

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: ‑