Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,889	A→G	P326P (CCA→CCG)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,889	0	A	G	92.5%	98.0 / 0.4	40	P326P (CCA→CCG)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (17/20);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-01

GGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTG‑ATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCC  >  NC_000913/4626644‑4627130                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ggCGGTGCCGACTCCCGTTTGCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTGTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGcgccg                                                                                                                                                                                                                                               <  2:179219/251‑1 (MQ=255) ggCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGcgccg                                                                                                                                                                                                                                               >  1:300843/1‑251 (MQ=255)         cGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGc                                                                                                                                                                                                                                         >  1:218609/1‑249 (MQ=255)            cTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGACTGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGCAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGCATCCGCGCCGCGTGGCAGTgg                                                                                                                                                                                                                                    >  1:169965/1‑251 (MQ=255)            cTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGATGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTACGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGAGCCGCGTGGCAgggg                                                                                                                                                                                                                                    >  1:37820/1‑248 (MQ=255)                       ggacgttacgcgtccatgaatgacGCTTCGGCGCGGTGAATTAGCCGGGCGTCCGCGCGATGTCCGTGCAGGGGCTGCGGGAAGTCATCAACTCATCGGCATTTTTCTTCTGTCGCGGAGATGAAGTGACCTCCGGCTGGTC‑GTGGATGTTGGTGTTGGAAGGGGCGCGTCCGTTGCTGGGTGAGATTCAGGCGCTGGTCGATCACTCGATTATTGCGAATCCGCTCCGCGTGGCAGTGGGGCTGGAACaa                                                                                                                                                                                                                         <  1:22350/228‑1 (MQ=255)                             taccctgccagaaaacgcGCTTCGGCGCGGTGAATGAGCTGTGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTGCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACGCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCg                                                                                                                                                                                                                    <  2:218609/233‑1 (MQ=255)                                  aGCCATAAAAACCGCTTCGGCGCGGTGAATGCGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGAAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGc                                                                                                                                                                                                              <  2:47472/251‑1 (MQ=255)                                      aTAAAAACCGCTTCGGCGAGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCTTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGCCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATc                                                                                                                                                                                                          <  1:224787/251‑1 (MQ=255)                                       tAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGGCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGAGGGCGACTCCGCGCCGCGTGGCGGTGGGGCTGGAACAAAACCGTCTGGCAATcc                                                                                                                                                                                                         >  2:326363/1‑251 (MQ=255)                                                  tCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCAGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCtgt                                                                                                                                                                                               <  1:202647/250‑1 (MQ=255)                                                      cgcgGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTTGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGAGTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGAAATCCTGCTGGCTGTCTTGc                                                                                                                                                                                          <  2:263840/251‑1 (MQ=255)                                                         ggTGAATGAGCTGGGCGTCTTCTCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTGCTGAAGTCGCGGAGATGAAGTGACCTCCGGCTGTTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGTTGGCTGTCTTGCAcc                                                                                                                                                                                       <  2:275427/251‑1 (MQ=255)                                                             aaTGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCa                                                                                                                                                                                   <  1:287635/251‑1 (MQ=255)                                                              aTGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCAc                                                                                                                                                                                  >  1:168623/1‑251 (MQ=255)                                                                    tGGGCGTCTTCGCGATAACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGACGCGGAGATGAAGAGACCTCCGGCAGCTCCGTG‑ATGGAGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGGGATTAAGGCGCTGGGCGATCACTCGATGATGGCGAATCCGAGCCGCGTGGCAGTGGGGCTGGACCAAAACCGTCTGGACATCCTGATGGCTGTCTTGACCCGTCACggtgt                                                                                                                                                                             >  2:334276/1‑249 (MQ=255)                                                                         gTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGGCGCGGAGATGAAGTGACCTCCGGCAGTTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCa                                                                                                                                                                       <  1:260807/251‑1 (MQ=255)                                                                                     aCCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCGCTGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGTCAATCCTGCTGGCTGTCTTGCTCCGTCACGGTGGTCTGCAAATGGCCGTTCa                                                                                                                                                           <  1:136213/251‑1 (MQ=255)                                                                                          aCAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGAtg                                                                                                                                                      >  1:77316/1‑251 (MQ=255)                                                                                               ggCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGttt                                                                                                                                                 <  1:201249/251‑1 (MQ=255)                                                                                                   gCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTtgtg                                                                                                                                              <  1:326363/250‑1 (MQ=255)                                                                                                                         gCAATTTTCTTAAGTCGCGGAGATGAAGTGCCCTCCGGCATCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTTCTGGTGGAGATTCAGGCTCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGTGGCTGGACCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGt                                                                                                                       <  1:14319/251‑1 (MQ=255)                                                                                                                             ttttCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGGATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGAAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAAcc                                                                                                                   >  1:163267/1‑251 (MQ=255)                                                                                                                                           ggcgATGAAGTGACCTCCGGTAGCTCAGTG‑ATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTgggg                                                                                                                                                                                                                                  >  1:16246/1‑126 (MQ=255)                                                                                                                                           ggcgATGAAGTGACCTCCGGTAGCTCAGTG‑ATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTgggg                                                                                                                                                                                                                                  <  2:16246/126‑1 (MQ=255)                                                                                                                                                 gAAGTGACCTACGGCAGCTCCGTG‑ATGGTTGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAAGCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCGGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGc                                                                                              <  2:168623/251‑1 (MQ=255)                                                                                                                                                    gTGACCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGtt                                                                                           >  1:129337/1‑251 (MQ=255)                                                                                                                                                       aCCTCCGGCAGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTAct                                                                                        >  1:24155/1‑251 (MQ=255)                                                                                                                                                            cggcaGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGGCTTGCAACGTCCCGGAGGGCTGCAAATGGCCGATCAGGAGGTGTTTGGTAAAGTGGTCGGCGGCGTGTAGGTAACCACAACACA‑AGCCGGCTTAGCGATACTGCTgg                                                                                   >  2:10125/5‑251 (MQ=255)                                                                                                                                                               caGCTCCGTG‑ATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGACGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTACTGCTGGCgg                                                                                >  1:232332/2‑250 (MQ=255)                                                                                                                                                                                         ggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAcc                                                       <  1:128710/251‑1 (MQ=255)                                                                                                                                                                                             cgcgTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTg                                                   <  1:258667/251‑1 (MQ=255)                                                                                                                                                                                                      tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGAt                                          >  1:80088/1‑251 (MQ=255)                                                                                                                                                                                                       gctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATc                                         <  1:259842/251‑1 (MQ=255)                                                                                                                                                                                                       gctgGTAGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGAGTGGCAGCGGGGCGGGACCACAAACGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTCGGCTGCAAATGACCGATCAGGATGTGTATGTGAACGTGGTCGACGGCGGCAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTAACTACTAGCGATGGTTTCCAGACTGCGCGACAGACCgatggcacaggccc                                         >  2:116860/1‑238 (MQ=255)                                                                                                                                                                                                           gtggAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATACGCGCCGCGGGGCAGTGGGGCTGGAACAAAACCGACTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAAAGGCCGATCAGGATGTGTTTGAGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTACTGCTGGCGCTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGAACtgg                                      >  1:374905/1‑250 (MQ=255)                                                                                                                                                                                                              gagaTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTACTGCTGGCGAAGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCtggtggg                                  >  1:180558/1‑250 (MQ=255)                                                                                                                                                                                                                                gATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCt                <  1:278325/251‑1 (MQ=255)                                                                                                                                                                                                                                     cTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCATGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGTGCCGACTTAGCGTTACTGCTGGAGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGc             >  1:251894/1‑249 (MQ=255)                                                                                                                                                                                                                                              ggCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATcc  <  1:10125/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTG‑ATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAC‑CAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCC  >  NC_000913/4626644‑4627130

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