Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,150	A→G	E183G (GAA→GGA)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,150	0	A	G	96.0%	220.7 / ‑4.2	74	E183G (GAA→GGA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (33/38);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAG‑AAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAC  >  NC_000913/4619910‑4620388                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               gggAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAAcctttg                                                                                                                                                                                                                                        <  1:63729/249‑6 (MQ=255)   gaaatcgacGGTGTCCAGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCAGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGACGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGCTTGCCTGCCACGAAG‑GAAcctttggcct                                                                                                                                                                                                                                    >  2:542928/9‑242 (MQ=255)     aatcgCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAAcctttggcctgg                                                                                                                                                                                                                                  >  1:544973/5‑240 (MQ=255)         gCCGTTGTCCCGGTTCTGTTGGAGTGGGGTGATTCCCCCGACCATGAAAACAGCTCCCCGCAAGTGCTGCTGGACAAACTGTCCGAACGCGCTAATCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAAccgttggcctggataa                                                                                                                                                                                                                              <  2:134289/251‑16 (MQ=255)         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                                                                                 <  2:254278/198‑1 (MQ=255)                                                                                                                                                                                                           cTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGa                                                                                 >  1:254278/1‑198 (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)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               GGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAG‑AAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAC  >  NC_000913/4619910‑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: ‑