Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	114,203	C→T	6.0%	R254C (CGC→TGC)	guaC →	GMP reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	114,203	0	C	T	6.0%	219.3 / 2.9	67	R254C (CGC→TGC)	guaC	GMP reductase
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base T (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

TGCCGATGCTGCGCACGGTCTGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGA  >  NC_000913/114058‑114348                                                                                                                                                  |                                                                                                                                                  tgacGATGCTGCGCACGGTCTGGGCGGAATGATCGTCAGCGATGGTGGCAGCACCACGCCGGGCGATGTGGAGACAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCTGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCg                                                                                                                                                  >  1:189250/4‑147 (MQ=255)      aTGCTGCGCACGGTCTGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGt                                                                                                                                            <  4:84100/148‑1 (MQ=255)         ctgcGCACGGTCTGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCGGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTg                                                                                                                                          <  1:263622/147‑1 (MQ=255)         ctgcGCACGGTCTGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTg                                                                                                                                          <  3:195645/147‑1 (MQ=255)              cACGGTCTGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTgaggag                                                                                                                                     >  3:343640/1‑147 (MQ=255)              cACGGTCTGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTgaggag                                                                                                                                     >  3:133534/1‑147 (MQ=255)                 ggTCTGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACg                                                                                                                                 <  6:331011/148‑1 (MQ=255)                 ggTCTGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACg                                                                                                                                 >  1:323694/1‑148 (MQ=255)                     tGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCga                                                                                                                             <  4:343640/148‑1 (MQ=255)                       ggCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCgag                                                                                                                            >  2:399772/1‑147 (MQ=255)                           gAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAAtt                                                                                                                       <  3:130148/148‑1 (MQ=255)                            aaTGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCTAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCTGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAAttt                                                                                                                      <  3:143675/148‑1 (MQ=255)                               gATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATg                                                                                                                   >  1:269942/1‑148 (MQ=255)                               gATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATg                                                                                                                   >  1:311344/1‑148 (MQ=255)                                           tggtggCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTg                                                                                                                                          <  3:220644/113‑1 (MQ=255)                                           tggtggCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGc                                                                                                       <  1:399772/148‑1 (MQ=255)                                            ggtggCTGAACCACGCCGGGCGATGTCGCGAAAGCCTATGGCGGCGGTGCCGATTACGTCATGCTTGGCGTCAAGCTGGCGGGCCACGAAGAGAGCGGCGTTCGCATCGTTga                                                                                                                                         >  4:220644/1‑113 (MQ=255)                                                    caccacTCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAATAGAGCGGCGGTCGCATCTTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTc                                                                                               <  4:225788/147‑1 (MQ=255)                                                    caccacGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTcc                                                                                              >  2:328387/1‑148 (MQ=255)                                                      ccacGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGa                                                                                                   >  7:308924/1‑141 (MQ=255)                                                      ccacGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGa                                                                                                   <  8:308924/141‑1 (MQ=255)                                                       cacGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGgaga                                                                                                                                    <  1:102984/107‑1 (MQ=255)                                                       cacGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGgaga                                                                                                                                    >  2:102984/1‑107 (MQ=255)                                                              ggCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGa                                                                                    <  2:269942/148‑1 (MQ=255)                                                               gCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCgtg                                                                                   >  3:23166/1‑146 (MQ=255)                                                                  aTGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGaa                                                                                <  1:328387/148‑1 (MQ=255)                                                                                  tGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGCAACGTCACGTTggcggc                                                                >  3:54238/1‑148 (MQ=255)                                                                                            gCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGc                                                                                                       <  4:55170/99‑1 (MQ=255)                                                                                            gCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGc                                                                                                       >  3:55170/1‑99 (MQ=255)                                                                                             ccGAGTTCGTCATGCTTTGCGGGATGCTGGCGTGCCACGAAGAGAGCGGCGGTCGCATCTTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTAGGGCATGAGCTCCGAGTCTGCGATGAAACGTGACGTTGGCGGCGTTGCGGAata                                                     <  2:189250/148‑1 (MQ=255)                                                                                                 tttCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATcgc                                                 >  7:122848/1‑148 (MQ=255)                                                                                                        aTGCTTGGCGGCATGCTGGCGGGCCCCGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGCGCTCCGAGTCTGCGATGAACCGTCACGGTGGCGGCGTTGCGGTATATCGCgcag                                             >  6:41640/1‑145 (MQ=255)                                                                                                          gCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGgaga                                                                                                                                    <  2:129926/56‑1 (MQ=255)                                                                                                          gCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGgaga                                                                                                                                    >  1:129926/1‑56 (MQ=255)                                                                                                          gCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGcagcag                                          >  8:22351/1‑146 (MQ=255)                                                                                                           cTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa                                        <  2:323694/147‑1 (MQ=255)                                                                                                              ggcggcATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGt                                     >  2:275324/1‑147 (MQ=255)                                                                                                                 ggcATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTaaaa                                 <  2:311936/148‑1 (MQ=255)                                                                                                                      gCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAAttt                                                                                                                      <  8:95911/58‑1 (MQ=255)                                                                                                                      gCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTc                                                                           >  8:155536/1‑101 (MQ=255)                                                                                                                      gCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGt                             >  4:321740/1‑147 (MQ=255)                                                                                                                      gCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGt                             >  6:317795/1‑147 (MQ=255)                                                                                                                      gCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGt                             >  7:39667/1‑147 (MQ=255)                                                                                                                      gCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGAGAGAAAttt                                                                                                                      >  7:95911/1‑58 (MQ=255)                                                                                                                      gCTGGCGGGCCACGAAGAGAGCGGCGGGCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTc                                                                           <  7:155536/101‑1 (MQ=255)                                                                                                                       cTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGtt                            <  8:39667/147‑1 (MQ=255)                                                                                                                       cTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTa                           <  2:311344/148‑1 (MQ=255)                                                                                                                           cGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCt                       >  4:25669/1‑148 (MQ=255)                                                                                                                            gggCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTg                      <  4:133534/148‑1 (MQ=255)                                                                                                                              gCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTTTGCGATGAAACGTCACGTTGGCGGCCTTGCGGAATATCGCGCAGCAGAAGGTAAACCCGTTAAGCTGcc                    >  6:366050/1‑148 (MQ=255)                                                                                                                                 aCGAAGAGAGCGGCGGTTGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAg                         >  8:49846/1‑140 (MQ=255)                                                                                                                                 aCGAAGAGAGCGGCGGTTGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAg                         >  8:43553/1‑140 (MQ=255)                                                                                                                                 aCGAAGAGAGCGGCGGTTGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAg                         <  7:49846/140‑1 (MQ=255)                                                                                                                                 aCGAAGAGAGCGGCGGTTGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAg                         <  7:43553/140‑1 (MQ=255)                                                                                                                                   gAAGAGAGCGGAGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGcagcag                                          >  1:115434/1‑121 (MQ=255)                                                                                                                                   gAAGAGAGCGGAGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGcagcag                                          <  2:115434/121‑1 (MQ=255)                                                                                                                                    acgagagCGGCGGTCGAATCGTGTTGCAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCAAGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCg              <  4:54238/146‑1 (MQ=255)                                                                                                                                    aaGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGATAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCg              <  5:41640/148‑1 (MQ=255)                                                                                                                                    aaGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCg              <  1:275324/148‑1 (MQ=255)                                                                                                                                    aaGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCg              <  8:122848/148‑1 (MQ=255)                                                                                                                                    aaGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCg              <  3:25669/148‑1 (MQ=255)                                                                                                                                    aaGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCg              <  3:321740/148‑1 (MQ=255)                                                                                                                                      gagagCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAg            >  1:235287/1‑148 (MQ=255)                                                                                                                                             gcggTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTaa                          >  2:267217/1‑127 (MQ=255)                                                                                                                                             gcggTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTaa                          <  1:267217/127‑1 (MQ=255)                                                                                                                                                gTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGa  <  7:22351/148‑1 (MQ=255)                                                                                                                                                 tCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGa  >  3:95845/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TGCCGATGCTGCGCACGGTCTGGGCGGAATGATCGTCAGCGATGGTGGCTGCACCACGCCGGGCGATGTGGCGAAAGCCTTTGGCGGCGGTGCCGATTTCGTCATGCTTGGCGGCATGCTGGCGGGCCACGAAGAGAGCGGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGA  >  NC_000913/114058‑114348

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