Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	386,290	G→A	E28E (GAG→GAA)	tauB →	taurine ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	386,290	0	G	A	97.4%	129.3 / NA	39	E28E (GAG→GAA)	tauB	taurine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/21);  minor base C (1/0);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACGCCGCAGCAACAAACGGCAGAACTGACCGGACCGGTGAACAAAGCGATCATCGACACCGCGCAGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCAC  >  NC_000913/386048‑386536                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       cTGACGCCGCAGCAACAAACGGCAGAACTGACCGGACCGGTGAACAAAGCGATCATCGACACCGCGCAGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACtggg                                                                                                                                                                                                                                                >  1:212655/1‑250 (MQ=255)   gACGCCGCAGCAACAAACGGCAGAACTGACCGGACCGGTGAACAAAGCGATCATCGACACCGCGCAGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACtgg                                                                                                                                                                                                                                                 <  1:41306/248‑1 (MQ=255)                          aaCTGACCGGACCGGTGAACAAAGCGATCATCGACACCGCGCAGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACTCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCGGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGt                                                                                                                                                                                                                       <  1:395103/251‑1 (MQ=255)                                   gaccggTGAAAAAAGCGATCATCGACACCGCCCAGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACCGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAaccacc                                                                                                                                                                                                              >  1:360729/1‑251 (MQ=255)                                     ccggTGAACAAAGCGATCATCGACACCGCGCAGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCct                                                                                                                                                                                                            <  2:286367/251‑1 (MQ=255)                                       ggTGAACAAAGCGATCATCGACACCGCGCAGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCctgc                                                                                                                                                                                                          >  2:41597/1‑251 (MQ=255)                                                       atcGACACCGCGCAGTTTTTGAAAGAGCAGGGCAAGGTCCCTGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTGCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCGCTCATCGTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCGGACCCTTGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCt                                                                                                                                                                                                                            <  1:249211/217‑1 (MQ=255)                                                       atcGACACCGCGCAGTTTTTGAAAGAGCAGCGCAGGGTCCCGGCTGTAGCGAATGATTACAGCAAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCCAACCGGCACTAGAAGATATCAAACTGACGCTGGAAACCGGCGAACTACTGGTGGTGCTGGGGCCCTCCGGCt                                                                                                                                                                                                                            >  2:249211/1‑217 (MQ=255)                                                              ccGCGCAGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTtgtg                                                                                                                                                                                   <  1:26706/251‑1 (MQ=255)                                                                    aGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTAt                                                                                                                                                                             >  1:312633/1‑251 (MQ=255)                                                                        tttGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAg                                                                                                                                                                          <  2:99285/250‑1 (MQ=255)                                                                          tGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGCTTGCCGGTTTTGTGCCTTATCAGCAt                                                                                                                                                                       >  1:237352/1‑251 (MQ=255)                                                                             ccgagCAGGGCCAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCGGGAAAACGGCGACCTACTGGGGGGGTTGTGGCCGTTCGGCTGCGGTAAAAACACCCTGCTGAATCTGATTGCCGGTTTTTTGCCTTATccgcctgtc                                                                                                                                                                    >  2:213510/3‑243 (MQ=255)                                                                                  cagtgCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGTACTGGAAGATATCAACCTGAAGCTGGAAAGAGGGGAACTACTGGTGGGGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGGCTTATCAGCATgtaagcat                                                                                                                                                               >  2:291046/5‑244 (MQ=255)                                                                                      gCAAGGTCCCGGCTGTAGAGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGGAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCaa                                                                                                                                                           <  2:284742/251‑1 (MQ=255)                                                                                                ggCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTTTCAGCATGGCAGCATTCAACTGGCggggg                                                                                                                                                 >  2:44767/1‑249 (MQ=255)                                                                                                ggCTGTAGCGAATGATTACAGCCAGTACGGTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCATCGGGCGGGTa                                                                                                                                                 <  1:77458/251‑1 (MQ=255)                                                                                                       gCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGttt                                                                                                                                          >  2:381115/1‑249 (MQ=255)                                                                                                                  cAGCCAGTACGTTACCTCTCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGG‑AAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCggg                                                                                                                              <  2:383101/251‑1 (MQ=255)                                                                                                                              tACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGt                                                                                                                    <  1:347525/250‑1 (MQ=255)                                                                                                                                 cTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAAACTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGTGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCTGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCg                                                                                                                >  1:261944/1‑251 (MQ=255)                                                                                                                                     cgcTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGCAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGAGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCAGGTTTTGTGCCTTATCAGCCTGGCAGCATTCAACTGGCGGGGAAGCGTATTGAGGGACCGGGAGCGGAGCGGGGCGtcgt                                                                                                            >  1:255507/1‑248 (MQ=255)                                                                                                                                             gCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGTGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCTGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAg                                                                                                     <  2:261944/250‑1 (MQ=255)                                                                                                                                               aaTAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGGGGCGTAGTTTTTCAggat                                                                                                  >  1:19386/1‑248 (MQ=255)                                                                                                                                                            gcAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCtactac                                                                                     <  2:312633/251‑1 (MQ=255)                                                                                                                                                                  gctgcAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGgc                                                                               <  2:237352/251‑1 (MQ=255)                                                                                                                                                                          tctctcATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGt                                                                        >  1:85538/1‑250 (MQ=255)                                                                                                                                                                             ctcATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTa                                                                                                                                                                                                                      >  1:290454/1‑105 (MQ=255)                                                                                                                                                                             ctcATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTa                                                                                                                                                                                                                      <  2:290454/105‑1 (MQ=255)                                                                                                                                                                             ctcATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGt                                                                        <  1:153429/247‑1 (MQ=255)                                                                                                                                                                                  cTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAAc                                                                                                                                                          <  2:360253/160‑1 (MQ=255)                                                                                                                                                                                  cTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAAc                                                                                                                                                          >  1:360253/1‑160 (MQ=255)                                                                                                                                                                                      aCGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTg                                                           <  2:118948/251‑1 (MQ=255)                                                                                                                                                                                                                 tGGAAGATATCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGAAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGGCCGGGAGCAGAGGGTGGCGGAGGTTTTCAGAATGGGGGGCTACTACCGTGGCGCAATGTGCAGGAACACGCAGGGTTCGGCCCGCAATTGGGAggcatt                                >  2:304857/1‑247 (MQ=255)                                                                                                                                                                                                                     aGATAGCAACCTGACGCTGGAAAGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGTACCGGGAGCAGAGCGTGGCGTAGTTTTCCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAgaga                            <  2:212655/251‑1 (MQ=255)                                                                                                                                                                                                                                           aGCGGCGAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATcgc      <  1:404023/251‑1 (MQ=255)                                                                                                                                                                                                                                            gcggcgAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATcgcg     <  1:336879/251‑1 (MQ=255)                                                                                                                                                                                                                                              ggcgAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCgcgcc   >  1:132235/1‑250 (MQ=255)                                                                                                                                                                                                                                               gcgAACTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCAc  <  1:304857/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       CTGACGCCGCAGCAACAAACGGCAGAACTGACCGGACCGGTGAACAAAGCGATCATCGACACCGCGCAGTTTTTGAAAGAGCAGGGCAAGGTCCCGGCTGTAGCGAATGATTACAGCCAGTACGTTACCTCGCGCTTCGTGCAATAAAAGGAGGCGCAGATGCTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCAC  >  NC_000913/386048‑386536

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