Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,886	C→T	N325N (AAC→AAT)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,886	0	C	T	97.7%	129.5 / NA	44	N325N (AAC→AAT)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (23/20);  minor base G (0/1);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTC  >  NC_000913/4626644‑4627111                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  ggCGAGGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGcgcc                                                                                                                                                                                                                            <  2:204839/250‑1 (MQ=255)   cGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATccgcgccgcg                                                                                                                                                                                                                         >  1:13818/1‑251 (MQ=255)    gATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAGTCCGCGCCGCGt                                                                                                                                                                                                                        >  1:76199/1‑251 (MQ=255)         cGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGCAGTGACCTCCGGCAGCTCCGTGCTGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGCGGc                                                                                                                                                                                                                     >  2:228019/1‑249 (MQ=255)              cccGTTTTCGCTCTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCATCCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGGACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTgggg                                                                                                                                                                                                              <  1:39796/251‑1 (MQ=255)              cccGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGATCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTggg                                                                                                                                                                                                               <  1:190543/250‑1 (MQ=255)                  ttttCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTgg                                                                                                                                                                                                          >  1:64895/1‑251 (MQ=255)                         actttacgcAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACaaa                                                                                                                                                                                                    >  2:98657/7‑250 (MQ=255)                         actttacgcAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACaaa                                                                                                                                                                                                    >  2:82887/7‑250 (MQ=255)                          ctttacgcAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTCAGCGTCTTCGCGATGAACGAACAGGGGCTGCGTGAAATCAACAACCCTTCGGCAATTTTCTTACGTCGCGGAGATGAAGTGACATCCGGCAGCTCCGTGATGGTGGTGGGGGAAGGAACGAGTCCGTTGCTGGTGGCGATTCAGGCGCTGGTCGATCACTAGATGATGGCGCATCCGCGCCGCGTGGCAGTGGGGCTGGACCAAAAc                                                                                                                                                                                                  >  1:20572/6‑251 (MQ=255)                                gcAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGCACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGCGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTAAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGATGGAACACAACCGTCTg                                                                                                                                                                                            >  2:5262/1‑251 (MQ=255)                                                       gcgGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGGCGATCACTCGAGGAGGGCGAATCCGCGCCGCGTGGCAGGGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCa                                                                                                                                                                     <  1:128095/251‑1 (MQ=255)                                                              aTGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCAc                                                                                                                                                              >  2:27971/1‑251 (MQ=255)                                                              aTGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGGCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGCGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCAc                                                                                                                                                              <  2:10424/251‑1 (MQ=255)                                                                 aGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACggt                                                                                                                                                           <  1:179730/251‑1 (MQ=255)                                                                   cTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGAGGGGCTGGAACAAAACCGTCGGGCAATACTGCTGGCGGTCTTGCACCGTCACggggg                                                                                                                                                         >  2:160873/1‑248 (MQ=255)                                                                        cGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGCAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGACCAAAACCGTCTGGCACTCCTGCTGGCTGTCTTGCCCCTTCACGTGGGTCTGc                                                                                                                                                    >  1:165010/1‑251 (MQ=255)                                                                         gTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCGTCGGCAATTTGCTGAAGTCGCGGAGATGAAGTGACCGCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGTCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGTTGTCTTGCACCGTCACGTTGGTCTGCa                                                                                                                                                   <  2:226061/251‑1 (MQ=255)                                                                           cTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGAAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGCGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCaa                                                                                                                                                  >  1:34466/1‑250 (MQ=255)                                                                               gcgATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGc                                                                                                                                             >  2:96571/1‑251 (MQ=255)                                                                                       cGAACAGCGGCTGCGTGAAGTCAGCAAACCTTCGGCAATTTTCTTAAGTCGCGGAGCTGCAGTGACCTCCGGCAGCTCCGTGATGGGGGCGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAgg                                                                                                                                     >  2:157466/1‑251 (MQ=255)                                                                                                               cAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGATCCGTGATGGTGGTGCGGGAAGGAACGCGTCCGTTGCTGGTGGCGATTCAGGCGCAGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTcggcg                                                                                                             >  1:183409/1‑251 (MQ=255)                                                                                                                 aCCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGAGTTTGTGAACGTGGTcggcgg                                                                                                            >  1:174006/1‑250 (MQ=255)                                                                                                                  cccTTCGGCAAGTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCGCCGTGATGGGGGTGTGGGAAGGAACTCGTCCGTTGCTGGTGGAGAGTCAGGCGCTGGTCGATCACGCGATGATGGCTAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGgcggcg                                                                                                          <  1:141663/251‑1 (MQ=255)                                                                                                                    cTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGCACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTg                                                                                                        >  2:103437/1‑251 (MQ=255)                                                                                                                           aaTTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTC‑AGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAAc                                                                                                <  2:213291/251‑1 (MQ=255)                                                                                                                                        cgcgGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAgcgc                                                                                    >  1:252294/1‑248 (MQ=255)                                                                                                                                            gagATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGCGATTCAGGCGCTGGTCGATCACTCGATGATGGCGACTCCGCGCCGCGTGGCAGTGGGGCTGGCACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGAc                                                                                >  2:128559/3‑251 (MQ=255)                                                                                                                                             agATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACt                                                                               <  1:82887/250‑1 (MQ=255)                                                                                                                                               aTGAAGTGACCGCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGGTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTa                                                                             <  2:17141/251‑1 (MQ=255)                                                                                                                                                      gACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGGCCGTTGCGGGTGGAGATGCAGGCGATGGGCTATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTATCGTTa                                                                       <  1:49111/250‑1 (MQ=255)                                                                                                                                                           ccggcagcgccgTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTActgctg                                                                 <  2:211858/240‑1 (MQ=255)                                                                                                                                                                        gATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGTATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGTAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCa                                                    <  1:160873/251‑1 (MQ=255)                                                                                                                                                                         aTGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAg                                                   <  1:27971/251‑1 (MQ=255)                                                                                                                                                                            gtggtgTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCt                                                <  1:98657/251‑1 (MQ=255)                                                                                                                                                                                 gtGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAGCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTgcgtg                                           <  1:47408/250‑3 (MQ=255)                                                                                                                                                                                            cgcggccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGATCGCTg                                <  2:13818/242‑1 (MQ=255)                                                                                                                                                                                               gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCCAGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGccc                             >  2:234857/7‑250 (MQ=255)                                                                                                                                                                                                tccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGCTGGCGCATCAGCGAAGCGTGGCAGTGGCGCTGGCACAAAACCGTCTGGCACTCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCCAATGGCCGAACAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGAAACCGGAACCAGCGCCGACTTAGCGGTAATGCTGGCGATGGTTTCCAGCCTGCGTGACCGGCCGCTGccac                            >  2:242046/6‑249 (MQ=255)                                                                                                                                                                                                 cagttgctgGTGCAGATTCAGGCGAGCGTCGAGCACTCGATGATTGCGAATCCGCGCCACGTGGCAGTTGGGCTGGAACAACACCGTCGTGCAAGCCGGCTGGCTGTCGTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTGGGTGGCGGGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTCCTGGCGATGGTTTCCAGCCTGCTTGACAGACCGCTgactc                            <  2:16063/246‑5 (MQ=255)                                                                                                                                                                                                   gttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTATTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAgg                         >  1:152773/3‑251 (MQ=255)                                                                                                                                                                                                       ctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGAGGTCTGCAAATGGCCGATCAGGATGGGTTTGTGACCGTGTTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGAGTTACTGCTGGCGGTGGTTTCAGGCCTGCGGGACAGACCGCTGCCCCAGGGTCt                     >  2:42788/1‑251 (MQ=255)                                                                                                                                                                                                             gagaTTCAGGCGCTGGTCGATCACTCGATGATGGCGAAGCCGCGCCGCGTGGCAGTGGGGCGGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTggtggt               <  2:200690/251‑1 (MQ=255)                                                                                                                                                                                                                           ggTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTc  <  1:259074/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  GGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTC  >  NC_000913/4626644‑4627111

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