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	88.9%	87.4 / 3.3	37	G28G (GGT→GGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (15/17);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

CGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTG  >  NC_000913/4618078‑4618542                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        cGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTa                                                                                                                                                                                                                         <  1:70211/251‑1 (MQ=255)   gTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAGGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTAtc                                                                                                                                                                                                                       >  1:116576/1‑251 (MQ=255)    tAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCACCGGCATTCGCGCCAAAACTAtct                                                                                                                                                                                                                      >  1:50402/1‑251 (MQ=255)                           taagATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTCCCTTGTTTCTCGCACAAGAAATTATTC‑GTACAAAACGTGATGGTCATGCGATGAGCGATGAAGAAATTCGTTTCTTTATCAACAGAATTCGCGCCAACACTCACTCCGATGGGCAGATTGCCGCCCTc                                                                                                                                                                                               >  2:197277/1‑251 (MQ=255)                            aagATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTcg                                                                                                                                                                                              <  1:78225/251‑1 (MQ=255)                                        gCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCAttt                                                                                                                                                                                  >  1:17340/1‑251 (MQ=255)                                           gtgtgctctgcgctTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGGAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCTTCGCGATGACCAGTttct                                                                                                                                                                               <  2:231664/244‑1 (MQ=255)                                                      gctTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTAtt                                                                                                                                                                                                                                      <  1:6593/185‑1 (MQ=255)                                                      gctTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTAtt                                                                                                                                                                                                                                      >  2:6593/1‑185 (MQ=255)                                                           ccaatcataagagtcGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACCCTATCTCCGAAGGGCAGATTGCCGCACTCGCGATGCCCATTTTCTTCCACGATATGACAAt                                                                                                                                                               <  2:161852/237‑1 (MQ=255)                                                                                           ggATAAGGCGTAGACGCATTCGGCAAAAGCCGCCTCATACTCTTGTCCTCGGGAGGTTACCGTGATGCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCCCTGAGCGATGAAGAAATGCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGTCCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGc                                                                                                                               <  2:108886/251‑1 (MQ=255)                                                                                                               cGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGATCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTtctc                                                                                                          <  1:235992/251‑1 (MQ=255)                                                                                                                         cgtcTCATACTCTTGTCCTCGGGAGGTTGCCTTGTTTCTCTAGCAAGAAATTTTTCTGTAAAAT‑CGTGATGGTCATTCACTGAGCGATGAAGAAATTCGATTCTTTATCTACGGCATTCGCGACAATACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCGGAGCGTGTCTCGCTGACCATGGCGATTCGAGA‑TTCAGGATCCGTTCTCTCCTGGaaaa                                                                                                <  2:257366/248‑1 (MQ=255)                                                                                                                              cATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCt                                                                                            <  1:282715/250‑1 (MQ=255)                                                                                                                              cATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTg                                                                                           <  1:104507/251‑1 (MQ=255)                                                                                                                               ctACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAATACCCGGGATTACAATCACCTGAGCGATGAAGAAATTCATGTCTTTAGAAACGGCATTCGCGACAACACTGTCTCCTAAGGGCAGGTTGACGCCATAGCGAGTACCATTTTCTTCCACGATAAGACCATGCCCGCGCGAGTCTCGCTAACATTGGCGATAC‑AGACTTCAGGAACCGGTCTCGACTGGAAAAGccgg                                                                                           >  2:88629/2‑248 (MQ=255)                                                                                                                                tACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGc                                                                                          >  1:189315/1‑250 (MQ=255)                                                                                                                                         ccTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGCCCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATgg                                                                                >  1:33187/1‑251 (MQ=255)                                                                                                                                         ccTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATgg                                                                                >  1:181613/1‑251 (MQ=255)                                                                                                                                            cGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTCTCTCCGACGGGAAGCTTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAAGGCCTGCACTTGTCTCGCTGACCTTGGCGATGCGAGC‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCAGCTGCATGGccc                                                                             >  2:99295/1‑251 (MQ=255)                                                                                                                                                    tACCGTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGAttgtt                                                                      <  1:77393/250‑1 (MQ=255)                                                                                                                                                     aCCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCAGGCACTGAGCGATGAAGAAATTCGTTTCTTTTTCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGAttgttg                                                                     <  2:134946/250‑1 (MQ=255)                                                                                                                                                                cGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTcc                                                          <  1:374673/250‑1 (MQ=255)                                                                                                                                                                   acaAGAAATTATTC‑GTAAAAAACGTGATGGGCATGCGCTGAGCGATGAGGAGATTCGTTTCTTTATCAACGGTATCCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACAATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACAATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAACCCCTGAATCTGAATGGCCCGATTGTTGATAAACACTCCACCg                                                      >  2:152873/1‑251 (MQ=255)                                                                                                                                                                            tattC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTc                                              <  1:92958/250‑1 (MQ=255)                                                                                                                                                                            tattC‑GTAAAAAACGTGAGGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTAGCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTc                                              <  2:181613/250‑1 (MQ=255)                                                                                                                                                                                   tAAAAAACGTGATGGTCATGCACTGAGAGATGAAGCAATTCGTTTCTTTATCAAAGGCATTCGCGACAACACTATCTCCAAAGGGCAAATTGACGCCCTCGCGATGACCATTTTCTTCCACGATATGAACATGCCTGAGCGTGTCTCGCTGACCATGACGATGCGAGA‑TTCAGGACCCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGGATGTTGATAAACATTCCACTGGTGTCGTAGGCgctg                                       >  2:384639/1‑248 (MQ=255)                                                                                                                                                                                      aaaaCGTGATGGTCTTTCGCTGAGCTATTAAGAAATGCTTTTCTTAATAAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCAGTTTCTTCCACGATATGACAATGCCTGAGCGTGCCGCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGa                                    <  2:110091/251‑1 (MQ=255)                                                                                                                                                                                      aaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAA‑GCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGAc                                   >  1:63836/1‑251 (MQ=255)                                                                                                                                                                                       aaaCGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGAc                                   <  1:135826/251‑1 (MQ=255)                                                                                                                                                                                       aaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGc                                                 >  2:178028/1‑237 (MQ=255)                                                                                                                                                                                       aaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGc                                                 <  1:178028/237‑1 (MQ=255)                                                                                                                                                                                         aCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACtt                                 >  1:197205/1‑251 (MQ=255)                                                                                                                                                                                                  tCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGt                        >  1:348669/1‑251 (MQ=255)                                                                                                                                                                                                        acTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGc                  <  1:269412/250‑1 (MQ=255)                                                                                                                                                                                                             gCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATg             >  1:368499/1‑251 (MQ=255)                                                                                                                                                                                                                        aTTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTg  >  1:381655/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        CGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTG  >  NC_000913/4618078‑4618542

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