Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,312	T→C	G28G (GGT→GGC)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,312	0	T	C	95.9%	214.2 / ‑1.6	74	G28G (GGT→GGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (34/36);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATT  >  NC_000913/4618067‑4618555                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ggccACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTcg                                                                                                                                                                                                                                                <  1:116612/248‑1 (MQ=255)    cACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTcgcg                                                                                                                                                                                                                                              >  1:309416/1‑250 (MQ=255)     aCGGTGACGATAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACACGAAGGCCGGCTAAGGCGTAGACGCATCCGACAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGATACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACATGATGGACATGCACTGAGAGATGAAGACATTAGATTCTTTATCCACGGCATAcgcgcc                                                                                                                                                                                                                                            >  2:681428/1‑249 (MQ=255)               tAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTCGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGGGATGGTCATGCACTGAGCGCTGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTAtct                                                                                                                                                                                                                                  >  2:546175/1‑251 (MQ=255)                cagagCGCCAGCAGCTACTAAGAAAGATGCTTTACTCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCAGCCTCATACTCTTTTCCTCGGGCGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAAAGTGAGGGTCATGCACTGAGCGTTGAAGAAATTCGTTTCTTTATAAACGGCACTCGAGACAACACtctatc                                                                                                                                                                                                                                 >  2:532936/2‑246 (MQ=255)                aaGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTAtctc                                                                                                                                                                                                                                 <  1:155270/251‑1 (MQ=255)                           catcgtataagtaagATGCTTTACGCCTGATTCGCTGCGCTTATCAGG‑CCTACTGGACGTATTTACCCTTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGTCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGACGAAGAAATTCGTTTCTTTAGCAACGGCATTCGCGACAACACTATCGCCGAAGGGCAGa                                                                                                                                                                                                                      <  2:643777/244‑1 (MQ=255)                                  taagtaagATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTgccgc                                                                                                                                                                                                               >  1:184547/1‑251 (MQ=255)                                   aagtaagATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTgccg                                                                                                                                                                                                                <  1:364265/249‑1 (MQ=255)                                      taagATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CATACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTc                                                                                                                                                                                                           >  1:178108/1‑251 (MQ=255)                                       aagATGCTTTACGCCTGATTCGTGGGGCTTATCAGG‑CCTTCGGGCCGTATCTTCCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGTCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTc                                                                                                                                                                                                           <  2:400001/250‑1 (MQ=255)                                           tcctttaaggcTGATGCGCTGCGCTTAT‑AGGACCTAGGAAACGTATCTACCCGTCGGCCGGATAAGGCCTAGACGCATCAGGCAGAGGCCGCATCATACTCGTTTCATCGGTAGGGTCCCGTGGTTGTCGGACTAGAACTTAGTCGTAATAAACGTGATGGTCAAGCGCTGAGCGCGGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGGCAACAATGTCTCCGAAGGGCCGGTTGCAGCACTcgcaat                                                                                                                                                                                                      >  2:713/11‑248 (MQ=255)                                                  cGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCAtt                                                                                                                                                                                               >  1:667498/1‑251 (MQ=255)                                                   gCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGCTGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGCCATTCGCGACAACACTATCTCCGCAGGGCAGATTGCCGCCCTCGCGATGACCAttt                                                                                                                                                                                              >  1:314432/1‑251 (MQ=255)                                                                   tatgctgg‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTATACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAAGTATTCGTAAAAAACGTGATGGTCATGCTCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACTATTTTCTTCCACGATATGa                                                                                                                                                                               <  2:416489/244‑1 (MQ=255)                                                                       cAGG‑CCTACGAGACGTATCTCCCAGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTCAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCACCGGAATTAGCGACAACACTATCTCCGAAGGGAAGATTGCCGCACTCGCGATGACCATTTTCTTCCCCGATATGACAATg                                                                                                                                                                          >  2:680775/1‑251 (MQ=255)                                                                        aGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGc                                                                                                                                                                         >  1:164723/1‑251 (MQ=255)                                                                          g‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCt                                                                                                                                                                       <  1:6262/251‑1 (MQ=255)                                                                                cGAGACGTATCTACCCGTAGGCCGGATAAGACGTCGTCCCATCCGGCCCACGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAACTTCGATTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGAGCAGATTGCCGCCCTCGCGTGGACCATTTTCTTCCACGATCTGACAATGCCTGGGCg                                                                                                                                                                  >  1:602510/1‑251 (MQ=255)                                                                                 gagaCGTAACTACCCGGTGGCCGGATAAGGCGTAGACGCAGCCGGCAAAAGCCTCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCGCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCgt                                                                                                                                                                 <  1:563166/251‑1 (MQ=255)                                                                                      gTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATCCTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCTCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGtctc                                                                                                                                                            <  1:365161/251‑1 (MQ=255)                                                                                            aCCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGAc                                                                                                                                                      <  1:24585/251‑1 (MQ=255)                                                                                              ccGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGAcc                                                                                                                                                     >  1:554177/1‑250 (MQ=255)                                                                                                gTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATg                                                                                                                                                  >  1:192080/1‑251 (MQ=255)                                                                                                       ggATTAGGCGTAGACTCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGc                                                                                                                                           <  2:597324/251‑1 (MQ=255)                                                                                                         aTAAGGCGTAGACGCATCCGGCAAAAGCAGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGAAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCg                                                                                                                                          >  1:483595/1‑250 (MQ=255)                                                                                                                   gACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGa                                                                                                                                >  1:11353/1‑250 (MQ=255)                                                                                                                              cAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTc                                                                                                                                                                                                           <  1:402146/164‑1 (MQ=255)                                                                                                                              cAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGTTTCAGGACCCGTTCTCGAc                                                                                                                    >  2:443023/1‑251 (MQ=255)                                                                                                                                   gccgccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGaa                                                                                                               <  1:620202/251‑1 (MQ=255)                                                                                                                                     cgccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCACCGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCTCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGaaa                                                                                                              >  2:547070/1‑250 (MQ=255)                                                                                                                                     cgccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGaaa                                                                                                              >  1:600509/1‑250 (MQ=255)                                                                                                                                       ccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGc                                                                                                                                                <  1:580515/214‑1 (MQ=255)                                                                                                                                       ccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGc                                                                                                                                                >  2:580515/1‑214 (MQ=255)                                                                                                                                             aCTCTTTTCATCGTGAGGTTTCCTTGTTTCTCGCACAAGAAATTAGTCGTAAAAAACGGGATTGTCAGGCACTGAGAGATGAAGAAATTCGTTTCTTTATCAACGGCAGTCGAGACAATACTATCTCCGAAGGGCAGGTTGCCGCCCTCGCGATGAACATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGATCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAGGCCTGCa                                                                                                      <  2:597892/250‑1 (MQ=255)                                                                                                                                                     ccTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATg                                                                                              <  1:658722/250‑1 (MQ=255)                                                                                                                                                        cGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATTACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGccc                                                                                          <  1:681428/251‑1 (MQ=255)                                                                                                                                                         gggAGGTTTCCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAAAGTGATGGTCATGCACTGAGCGATGAAGAAATTCTTTTCTTTATCAACGGCATGCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGGGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTGCTCGACTGGAAAAGCCTGCATCTGAATGGCCCg                                                                                         <  2:45776/251‑1 (MQ=255)                                                                                                                                                          ggAGGTTACCTTGTTTCTCGCACAAGAAATTAGTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGTAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAATGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAATCTTTCTCGACTGGAAAAGCCTGCATCTGAATGGcc                                                                                           <  2:658054/248‑1 (MQ=255)                                                                                                                                                               ttACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTcgcg                                                                                                                                                                                                        >  2:375624/1‑134 (MQ=255)                                                                                                                                                               ttACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTcgcg                                                                                                                                                                                                        <  1:375624/134‑1 (MQ=255)                                                                                                                                                                 tccTTGTTGCCCGCACAAGAAAATATTCGTAAAAACCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCGCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGAttgttg                                                                                  <  2:361568/249‑1 (MQ=255)                                                                                                                                                                  ccTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGACGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACAATATCTCCGAAGGGCAGATTGCCGACCTCGCGCTGACCATTTTCTTCCACGATCTGACAATGCCTGAGCGTGTCGCGATGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGa                                                                                 >  2:673803/1‑250 (MQ=255)                                                                                                                                                                      gTTTCGCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTTATaaa                                                                             <  2:600509/250‑1 (MQ=255)                                                                                                                                                                        ttCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAaca                                                                           <  1:599107/250‑1 (MQ=255)                                                                                                                                                                            cGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAacatt                                                                         <  1:288046/248‑3 (MQ=255)                                                                                                                                                                            cGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAacatt                                                                         <  1:649906/248‑3 (MQ=255)                                                                                                                                                                            cGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTcc                                                                       <  2:309416/250‑1 (MQ=255)                                                                                                                                                                              cacaAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTtctc                                                                                                                       <  1:395056/200‑1 (MQ=255)                                                                                                                                                                              cacaAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTtctc                                                                                                                       >  2:395056/1‑200 (MQ=255)                                                                                                                                                                               acaAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCActg                                                                   >  1:346321/1‑249 (MQ=255)                                                                                                                                                                               acaAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATAGCCCGATTGGTGATAAAAATTCCActg                                                                   >  2:465840/1‑249 (MQ=255)                                                                                                                                                                                   gAAATTATGCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTAGCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGCCCATGGCGATGCGAGATTCAGGAACCGTtc                                                                                                                         <  1:70776/193‑1 (MQ=255)                                                                                                                                                                                     aaTTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATGGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCg                                                             <  2:361341/251‑1 (MQ=255)                                                                                                                                                                                      attattCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGt                                                            <  1:581628/251‑1 (MQ=255)                                                                                                                                                                                         attCGTAACAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGTAGAGCAGATTGCCGCCCTCGCGATGCCCATTTTCTTCCCCGATATGACAATGCCTGAGCGTGTCTCGCTGACAATGGCGATGCGAGATTCAGGCACCGTTCTCGCCTGGCAAAGCCTGAATCTAGATGACCCGCTTGGTTATAAacattcacatggtggcgtagg                                                         >  2:518031/1‑233 (MQ=255)                                                                                                                                                                                         attCGTAAAAAACGGGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCg                                                          <  2:164723/250‑1 (MQ=255)                                                                                                                                                                                          ttCGTAACAAAAGTGAAGGTCATGCACTACGCGATGAAGAAATTCGTTTATTTATCAACGGCATTCGAGCCAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCATGTCGCGCTGACCATGGCGATACGAGATTCAGGAACGGTTCTAGACTTGAAAAGCATGCATGTGAATGGCACGCATGGTGATAAATATTCCACTGGTGGCGTCGGc                                                        >  2:72791/1‑251 (MQ=255)                                                                                                                                                                                          ttCGTAAAAAACGTGATGGTCATGCCCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGAGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGGATGGCCCGCTTGTTGATAAACATTCCACTGGTGGCGTCGGc                                                        >  1:492518/1‑251 (MQ=255)                                                                                                                                                                                          ttCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGc                                                        >  1:149719/1‑251 (MQ=255)                                                                                                                                                                                           tCGTAAAAAACGGGATGGGCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGc                                                        <  2:667498/250‑1 (MQ=255)                                                                                                                                                                                                aaaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGAtgtg                                                  >  1:173055/1‑251 (MQ=255)                                                                                                                                                                                                  aaaCGTTATGGTCATGCACTGAGCGATGAAGAAATTCGGTTCTTTATCAACGGCATTCGCGACAACTCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGAGGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGaa                                                <  2:466643/251‑2 (MQ=255)                                                                                                                                                                                                      gTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCAATGGTGGCGTCGGCGATGTGACTTCg                                            >  2:567911/1‑251 (MQ=255)                                                                                                                                                                                                      gTGATGATCATGCACTGAGCGATGAAGCAATTCGTTTCTTTATCAACGGCATTCGAGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCCTTTTCTTCCACGATATGCCAATGCCTGAGCGTGTCTCGCTGCCCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGCTTGTTGATAAACATTCCACTGGGGGCGTCGTCGATGTGACTTCg                                            >  2:98009/1‑251 (MQ=255)                                                                                                                                                                                                                      gagagATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACCCTAGCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGcc                              <  2:192080/245‑1 (MQ=255)                                                                                                                                                                                                                             gaagaaATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCAGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGCTGTGCCTTCGCTGCTGTTGGGGCCGATGGTCGc                     >  2:128305/1‑251 (MQ=255)                                                                                                                                                                                                                               agaaATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGc                     >  1:190156/1‑249 (MQ=255)                                                                                                                                                                                                                                     tCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGCTGGTCGCAGCCTgcg             >  2:336345/1‑251 (MQ=255)                                                                                                                                                                                                                                        tttCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAAcc                                                                                                                             >  1:357473/1‑136 (MQ=255)                                                                                                                                                                                                                                        tttCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAAcc                                                                                                                             <  2:357473/136‑1 (MQ=255)                                                                                                                                                                                                                                                 cAACGGCATACGCGACAACACTATCTCCGAAGGGCCGAGTGCCGCCCTCGCGATGACCATTTTCTACCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATGCAGGAACCGTGCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGGCGGCGATGTGACTGCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATAtt  <  2:158222/250‑1 (MQ=255)                                                                                                                                                                                                                                                     ggTATTCGCGACAACCCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATaaccaa                                                                          >  2:260729/1‑170 (MQ=255)                                                                                                                                                                                                                                                     ggTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAacac                                                                          <  1:260729/174‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    GGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATT  >  NC_000913/4618067‑4618555

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