Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,069	C→T	R386R (CGC→CGT)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,069	0	C	T	94.9%	103.0 / ‑4.3	40	R386R (CGC→CGT)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (14/23);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATG  >  NC_000913/4626827‑4627317                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         ggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAcc                                                                                                                                                                                                                                                  <  1:128710/251‑1 (MQ=255)     cgcgTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTg                                                                                                                                                                                                                                              <  1:258667/251‑1 (MQ=255)              tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGAt                                                                                                                                                                                                                                     >  1:80088/1‑251 (MQ=255)               gctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATc                                                                                                                                                                                                                                    <  1:259842/251‑1 (MQ=255)               gctgGTAGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGAGTGGCAGCGGGGCGGGACCACAAACGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTCGGCTGCAAATGACCGATCAGGATGTGTATGTGAACGTGGTCGACGGCGGCAA‑GGTAACCGAAACCAGCGCCGACTTAGCGTAACTACTAGCGATGGTTTCCAGACTGCGCGACAGACCgatggcacaggccc                                                                                                                                                                                                                                    >  2:116860/1‑238 (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)                                                                                                                                                                                                                                        ttCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTTTAATGAGGTATACgg             >  1:373951/1‑250 (MQ=255)                                                                                                                                                                                                                                           cAGCCTGCGTGACAGACCGCGGCCACAGGATCTGGTGGTGTTTGGTGAAGGCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACggaggg         <  2:77316/251‑1 (MQ=255)                                                                                                                                                                                                                                                  cgcgACAGCCCGCGGCCGCAGGATCTGGTGGTGTTTGGTGAAGGCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCACTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGTTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGATTAGCGTGTTCGACGACTTATAAGGAGATATACGGAGGGAGATATg  <  2:251894/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         GGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAA‑GGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATG  >  NC_000913/4626827‑4627317

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