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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,084	0	G	A	92.9%	107.3 / ‑3.1	43	P391P (CCG→CCA)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (17/22);  total (18/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.96e-01

TGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑ATATACGGAGGGAGATATGTCGTCATTTGAT  >  NC_000913/4626840‑4627329                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGAt                                                                                                                                                                                                                                                  >  1:80088/1‑251 (MQ=255)  gctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATc                                                                                                                                                                                                                                                 <  1:259842/251‑1 (MQ=255)      gtggAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATACGCGCCGCGGGGCAGTGGGGCTGGAACAAAACCGACTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAAAGGCCGATCAGGATGTGTTTGAGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGCTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGAACtgg                                                                                                                                                                                                                                              >  1:374905/1‑250 (MQ=255)         gagaTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGAAGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCtggtggg                                                                                                                                                                                                                                          >  1:180558/1‑250 (MQ=255)                           gATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCt                                                                                                                                                                                                                        <  1:278325/251‑1 (MQ=255)                                cTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCATGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGAGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGc                                                                                                                                                                                                                     >  1:251894/1‑249 (MQ=255)                                         ggCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATcc                                                                                                                                                                                                          <  1:10125/251‑1 (MQ=255)                                                      cgcgTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAg                                                                                                                                                                                             <  1:314535/251‑1 (MQ=255)                                                         gTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCgg                                                                                                                                                                                          >  1:149417/1‑251 (MQ=255)                                                          tGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGACACCAGCGCCGACTTAGCGTTGCTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGGGGGGGTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCACAGCGGt                                                                                                                                                                                         >  1:63024/1‑251 (MQ=255)                                                                 ggggCTGGAACAAAACCGGCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGaa                                                                                                                                                                                   <  1:337782/250‑1 (MQ=255)                                                                     cTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATAAGGCTATGTTTGTGAACTAGGTCGGCGGCGTGAA‑GGTAACCGCAACCAGCGACGACTTAGCGATACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCGGCCACAGGATCTGGAGGTGTTTGGTGAATTCGGGCTGGCCGGGGAGATCCGCCAGGTGCCCAGCGGTCAGGCACgca                                                                                                                                                                               >  2:371700/1‑248 (MQ=255)                                                                       ggAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtt                                                                                                                                                                             <  1:77150/250‑2 (MQ=255)                                                                         aaCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttca                                                                                                                                                                          <  1:362152/251‑5 (MQ=255)                                                                         aaCAAAACCGGCTGGCAATCCTGCTGGCGGTCTTGCACCGTCACGGTGGTCTGCAACTGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGGGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttca                                                                                                                                                                          <  1:382317/251‑5 (MQ=255)                                                                             aaaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAg                                                                                                                                                                      <  1:265240/251‑1 (MQ=255)                                                                                   tCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGTTGCGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACACGCGCCGACTTCGCGTTACTGCTGGCGCTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGTTCTGGCGGTGTTTGGTGAAGTCGGGCTGGCAGGGCAGATCCGCCCGGTGCCCCGTGGCCAGGCACGCATTTCAGAAgcggcg                                                                                                                                                                 >  1:54870/1‑250 (MQ=255)                                                                                            ccTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGtt                                                                                                                                                       >  1:191799/1‑251 (MQ=255)                                                                                                      gtcttGCACCGTCACGGTGGTCGGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGGGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGTCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCgg                                                                                                                                               <  1:350029/246‑1 (MQ=255)                                                                                                                             gccatgggacggtcggattggttttGTGAACGTGGTTGGCGG‑ATGAATGGTAACCAAATCGAGGCCCGGCTTGGCGTTAGTGCGGGCAATGGTTTCCGGCCTGCGGTCCAGTCCGCTTTCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCCGGGAGGGGAGATCCGCCCGGTGCCCAGCGGTGAGGAACGAATTTCAGAAGCGGCGAAACACCGTTTTCGCCGGGTGATTGTTCCGGCGGCTAACGtac                                                                                                                      <  2:3608/229‑3 (MQ=255)                                                                                                                                   ggctgATCATGATGTGTTGGTGACCGTGTTCGGCGGCTGGAG‑GGTACACGAAACCAGCGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGCCAGACCGCTGCCCCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGACCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCTATTGTTCCGGCGGCTAACGTACCGaaa                                                                                                                 <  2:301318/246‑1 (MQ=255)                                                                                                                                   ggCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGaaaa                                                                                                                <  1:301528/251‑1 (MQ=255)                                                                                                                                       gATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgc                                                                                                            >  1:379366/1‑251 (MQ=255)                                                                                                                                         tCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgcgc                                                                                                          >  1:188408/1‑251 (MQ=255)                                                                                                                                                 tgtTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCCGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCAACAGGATCTGGGGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAAATCCGCCCGGAGCCCAGCGGTCAGGCACGACTTTCAGAAGCGGCGAAACACGGTTTTCGCCGGACGATTGTACCGGCGGCTaaagaacagaaaaatgcgcaggacggg                                                                                                  >  1:307450/1‑226 (MQ=255)                                                                                                                                                   tttGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAAc                                                                                                                                                                                  >  1:80629/1‑169 (MQ=255)                                                                                                                                                   tttGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCGGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCTGTGCCCAGCGGTCAGGAAc                                                                                                                                                                                  <  2:80629/169‑1 (MQ=255)                                                                                                                                                     tgtgAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATg                                                                                               <  1:151472/250‑1 (MQ=255)                                                                                                                                                                      gcgTGAA‑GGTAACTGAACTCAGCGCCGCCTTTGCGTTACTGCTGGCGATGTTGTCCAGCCTGCGTGACAGCCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTGCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGTATGCAGATTTTTGGCGTTaaa                                                                             <  2:129337/251‑1 (MQ=255)                                                                                                                                                                       ggTGAA‑GGCAACCGAACACAGCGCCGCCTTTGCGTGTCTGCTTGCGATGGTTTCCAGCCTTCGTGACATACCGCTGACACAGGATCTGGGGGTGTTTGGGGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTaaaa                                                                            <  2:300843/250‑1 (MQ=255)                                                                                                                                                                               gTAACCGAAACCAGCGCCGATTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACtctc                                                                     <  2:191799/251‑1 (MQ=255)                                                                                                                                                                                      aaaCCAGCGCCGCCTTAGTGTTACTGTTGGCGATGGTGTCCAGCCTGCGTGACAGACCGCTGCCACAGGATTTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGCATTTCAGAAGCGGCGAAACTCGGTTTTCGCCGGGCGATTGTGCCGGCGGCTAACTTGCCTAAAAAAGTGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAcgcg                                                              <  1:205088/251‑1 (MQ=255)                                                                                                                                                                                                 gACTTTGCGTTACTGCTGGCGATGGTTTCCAGCCTTCGTGACAGGCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGtt                                                   <  1:312889/251‑1 (MQ=255)                                                                                                                                                                                                         gTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACt                                           <  1:72450/251‑1 (MQ=255)                                                                                                                                                                                                          ttACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAACACTCTCCGAAGAGCTTAGCGTGTGCGACGACtt                                          >  2:51748/1‑251 (MQ=255)                                                                                                                                                                                                                           ttCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTTTAATGAG‑GTATACgg                         >  1:373951/1‑250 (MQ=255)                                                                                                                                                                                                                              cAGCCTGCGTGACAGACCGCGGCCACAGGATCTGGTGGTGTTTGGTGAAGGCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACggaggg                     <  2:77316/251‑1 (MQ=255)                                                                                                                                                                                                                                    gcgtgACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGGGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGGAGGGATGCAGGTTTTTGGTGTTAAAGAACTCTCCGTCGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GCATCCGGAGGGAGatat               >  1:177934/5‑251 (MQ=255)                                                                                                                                                                                                                                     cgcgACAGCCCGCGGCCGCAGGATCTGGTGGTGTTTGGTGAAGGCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCACTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGTTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGATTAGCGTGTTCGACGACTTATAAGGAG‑ATATACGGAGGGAGATATg              <  2:251894/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATgtcgt          <  1:65750/251‑1 (MQ=255)                                                                                                                                                                                                                                           tgACCGCTGCCACAGGATCTGGTGGAGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGGAACCCGGTTTTCGCCGGGCGCTATTTCCGGCGGCTAAAGGACCGAAAACAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTTTCCGACGCGCTTAACGTGTTAGACGCCTTATAATGAGGAT‑TACGGGGGGAGATATgtcgtct        >  2:341315/2‑250 (MQ=255)                                                                                                                                                                                                                                             aCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCAtt      <  1:201289/251‑1 (MQ=255)                                                                                                                                                                                                                                                 cTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGGAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGCGATATGTCGTCATTTGAt  >  2:326865/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       TGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑ATATACGGAGGGAGATATGTCGTCATTTGAT  >  NC_000913/4626840‑4627329

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