Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,213	G→A	V434V (GTG→GTA)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,213	0	G	A	97.8%	135.1 / NA	45	V434V (GTG→GTA)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/28);  minor base T (0/1);  total (16/29)
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.75e-01

AAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTG  >  NC_000913/4626966‑4627458                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      aaaTGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCg                                                                                                                                                                                                                                                    <  2:217413/251‑1 (MQ=255)            cAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGC‑ACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGAGCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGGTTGTTCCGGCGGCTAACGTTCCGAAAAAAGCGCCg                                                                                                                                                                                                                                        <  1:76894/251‑1 (MQ=255)             aGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCg                                                                                                                                                                                                                                        <  1:244410/251‑1 (MQ=255)                     tttGTAAAAGTGGTCGGCGGCGTGCAGGTAACCGCAAACAGAGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGACACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTACAGAAGCGGCGAAACACGGCTTTAGCCGGGAGATAGTTCCGGCGACTAACGTAACGAAAAAAgacacggaaggggt                                                                                                                                                                                                                                >  2:186850/1‑238 (MQ=255)                            ccGTGGTCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAtt                                                                                                                                                                                                                         <  2:166032/250‑1 (MQ=255)                               tGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAttttt                                                                                                                                                                                                                      >  1:79245/1‑251 (MQ=255)                                gggcggcggCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTGACTGCTGGCGATGGTTTCCAGCCTGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGGTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGTCCGGGGGCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTGCGGCGTCTAACGGACCTAAAAAAGCGCCGGAAGGGATGCAGAttgtt                                                                                                                                                                                                                      <  2:243160/247‑4 (MQ=255)                                  gcggcggCTTGAATGTAACCGAATCCTGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTTACAGACCTCTGCCACAGGATCTGGTGGTGTGTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCATATTTTTGGc                                                                                                                                                                                                                   <  1:96942/250‑1 (MQ=255)                                    ggcggcGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGt                                                                                                                                                                                                                 >  2:165583/1‑251 (MQ=255)                                    ggcggcGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGt                                                                                                                                                                                                                 >  1:57114/1‑251 (MQ=255)                                          gTGAGGGTAACCGAATCCAGCGCCGACTTAGCGTTACTGCTGGCGAGGGTTTCCAGGCTGCGTGACGGACCGCTGCCACAGGAGCTGGTGGTGTTTGGTGAAGTCGTGTTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAGTTCAGAAGCGGCGAAACACGGTTTTCGCGGGGGGATTGTTCCGGCGGCGAACGTACCGAAAAAATCGCCGGAAGGGATGCAGAGTTGTGGCGGTaaaaa                                                                                                                                                                                                           <  2:34466/251‑1 (MQ=255)                                               ggTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTAGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTAAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACAGAAAAAAGCGCCGGAAGGGATGCAGAATTTTGGAGTTAAAAAActct                                                                                                                                                                                                      >  2:72463/1‑251 (MQ=255)                                                  aaCCGAAACCCGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGACACAGGATCGGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGCTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCg                                                                                                                                                                                                   <  2:152773/251‑1 (MQ=255)                                                  aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGTTCAGGAACGAATTTCAGAATCGGCGAAACACGGTTTGCGCCGGGGGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCg                                                                                                                                                                                                   <  1:242046/251‑1 (MQ=255)                                                     cGAAACCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAATCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAcg                                                                                                                                                                                                <  1:228019/251‑1 (MQ=255)                                                     cGAAACAAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAACGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAcg                                                                                                                                                                                                >  2:235801/1‑251 (MQ=255)                                                              cgCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCg                                                                                                                                                                                       <  2:115028/250‑1 (MQ=255)                                                                 cGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCgtgt                                                                                                                                                                                    <  1:234857/251‑1 (MQ=255)                                                                 cGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCgtgt                                                                                                                                                                                    >  1:91523/1‑251 (MQ=255)                                                                 cGACTGAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCgtgt                                                                                                                                                                                    <  1:172839/251‑1 (MQ=255)                                                                  gACTTAGCGTGACTGCTGGCGAGGGTTTCCAGCCTGCGTGACAGTCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGtt                                                                                                                                                                                   <  2:76199/251‑1 (MQ=255)                                                                      tAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTcg                                                                                                                                                                                 <  1:128559/249‑1 (MQ=255)                                                                       aGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGCCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGTCTGGCAGGGTAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTcgacg                                                                                                                                                                              <  1:120068/251‑1 (MQ=255)                                                                           ttACTGCTGGCGAGGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCgacgac                                                                                                                                                                            <  2:10762/249‑1 (MQ=255)                                                                                            ttCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGTGATGCATATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCTTGTTCGACGACATATAAttgggaatacgg                                                                                                                                                          <  2:183409/250‑12 (MQ=255)                                                                                                  ccTGCGTGACAGACCTGTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGgaga                                                                                                                                                   <  1:103437/251‑1 (MQ=255)                                                                                                    tGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGat                                                                                                                                                  >  2:29019/1‑250 (MQ=255)                                                                                                               ccGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCCGGGGAGATCCGCCCGGTGCCCCGCGGTAAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCAGGCGGCTAACGTACCGACAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTAGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTAAttt                                                                                                                                      >  2:120269/1‑251 (MQ=255)                                                                                                                        cAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTg                                                                                                                             >  1:164934/1‑251 (MQ=255)                                                                                                                                   ggtgTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCa                                                                                                                   <  1:157466/250‑1 (MQ=255)                                                                                                                                                   gggcgggcAGGGGAGATCCGCCCGGTGCCCAGCGGGCAGGAACTAATGTCAGAAGCGGCGAAACACGGTTTGCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATGTTTGGCGTTAAAAAACTCTCCGACGCTCTTAGCGTGTTCGACGACTTACAATGAGGAATACggaggg                                                                                                                                                      <  2:232327/199‑1 (MQ=255)                                                                                                                                                   gggCTGGCAGGGGAGATCCGCCAGGTGCCCAGCGGTCAGGAACGCATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAGGGGATGCAGATTTTTGGCGTTAAACAACTCTCCGACGCGCTTAGCGTGTTAGACGACTTGTAATGAGGAATACggaggg                                                                                                                                                      >  1:232327/1‑199 (MQ=255)                                                                                                                                                              ggAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTAcagca                                                                                       <  2:57114/251‑1 (MQ=255)                                                                                                                                                                                 aGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGt                                                                    <  1:165583/251‑1 (MQ=255)                                                                                                                                                                                         ggcaCGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGCCGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCa                                                             >  2:70696/4‑250 (MQ=255)                                                                                                                                                                                                           ggcgAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGtt                                          <  1:235801/251‑1 (MQ=255)                                                                                                                                                                                                               aaaCACGGTTTTCGCCGGGCGATTGATCCGGCGGCTACCGTACCGAAAAAAGCGCCGGAAGGGATGCAGCTTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACATGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCGGTTACTg                                      >  2:251701/1‑251 (MQ=255)                                                                                                                                                                                                                   acGGTTTTCGCCGGGCGATTGGTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGGTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACGTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATg                                  <  2:131679/251‑1 (MQ=255)                                                                                                                                                                                                                      gTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATTACCAAAGGGTATTTAAGCCAGTTACTGAATGcc                                >  2:80275/1‑250 (MQ=255)                                                                                                                                                                                                                           cGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGCCCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATc                          <  1:5262/251‑1 (MQ=255)                                                                                                                                                                                                                                 gCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGc                    >  1:62854/1‑251 (MQ=255)                                                                                                                                                                                                                                      tGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGGAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGGCTTTAAGCCAGTTACTGAATGCCAAAATCAAAAGccca                 >  2:6307/1‑248 (MQ=255)                                                                                                                                                                                                                                                gcTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAg     <  2:216232/251‑1 (MQ=255)                                                                                                                                                                                                                                                   aaCGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTg  <  1:42788/251‑1 (MQ=255)                                                                                                                                                                                                                                                   aaCGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTg  <  2:79245/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      AAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTG  >  NC_000913/4626966‑4627458

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