Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,619,906	C→A	93.0%	H102N (CAC→AAC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,906	0	C	A	93.0%	299.7 / 6.0	97	H102N (CAC→AAC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base A (48/42);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGT  >  NC_000913/4619760‑4620054                                                                                                                                                   |                                                                                                                                                     tGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAG‑TATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACt                                                                                                                                                   >  6:38462/1‑149 (MQ=255)  gCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACt                                                                                                                                                   >  1:148777/1‑149 (MQ=255)    cAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTgg                                                                                                                                                 >  3:84260/1‑149 (MQ=255)    cAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTgg                                                                                                                                                 >  4:53065/1‑149 (MQ=255)    cAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTgg                                                                                                                                                 <  1:99938/149‑1 (MQ=255)      aaTCTGACCCGTCTGGGGCTGGCGAAAGCACACGACGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCAGGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTggg                                                                                                                                                >  6:25338/1‑148 (MQ=255)           gACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGCTACCCCGTCTGGTAACTGGGAAATTg                                                                                                                                          >  5:230049/1‑149 (MQ=255)            aCCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTg                                                                                                                                          <  3:126952/148‑1 (MQ=255)              ccGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCg                                                                                                                                       <  1:203123/149‑1 (MQ=255)                    ggtgCTGGCGAAAGCACACGAAGGTCCTACCGGTGTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCAGGGGCGCACGAAATTTCATCCGGTAAAGATTCCCCGTCTGGTCACTGGGAAATTGCCGGTGTcc                                                                                                                                 <  2:328732/146‑1 (MQ=255)                     gggCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTcc                                                                                                                                 <  5:233206/148‑1 (MQ=255)                     gggCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGAGGGAATGGACGGCAACGCTGAAG‑TATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCg                                                                                                                               <  5:38462/149‑1 (MQ=255)                        cTGGCGAATGCACACTAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCGGAAGTTATCGGCGCGTACGCAGGGGCGCACGATATGTCATACGGGAAAGATACACCTTCTGGGAACTGGGAAATTGACGGTGTCCCGGt                                                                                                                             >  1:39846/1‑149 (MQ=255)                        cTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGt                                                                                                                             >  3:161216/1‑149 (MQ=255)                           gCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCt                                                                                                                          >  4:2611/1‑149 (MQ=255)                           gCGAAAGCACACGAAGGTTCTACCGGTTTCATTCAGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCt                                                                                                                          >  5:326260/1‑149 (MQ=255)                             gAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGt                                                                                                                        <  3:132121/149‑1 (MQ=255)                             gAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGt                                                                                                                        >  5:355977/1‑149 (MQ=255)                             gAAAGCACACCAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGAAACCCTGAAGTTATCGCCGCGTACGCATGGGCGCACGAAATGTCATCCTGAAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGt                                                                                                                        >  6:70621/1‑149 (MQ=255)                               aaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGttt                                                                                                                      <  2:148777/149‑1 (MQ=255)                               aaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGttt                                                                                                                      <  2:56888/149‑1 (MQ=255)                               aaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGttt                                                                                                                      <  3:210065/149‑1 (MQ=255)                               aaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGt                                                                                                                        <  7:135578/147‑1 (MQ=255)                               aaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGt                                                                                                                        >  8:135578/1‑147 (MQ=255)                                 gCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGttt                                                                                                                      >  7:337117/1‑147 (MQ=255)                                 gCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGttt                                                                                                                      <  8:337117/147‑1 (MQ=255)                                 gCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                                    <  3:92981/149‑1 (MQ=255)                                  cacaCGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGGCCCGGTTCTGTTTGAg                                                                                                                   <  3:32556/149‑1 (MQ=255)                                     acGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCAGCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                                >  6:131120/1‑149 (MQ=255)                                       gAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCATGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTggg                                                                                                               >  7:359291/1‑148 (MQ=255)                                       gAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTggg                                                                                                               >  7:29153/1‑148 (MQ=255)                                       gAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTggg                                                                                                               >  6:377891/1‑148 (MQ=255)                                           gTTCTACCGGTTTCATTCCGGCGGGAATTGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGTGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGatat                                                                                                          <  5:281959/149‑1 (MQ=255)                                               tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTtct                                                                                                      >  8:324224/1‑149 (MQ=255)                                               tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTtct                                                                                                      >  5:204764/1‑149 (MQ=255)                                               tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGACGGTGGCCCGGTTCTGTTTGAGTGGGGATATTtct                                                                                                      >  5:199809/1‑149 (MQ=255)                                                aCCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTtctc                                                                                                     >  1:338240/1‑149 (MQ=255)                                                 ccGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGTGGATATTtct                                                                                                      <  7:82874/147‑1 (MQ=255)                                                 ccGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTtctc                                                                                                     <  7:136491/148‑1 (MQ=255)                                                 ccGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTcc                                                                                                    >  5:22522/1‑149 (MQ=255)                                                  cGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAGAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTcc                                                                                                    >  1:248484/1‑148 (MQ=255)                                                        cATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAc                                                                                              <  6:355977/148‑1 (MQ=255)                                                         aTTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACg                                                                                             >  4:177559/1‑148 (MQ=255)                                                           tCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGAAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGaa                                                                                           >  8:269773/1‑148 (MQ=255)                                                             cggcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGaaaa                                                                                         <  5:18787/148‑1 (MQ=255)                                                             cggcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGaaaa                                                                                         <  5:377891/148‑1 (MQ=255)                                                             cggcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAAc                                                                                        <  2:88665/149‑1 (MQ=255)                                                               gcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAg                                                                                      <  4:84260/149‑1 (MQ=255)                                                                 gggAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCt                                                                                    >  4:325545/1‑149 (MQ=255)                                                                  ggAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCtt                                                                                   >  1:24915/1‑149 (MQ=255)                                                                   gAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTc                                                                                  <  4:33076/149‑1 (MQ=255)                                                                        gACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGc                                                                              >  8:83374/1‑148 (MQ=255)                                                                        gACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGc                                                                              >  1:56296/1‑148 (MQ=255)                                                                                 gCTGAAGTTATCGGCGCGTATGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgct                                                                    <  2:46658/149‑1 (MQ=255)                                                                                 gCTGAAGTTATCGGCGCGTACGCATGGGCGCATGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgct                                                                    <  8:359291/149‑1 (MQ=255)                                                                                   tGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgc                                                                     >  4:79882/1‑146 (MQ=255)                                                                                   tGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgc                                                                     <  3:79882/146‑1 (MQ=255)                                                                                         tATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACACCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAAc                                                            >  4:208411/1‑149 (MQ=255)                                                                                          aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCTGTGTCCCGGTTCTGTTTGAGTGGGTATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAAc                                                            <  7:142006/148‑1 (MQ=255)                                                                                          aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAAc                                                            <  8:29153/148‑1 (MQ=255)                                                                                          aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAAc                                                            <  3:325545/148‑1 (MQ=255)                                                                                          aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAAc                                                            <  2:24915/148‑1 (MQ=255)                                                                                          aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAAc                                                            <  6:332297/148‑1 (MQ=255)                                                                                           tCGGCGCGTTCGCATGGGCGCACGAAATGGGATCCGGTAAAGATTCCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTg                                                          <  5:130081/149‑1 (MQ=255)                                                                                           tCGGCGCGTCCGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAAC‑AGGAAATTGCCGGTGTCCCGGTTCTGGTTGAGTGGGGCTATTTCACCGATCACGAAACCAGCTTCCCGCAAGAGCTGCTGGATAACCTgg                                                         >  8:395500/1‑149 (MQ=255)                                                                                              gcgcGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGTGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTc                                                       >  2:143107/1‑149 (MQ=255)                                                                                              gcgcGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTc                                                       >  3:171647/1‑149 (MQ=255)                                                                                              gcgcGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTc                                                       <  3:2611/149‑1 (MQ=255)                                                                                                  gTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGGGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAc                                                   <  6:22522/149‑1 (MQ=255)                                                                                                   tACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcg                                                  >  4:352095/1‑149 (MQ=255)                                                                                                     cGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgcg                                                >  7:163759/1‑149 (MQ=255)                                                                                                       cATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACgcgc                                               >  3:287242/1‑148 (MQ=255)                                                                                                        aTGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATCTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCt                                              >  8:93349/1‑148 (MQ=255)                                                                                                           ggCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAAt                                           <  6:104804/148‑1 (MQ=255)                                                                                                           ggCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATc                                          <  1:315742/149‑1 (MQ=255)                                                                                                                aCGAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCg                                     <  5:25338/149‑1 (MQ=255)                                                                                                                  gAAATGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCggg                                   >  5:329275/1‑149 (MQ=255)                                                                                                                   aaaTGTCATCCGGTAAAGAACCCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                                  <  8:387197/149‑1 (MQ=255)                                                                                                                    aaTGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCTGTGTCCCGGTTCTGTTTGACTCGGCATATTTCTCCGATAACGAAAACAGCTTCCCGCAAGCGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                                  >  5:143997/1‑148 (MQ=255)                                                                                                                    aaTGTCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGtt                                 <  3:53065/149‑1 (MQ=255)                                                                                                                        tCATCCTGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTAcc                              <  8:347581/148‑1 (MQ=255)                                                                                                                        tCATCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCt                             >  8:201637/1‑149 (MQ=255)                                                                                                                          aTCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTTGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCg                           >  8:105407/1‑149 (MQ=255)                                                                                                                          aTCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTc                            <  7:83374/148‑1 (MQ=255)                                                                                                                          aTCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCg                           >  6:193128/1‑149 (MQ=255)                                                                                                                          aTCCGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCg                           >  6:186831/1‑149 (MQ=255)                                                                                                                            ccGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCgg                          >  8:88654/1‑148 (MQ=255)                                                                                                                             cGGTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGt                         <  7:324224/148‑1 (MQ=255)                                                                                                                              ggTAAAGATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGTTTCTGTTTGACTGGGGATATTTATCAGATCACGAAAACAGCTTCCCGCAAGAGCTGCTTGATAAACTGGTCGAACGCGCTAATCTGCCGGGATACCTCGGTaa                       >  6:300335/1‑149 (MQ=255)                                                                                                                                    gATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGcc                  >  1:169376/1‑148 (MQ=255)                                                                                                                                    gATACCCCGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa                 <  5:70621/149‑1 (MQ=255)                                                                                                                                          ccGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTc           >  3:75052/1‑149 (MQ=255)                                                                                                                                          ccGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTc           <  8:189258/149‑1 (MQ=255)                                                                                                                                           cGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACt                     <  1:133483/138‑1 (MQ=255)                                                                                                                                           cGTCTGGTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACt                     >  2:133483/1‑138 (MQ=255)                                                                                                                                                ggTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACAGGTc                                                       <  8:152916/99‑1 (MQ=255)                                                                                                                                                ggTAACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACAGGTc                                                       >  7:152916/1‑99 (MQ=255)                                                                                                                                                   cACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGt  >  8:229251/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGT  >  NC_000913/4619760‑4620054

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