Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,886	C→T	N325N (AAC→AAT)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,886	0	C	T	97.9%	150.1 / NA	48	N325N (AAC→AAT)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (23/24);  minor base G (0/1);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAG‑C‑TCAGTGATGGTGGTATG‑GGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCA  >  NC_000913/4626651‑4627120                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            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)                                                                                                                                                                cgTGATGGTGGTGTG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGCCTTAGCGTTACTGCTGGCGATGGttt                                                                >  1:326597/2‑251 (MQ=255)                                                                                                                                                                            gtG‑GGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAAAGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGGTTCCAGCCTgcgtg                                                    >  2:381675/2‑249 (MQ=255)                                                                                                                                                                                        cgcgTCCGTTGCTGGTGGAGATCCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCCCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTCCTGATGGCGGTGGGTTCCAGCCTGCGGGACAGACCGCTg                                         >  2:267524/1‑251 (MQ=255)                                                                                                                                                                                           gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCa                                                             <  1:332495/222‑1 (MQ=255)                                                                                                                                                                                           gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGcc                                       >  1:63405/7‑250 (MQ=255)                                                                                                                                                                                           gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑GGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCa                                                             >  2:332495/7‑228 (MQ=255)                                                                                                                                                                                           gtccgttgctgGTGGAGATTCAGGCGCTGGGCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGcc                                       >  1:264871/7‑250 (MQ=255)                                                                                                                                                                                                 tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAGCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGCCGGCGTGAAGGTAACCGTAACCAGCGCCTGCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGa                                 <  1:445766/250‑1 (MQ=255)                                                                                                                                                                                                          agaTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCAACAGGATATGgtggtg                       >  1:182654/1‑251 (MQ=255)                                                                                                                                                                                                                        gTCGATCACTCGATGATGGCTAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCgg         <  1:425647/251‑1 (MQ=255)                                                                                                                                                                                                                        gTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCgg         >  1:48938/1‑251 (MQ=255)                                                                                                                                                                                                                               acgcGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCa  <  1:404505/248‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            CCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAG‑C‑TCAGTGATGGTGGTATG‑GGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCA  >  NC_000913/4626651‑4627120

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