Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,943	G→C	V344V (GTG→GTC)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,943	0	G	C	98.5%	207.0 / NA	67	V344V (GTG→GTC)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (29/37);  minor base T (0/1);  total (29/38)
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

CGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACG  >  NC_000913/4626697‑4627178                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              cgcgGTGAATGAGCTGGTCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGc                                                                                                                                                                                                                                           <  1:519503/251‑1 (MQ=255) cgcgGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGCACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGc                                                                                                                                                                                                                                           >  1:480991/1‑251 (MQ=255)   cgGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTCATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCAc                                                                                                                                                                                                                                         <  1:348741/251‑1 (MQ=255)      ggAATGAGCTGGGCTTCGTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGt                                                                                                                                                                                                                                      <  2:558923/250‑1 (MQ=255)         aggAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGGGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACGCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCAc                                                                                                                                                                                                                                   <  2:310367/249‑1 (MQ=255)                gggCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTc                                                                                                                                                                                                                            <  1:355585/251‑1 (MQ=255)                 ggCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTGCGGCAATTTTCTGAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGGCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAAGCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCt                                                                                                                                                                                                                           <  2:338213/251‑1 (MQ=255)                    gTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGc                                                                                                                                                                                                                         >  1:18979/1‑250 (MQ=255)                      cTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCaa                                                                                                                                                                                                                       >  1:550322/1‑250 (MQ=255)                             aTGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCCCTCGATGATGGCGAATCCGCGCCGCGTGGCCGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGGCTTGCACCCTCACGGGGGTCTGCAACTGGCCg                                                                                                                                                                                                                >  2:674603/1‑250 (MQ=255)                              tGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGAt                                                                                                                                                                                                              >  1:505363/1‑251 (MQ=255)                                 ccGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAg                                                                                                                                                                                                           <  1:674603/251‑1 (MQ=255)                                 ccGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAg                                                                                                                                                                                                           <  2:18979/251‑1 (MQ=255)                                    aaCAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGGTGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAAAGCGTCCGTTGCTGGTGGGGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGAGCCGCGGGGCAGTGGGGCTGGAACAAAACCGTCAG‑GCAATCCTGCTGGCTGTCTTGCACCGTAACGGTGGTCTGCAAATGGCCGATCAGGAt                                                                                                                                                                                                        >  1:97828/1‑251 (MQ=255)                                           gCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGGCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTtg                                                                                                                                                                                                 <  2:269328/251‑1 (MQ=255)                                                       cAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGTCCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGATATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTc                                                                                                                                                                                     <  1:220549/250‑1 (MQ=255)                                                       cAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTcg                                                                                                                                                                                    >  1:81565/1‑251 (MQ=255)                                                             cccTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGgcggcg                                                                                                                                                                              <  1:130138/251‑1 (MQ=255)                                                                  cGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAg                                                                                                                                                                         <  1:458050/251‑1 (MQ=255)                                                                           tCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGCCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGGGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGa                                                                                                                                                                 >  1:34326/1‑250 (MQ=255)                                                                              tAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAAcc                                                                                                                                                             >  1:272108/1‑251 (MQ=255)                                                                                   cgcgGAGATGAAGTGACCTCCGGCCGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGCTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTTGCAGTGGGGCTGGAACAAATCCGTCTG‑GCAATACTGCTGGCTGTCTTGCACCGACACGGTGGTCTGCCAATGGCCGTTCAGGCTGTGTTTGTGAACG‑TGGTCGGCGGCGGGCAGGTAACCGaacacagagc                                                                                                                                                        >  2:123332/1‑243 (MQ=255)                                                                                      ggagATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCg                                                                                                                                                      >  1:438021/4‑250 (MQ=255)                                                                                      ggagATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGAGGAGATTCAGGCGCTGGACGATCACTCGATGATGGCGAATCCGCGCCGCGAGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGGAACGGTGGTCTGCACCTGGCCGATCAGTGTGTGTTTGTGAACG‑TGGTCGGCGGCGTGCAGGTAACCGCAAcccgcgccgg                                                                                                                                                     >  2:3872/4‑243 (MQ=255)                                                                                          aTGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACtt                                                                                                                                                  >  1:290367/1‑250 (MQ=255)                                                                                            gcaGTGACCTCCGCCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGTTGGAGATTCAGGCGTTGGTCGATCACTCGATGATGGCGACTCCGCGCCTCGTGGCAGTGGTGCTGGACCAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATTGCCGATCAGGATGTGTTTGTTAACG‑TGTTCGGCGGCTTGTAGGTAACCGAAACCAGCGCCGACTTAGc                                                                                                                                               <  2:348859/249‑1 (MQ=255)                                                                                             aaGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCg                                                                                                                                              >  1:79731/1‑251 (MQ=255)                                                                                                  cccTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑Act                                                                                                                                        <  2:550322/250‑1 (MQ=255)                                                                                                   ccTCCGGCAGCTCCGTGATGGTGGTGTGGGCAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGGCTGCAAATGGCCGTTCAGGATGTGTTTGTGACCG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCCGCGCCGACTTAGCGTT‑Actg                                                                                                                                       >  2:313528/1‑251 (MQ=255)                                                                                                      ccggcagctccggagtgggggtgtGGGAAGGATCGTGTCCGGTGCTGGGGGAGCATCAGGCCCTGGTCGAGCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTTGGGCTGGAACAAAACCGTCTG‑GCAATCCTGAGGGCGGTCTTGCACCGGCACGGTGGTCTGCAAATGGCCGCTCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAATCCTGCGCCGACTTAGCGTT‑Tctgct                                                                                                                                     <  2:527580/227‑1 (MQ=255)                                                                                                       cggcaGCTCCGTGATGGTGGTGTGGGAAGTAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAATCCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAATCCAGCGCCGACTTATCGGT‑Actgcttg                                                                                                                                   <  1:131991/247‑3 (MQ=255)                                                                                                          caGCTCAGTGATGGTGGGGTGGGAAGGAAAGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGAGGGGCTGGGACAAAACAGGCTG‑GCAATACTGCTGGCTGTCTTGCACCGGCACGGTGGTCTGCAACTGGCCGATCAGGATGTGTTTGTGAACG‑TAGTCGTCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTA‑CGTTTACTGCGGGCg                                                                                                                                 >  2:516594/2‑250 (MQ=255)                                                                                                           aGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTtg                                                                                                                                                                                                 <  1:614453/187‑1 (MQ=255)                                                                                                             ctccgTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCAGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCGGGGTGTGTTTGTGAACGGTGG‑CGGCGGCGGGAAGGTAACCGAAAACCGCTACGACTTAGCGTT‑ACTGCTGGGGATgg                                                                                                                             >  2:454887/5‑251 (MQ=255)                                                                                                                   gtttcgtgagtGGGATGGAACGCGTTCGTTGCTTGTGGCGAGTCAGGCGCTGGTCGGTCCCGCGATTAGGGCTCATCCGCGCGGCGGGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCGGCTGGCTGGCTTGCACCGTCACGGTGGGCTGCAAATGGCCGATCAGGAGGTGTCTGTGATCG‑TGGTCGGCGGCGTGAAGTTAACCGAAACCAGCGCAGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCa                                                                                                                       <  2:71762/241‑1 (MQ=255)                                                                                                                   gATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCa                                                                                                                       <  1:641681/251‑1 (MQ=255)                                                                                                                                  aggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAAAAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGACGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGGT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAc                                                                                                        >  2:523349/1‑251 (MQ=255)                                                                                                                                   ggaaCGCGTCCGTTGCTGGTGGAGATTCCAGCGCTGCTCCCTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAACCCGTATG‑GCAATCACGCTGGCTGTCTTGCACCGTCACGGAGGTCTGCAAATGGCCGATCAGGTTGTGTCTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAGCCAGCGCCGGCTTAGCGTT‑ACTGCTGGCGATGGTTTACCGCCTGCGTGGCAGAcc                                                                                                       >  2:660129/1‑251 (MQ=255)                                                                                                                                    gaacgctctctttttgggtggAGATTCAGGTGCTGGTCGATCACTGGATGATGGCGAATCCGCGCCGCGTGGCAGGGGGGCTGGAACAAACCCGTCTG‑GCACTCCTGCTGGCTGTCTTGCACCTTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑GGGTCGGCGGCGTGAAGGTAACCGAAACCAGCTCCGACTAAGCGTT‑TCTGCTGGCGATGGGTTCCAGCCTGCGTGACAGACCg                                                                                                      <  2:410471/235‑1 (MQ=255)                                                                                                                                       cggttccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGAGGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGCCTG‑GCAGTCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTg                                                                                                   <  2:648762/242‑1 (MQ=255)                                                                                                                                          gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGCCTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGACTGCGGGACAGACCGCTGccc                                                                                                >  2:236611/7‑250 (MQ=255)                                                                                                                                          gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGcc                                                                                                 >  1:545516/7‑250 (MQ=255)                                                                                                                                             cgttgctgGTGGAGATTCAGGCGCTGTTCGATCACGCGATGATGGCGAATCCGCGCCTCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAg                                                                                             <  1:97439/248‑1 (MQ=255)                                                                                                                                             cgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑Tgg                                                                                                                                                                                       >  1:332203/4‑162 (MQ=255)                                                                                                                                             cgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑Tgg                                                                                                                                                                                       <  2:332203/159‑1 (MQ=255)                                                                                                                                                tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGa                                                                                           <  1:50588/250‑1 (MQ=255)                                                                                                                                                tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGa                                                                                           <  1:594237/250‑1 (MQ=255)                                                                                                                                                     gtggAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGTTGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCtggt                                                                                     <  1:660129/251‑1 (MQ=255)                                                                                                                                                            ttCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGttt                                                                                                                                                                                                  >  1:319157/1‑137 (MQ=255)                                                                                                                                                            ttCAGGCGCTGGTCGATCACGCGATGATGGCGAAGCCGCGCCGCGTGGCAGGTGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCGGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGttt                                                                                                                                                                                                  <  2:319157/137‑1 (MQ=255)                                                                                                                                                             tCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGttt                                                                              <  1:584484/250‑1 (MQ=255)                                                                                                                                                                           aTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTg                                                               <  1:313626/251‑1 (MQ=255)                                                                                                                                                                              ccTCGATGATGGGGAATCCGCGCCGCGTGGCAGTGGGGCGGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCa                                                            <  2:438021/250‑1 (MQ=255)                                                                                                                                                                                  gatgatGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGACCAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTACGTGACAGGCCGCTGCCACCGGATATGGTGGTGTTTGGtagaggagggctggcggggg                                                        >  2:675795/1‑232 (MQ=255)                                                                                                                                                                                   atgatgGCGAATCCGCGCCGCGTGGCAGAGGGGCTGGCACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGGAACCGAAACCAGCGCCGACTTAGCGTT‑CCTGCTGGCGATGGTTTCCAGCCTGCGTGGCAGACCGCTGCCACAGGTTCTGGTGGGGTTTTGTGAAGTCGGGcgggacggggg                                                       >  2:31547/1‑241 (MQ=255)                                                                                                                                                                                    tgatgGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAgggg                                                        >  1:428578/1‑249 (MQ=255)                                                                                                                                                                                       gggCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATc                                                   <  1:323335/250‑1 (MQ=255)                                                                                                                                                                                            aatccgcgccgcGTGGCAGTGGAGATGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGCACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACCGACCGCTGCAACAGGTTCTGGTGGTGTTTGGGGAGGTCGGGCTGGCGGGGGAGATCCGccc                                              >  2:465289/7‑251 (MQ=255)                                                                                                                                                                                               ccgcgccgcGTGGCAGTGGGGCTGGGACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCTCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGt                                           <  1:606953/248‑1 (MQ=255)                                                                                                                                                                                                                ggggCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCGCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAAc                          <  1:646331/251‑1 (MQ=255)                                                                                                                                                                                                                    cTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAt                      >  1:83922/1‑251 (MQ=255)                                                                                                                                                                                                                      ggAACAAAACCG‑GTGTGCAATCCGGGTGGCTGTTTTGCTCCGTCACGGTGGTCTGCAAATGGCCGAGCAGGATGTGTTTGTGAACG‑TGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACGTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAt                      <  2:10935/249‑1 (MQ=255)                                                                                                                                                                                                                       gAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTCGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGACGGGGGGATCCGCCCGGCGCCCAGCgggatggtcagcaatgc                   >  2:446311/1‑236 (MQ=255)                                                                                                                                                                                                                       gAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttc                   <  1:651861/251‑4 (MQ=255)                                                                                                                                                                                                                           aaaaCCGTCGG‑GCAAACCTGCTGGCTGTCTGGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttcagaa               <  2:428578/251‑8 (MQ=255)                                                                                                                                                                                                                                      ggaaattctgctgTCTGTCTGGCACCGGCACGGGGGGCTGCAAATGGCCGATCAGGATGTGTTGGTGAACG‑TGGGCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTT‑ACTGCTGGCGTTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCACAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAac     <  2:272108/244‑1 (MQ=255)                                                                                                                                                                                                                                         aagccTGCTGGCTTTCTTGCCCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACg  <  2:386556/248‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              CGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACG  >  NC_000913/4626697‑4627178

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