Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,154	T→C	T184T (ACT→ACC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,154	0	T	C	98.6%	220.3 / ‑6.4	69	T184T (ACT→ACC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (32/36);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAG‑AAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAC  >  NC_000913/4619918‑4620388                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           gCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATa                                                                                                                                                                                                                               >  2:158806/1‑250 (MQ=255)   cGGTGTACCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACTGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAAc                                                                                                                                                                                                                            >  1:480110/1‑251 (MQ=255)      tgtCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAActct                                                                                                                                                                                                                         >  2:139523/1‑251 (MQ=255)      tgtCCCGGATCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGCGCACATGAAAACCGGCAAGCCGCTTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAACctct                                                                                                                                                                                                                         >  2:287228/1‑251 (MQ=255)               tCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACtgt                                                                                                                                                                                                                >  1:546333/1‑251 (MQ=255)                  gTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACtgt                                                                                                                                                                                                                <  1:343113/248‑1 (MQ=255)                    ttGAGTGGGGATATTTCCCCGATCACGAAACCAGCTTCCCGCAAGAGCTGCTGGATAAACGGGTCGAACGCGCAAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGa                                                                                                                                                                                                            <  1:118039/250‑1 (MQ=255)                      gAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACGCTACGAACTGTGCGaaa                                                                                                                                                                                                          <  1:169018/250‑1 (MQ=255)                       aGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATc                                                                                                                                                                                                        <  1:158806/251‑1 (MQ=255)                           gggATATTTCTCCTATCACGAAAACATCTTCCCGGAAAGGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGTTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCTAAATCGccc                                                                                                                                                                                                    <  2:310014/251‑1 (MQ=255)                            ggATATTTCTCCGATCAAGAAAAAAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTATGGCCTGGATAAACTCTACGAACTGTGCGAACTCGCCCg                                                                                                                                                                                                   >  2:129480/1‑251 (MQ=255)                                atTTCTCCGATCACGAACACAGCTTCCCGCAAGAACTTCTGGATAAACTGGTCGAACGCTCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCAGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAGAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACAAAG‑AAACCTTTGGCCTGGATACACTCGACAAACTGTGCGAAATCACCCGTGac                                                                                                                                                                                               >  2:554644/1‑250 (MQ=255)                                   tctcCGATCACGAAACCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAAATCTACGAACTGTGCGCAATCGCCCGTGAag                                                                                                                                                                                              >  2:420927/1‑249 (MQ=255)                                       cgatcaagaccacaacTTCCCGCAAGAGCTGCTGGATAAACTGGTAGAAAGCGCTAATCTGCCGGGTAACCTCGGTAACTGCCACTATTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACAGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCAAGATTGCCTGCCACGACG‑GAACCTTTGGCCTGGCTAAAGTCTACGAACTGTGCGAAATCGCCCGTGACGAGcggc                                                                                                                                                                                        >  2:91607/16‑248 (MQ=255)                                             cGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGTTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGTTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAATAGCTGa                                                                                                                                                                                        <  1:513637/245‑1 (MQ=255)                                              gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGCAACTGCCACTCTTCCGGCACGGTCATTCCGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGCACTGTGCGAAATCGCCAGTGAAGAGCTGACCAAc                                                                                                                                                                                   >  1:300382/1‑249 (MQ=255)                                                aaaCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcggcg                                                                                                                                                                               <  1:159167/251‑1 (MQ=255)                                                    aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTa                                                                                                                                                                           <  1:19774/251‑1 (MQ=255)                                                          ccGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACaacac                                                                                                                                                                     >  2:508074/1‑248 (MQ=255)                                                           cGCAAGAGCTGCTGGATAAACTGGGCAACCGCGCTAATCTGCCGGGTTGCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATc                                                                                                                                                                    <  1:508074/251‑1 (MQ=255)                                                                 agCTGCGGGATAAACTGGTCGAACGCGCTAACCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGATCACATGAAAACCGGCAAGCCGATTTTCGATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCgt                                                                                                                                                              <  1:97798/251‑1 (MQ=255)                                                                  gctgcggtaTAAACTGGGCGAACGCGCTAATCGGACGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACGCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCgtg                                                                                                                                                             <  2:180710/243‑1 (MQ=255)                                                                   ctgctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGAACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTCTACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGACCGTGAAGAGCTGACCAACCGCGGCTATAACATCAGTCgtgt                                                                                                                                                            >  1:56988/1‑251 (MQ=255)                                                                    ggctgaaaaaaCTGGTCGACCGCGCTAATCTGCCGGGTTCCCTCGGTACCTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGtt                                                                                                                                                           <  2:366976/243‑1 (MQ=255)                                                                      ctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCtcccacagcggtcgggtttt                                                                                                                                                         >  2:497009/1‑232 (MQ=255)                                                                        ggATAACCCGGTCCAACGCGCAAACTTGCCGGGTTCCCGCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCCCCGCGGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTGGGCCTGGATAAACTCTACGAACTATGTGAAATCGCCCGTGAAGAGCTGACCAACTGCGGCTACAACATCGGTCGTGTTAt                                                                                                                                                         <  2:630243/249‑1 (MQ=255)                                                                          aTAAACTGGTCGAACGCTCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTGCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCt                                                                                                                                                     <  1:285786/251‑1 (MQ=255)                                                                           tAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTc                                                                                                                                                    >  1:156083/1‑251 (MQ=255)                                                                            aaaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGGTAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTc                                                                                                                                                    <  1:252287/250‑1 (MQ=255)                                                                            aaaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCg                                                                                                                                                   >  1:130787/1‑251 (MQ=255)                                                                            aaaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGAAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTc                                                                                                                                                    <  1:452505/250‑1 (MQ=255)                                                                                     gAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGACGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GACCCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGATCGGGTAATCGCTCGTCCGTTTAt                                                                                                                                          >  1:56211/1‑251 (MQ=255)                                                                                      aaCGCGCTACTCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGCGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATc                                                                                                                                         >  1:173186/1‑251 (MQ=255)                                                                                       aCGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGCTCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGCCGTGTTATCGCCAGTCCGTTTATCg                                                                                                                                        >  2:614451/1‑251 (MQ=255)                                                                                           gcTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCggtga                                                                                                                                    <  1:493700/251‑4 (MQ=255)                                                                                             tAATCTGCCGGGTTAACTCGGAACCTGCCACTCTTCCGGTACGGTCATTCTGGATAAGCTGGGAGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTCTACATCAGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAACtttggcc                                                                                                                                                                                                                                      >  2:82943/1‑148 (MQ=255)                                                                                                           aCCTCGGTAACTGCCACTCTTCCGGTACGGACATTCTGGATCAGCTGGGCGAAGAGCACATGCAAACCGGCAAGCCGATTTTCTATACCTCCGCAGACCCCGTGTTCCAGCTTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACCTCGGTCGCGTTATCGCTCGACCCGTTATCGGCGCCAAATCCGGTTACTTc                                                                                                                    >  2:278604/1‑251 (MQ=255)                                                                                                                             ctTCCGGTACGGTCATTTGGGCTCAGCTGGGCGAAGAGCACATGAAATCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTTGTCCAGAGTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGTCCTACGGCGGCTACAACATCGGTCTTGTTATCGCTCGTCCGTTTATCGGTGACAAGGCCGGTAACTTCCATCGTACCGGTatcc                                                                                                    <  2:360515/249‑4 (MQ=255)                                                                                                                               tCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGt                                                                                                  <  1:137081/249‑1 (MQ=255)                                                                                                                                acGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAAAATCGGTCGTGTTATAGCTCGTCCGTTTTTCGGTGACCAATCCGGGAACTTCCAGCGTACCGGTAAACGTCa                                                                                                >  2:247827/2‑250 (MQ=255)                                                                                                                                                 gATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCCGCGTCCCGGTAACCGTCACGACCTGGCTGTTGAgcc                                                                              >  1:270021/1‑251 (MQ=255)                                                                                                                                                          ggCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTTAGACGCCCGCacc                                                                     >  2:559371/1‑251 (MQ=255)                                                                                                                                                            cGAAGAGCACAGGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACTGCGGCTACAACATCGGTCGTGTTATCGCGCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCaccga                                                                   <  1:369655/251‑1 (MQ=255)                                                                                                                                                               agaacacaTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCaccgaccg                                                                <  2:83932/247‑1 (MQ=255)                                                                                                                                                                   cacaTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgct                                                            <  1:295173/251‑1 (MQ=255)                                                                                                                                                                           aaCCGGCAAGCCGATTTTCTATACCTCCGCTGACTACGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGGTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTACCTGCCAGCGTACCGGGAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAAc                                                    <  2:321881/251‑1 (MQ=255)                                                                                                                                                                             ccGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAATCCGGTAAATTCAAGCGTACCGGTAACCGGCAAGACCTGGATGTTGAGCCGCCAGCACCGACCGTGCGGCAGAAACTg                                                  >  1:580229/1‑251 (MQ=255)                                                                                                                                                                               ggCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGt                                                <  1:599437/251‑1 (MQ=255)                                                                                                                                                                                 cAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGACCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACAAACGGCGGCTAAAACATCGGTCGTGTTATAGCTCGTCCGTTTATCGGTGACAAAACCGGTAACTTCCAGCGTACCGGTAACCGTCACAACCTGGCTGTTGAGCCGCCCGCACCGACCGTGCAGCAGAAACTGGTtg                                              >  1:497988/1‑251 (MQ=255)                                                                                                                                                                                  aaGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGACGCCCGCACCGACCGTGCTGCAGAAAATGGTtg                                              >  1:243332/1‑250 (MQ=255)                                                                                                                                                                                      cGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGTCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGaa                                         <  1:13379/251‑1 (MQ=255)                                                                                                                                                                                               ataCCTCCGCTGACTCCGTGTTACAGCTTGCCTGACACGAAGGAAACTTT‑GGCCTGGATAAACTCTACGACCTCTGCGAAATCGCCCGTGAAGAGCTGACCAACGCCGGCTACAACATCGGTCGTGATATCGGTCGTCCGTTTATCCGTGACAAAGCCGGTAACTTCCAAAGCACCGGTAACACTCACGCCCTGGCTGGTGTGCCGCCAGCACCGACCGTGCTGCAGACACTGGTTGCTGAAAAacagggg                                >  2:267927/1‑247 (MQ=255)                                                                                                                                                                                                 aCCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGGAACCGTCCCGAACTGGCTGTTTGGCCGCCAGCCCCGCCCGTGCTGCAGCAAATGGGTGATGAAAAACACggtc                               >  2:177778/1‑248 (MQ=255)                                                                                                                                                                                                  ccaccGCTGACCCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAg                             <  2:480110/248‑1 (MQ=255)                                                                                                                                                                                                   cTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGa                                                                                 >  1:254278/1‑198 (MQ=255)                                                                                                                                                                                                   cTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGa                                                                                 <  2:254278/198‑1 (MQ=255)                                                                                                                                                                                                     ccGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCACCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACAGTCACGACCTGGCTGTTGAGCCCCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtg                          >  1:61714/1‑251 (MQ=255)                                                                                                                                                                                                     ccGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAgg                            <  2:217192/249‑1 (MQ=255)                                                                                                                                                                                                     ccGCTGACTCCGTGTTCCAAATTGCCTGCCACGAAG‑GACCCTTTGGCCTGGATAAACTCTACGAACTGTGCTAAAGCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCATCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTTCTGCAGAAACTGGTTTATGAAAAACACGGTCAgg                            <  2:441674/249‑1 (MQ=255)                                                                                                                                                                                                       gCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTCCGACCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtggt                        <  1:512251/251‑1 (MQ=255)                                                                                                                                                                                                       gCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCCGTGACAAAGCCGGTAACTTTCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCCGCACCCACCGTGCTGCAGAAACTGGTTGATGACAAACACGGTCAggtggt                        >  2:557375/1‑251 (MQ=255)                                                                                                                                                                                                       gCTGACGCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGTATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACAAACGGCGGCTACAACATCGGTCGTGTTAACGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCATCACCGTCCGTGCTGCATAAACTGGTTGATGAAAAACACGGTCAggtggt                        <  1:205520/251‑1 (MQ=255)                                                                                                                                                                                                          gACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTc                     <  1:139523/251‑1 (MQ=255)                                                                                                                                                                                                           aCTCCGGGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACACCATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCt                    <  2:130787/251‑1 (MQ=255)                                                                                                                                                                                                                gtgtTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAAAATCGGTCGTGTTATCGATCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGTCgg               >  1:560342/1‑251 (MQ=255)                                                                                                                                                                                                                  gtTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGt              <  1:402771/250‑1 (MQ=255)                                                                                                                                                                                                                  gtTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTAGCGCGCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTCCCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGt              <  2:128811/250‑1 (MQ=255)                                                                                                                                                                                                                    tCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTAAAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGCTGAAAAACACGGTAAGGTGGTTTCTGCCGGTaa            >  1:360831/1‑250 (MQ=255)                                                                                                                                                                                                                       aGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAAtt        >  1:390006/1‑251 (MQ=255)                                                                                                                                                                                                                             cctgccACGAAG‑GAACCTTGGGCCTGGATAAACTCTCCGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCATCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGAc  <  2:399168/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           GCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAG‑AAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAC  >  NC_000913/4619918‑4620388

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