Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	749,268	C→T	G147E (GGA→GAA)	ybgP ←	putative fimbrial chaperone YbgP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	749,268	0	C	T	97.4%	123.3 / NA	38	G147E (GGA→GAA)	ybgP	putative fimbrial chaperone YbgP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/16);  minor base G (0/1);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTCGTCATGCGTAATGCGCCGTAATCATCCATGTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTCCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCGCTTTCGGTTCCAGACGTTGAATAGGGGGCAAAGCGACCAGAAGAGCA  >  NC_000913/749025‑749510                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   aGCGTCGTCATGCGTAATGCGCCGTAATCATCCATGTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCttt                                                                                                                                                                                                                                             >  1:326693/1‑251 (MQ=255)                aaTGCGCCGTCATCATCCATGTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGGACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGCCCTTTTCTTCTGCTTTCTTGCGTAATGCAGc                                                                                                                                                                                                                              >  2:388945/1‑251 (MQ=255)                aaTGCGCCGTAATCATCCATGTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTGGCGTAATGCAg                                                                                                                                                                                                                               <  2:444556/250‑1 (MQ=255)                   gcgcCGTAATCATCCATGTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCgg                                                                                                                                                                                                                           <  1:272808/251‑1 (MQ=255)                      ccGTAATCATCCATGTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGTTCg                                                                                                                                                                                                                        <  2:574919/251‑1 (MQ=255)                         tAATCATCCACGTAACC‑AGAAAACAACTGGCTGCCGCTATAAATTCAGGTGTTTGTGTTGAAGGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGTCGACCGCCTTTTTCCTTTCTCCCCAGATAGGCAATGGTCTGATAATACGAGGTAGGATTTTTCAGTGTCAGTCGGCTGCCCTGATGGCTGAACTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGcc                                                                                                                                                                                                                      >  2:257854/1‑250 (MQ=255)                             catcCATGTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAAc                                                                                                                                                                                                                 >  2:455675/1‑251 (MQ=255)                                ccATGTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCTGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTTTTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAAt                                                                                                                                                                                                              <  1:32018/251‑1 (MQ=255)                                    gTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTa                                                                                                                                                                                                           >  1:625178/1‑250 (MQ=255)                                        ccGAGATAA‑AACTTGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGTCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTTCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCACCAGAACAATTTAATAc                                                                                                                                                                                                       <  2:196817/250‑1 (MQ=255)                                                                ataACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGt                                                                                                                                                                                >  2:505495/1‑250 (MQ=255)                                                                   aCTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGACTAGCGACCTGAAGTATCg                                                                                                                                                                            >  2:432188/1‑251 (MQ=255)                                                                                    tGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCtttt                                                                                                                                                            >  2:310022/1‑250 (MQ=255)                                                                                           aCTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGa                                                                                                                                                     <  1:578421/250‑1 (MQ=255)                                                                                                               cGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCTGGAGCGGGAGGAATTTCgcgcttc                                                                                                                                >  2:20070/1‑248 (MQ=255)                                                                                                                     tGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTa                                                                                                                           <  2:317082/250‑1 (MQ=255)                                                                                                                                  aaCGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAGTTTCGCGCATATTGTAGAAAAATAGCGttt                                                                                                             >  2:484243/1‑251 (MQ=255)                                                                                                                                   aCGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTGCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTc                                                                                                            <  2:326693/251‑1 (MQ=255)                                                                                                                                           ttCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGc                                                                                                    >  1:113513/1‑251 (MQ=255)                                                                                                                                              aTTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGgga                                                                                                                                               <  1:608661/205‑1 (MQ=255)                                                                                                                                              aTTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGgga                                                                                                                                               >  2:608661/1‑205 (MQ=255)                                                                                                                                                    ccccAGATAGGCAATGGTCAGATAAGACGCGTTAGGATTTTTCAGTGTCAGTTGGTTGCCCCGCTGGCTGACCGGTAACTGCATTTCGGCCTTTTCTTCTGCTTTCTTGCGTAATGCGGCCGGTCGCCAGAACAATTGAATACGTGTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCGTTTATCCGGAGCGGTAGGAATTGCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTa                                                                                              <  1:622367/248‑1 (MQ=255)                                                                                                                                                             ggCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTg                                                                                  >  1:288944/1‑251 (MQ=255)                                                                                                                                                             ggCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTg                                                                                  >  1:369567/1‑251 (MQ=255)                                                                                                                                                                      cAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGtt                                                                         >  1:422246/1‑251 (MQ=255)                                                                                                                                                                                              cAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCg                                                 >  2:650744/1‑251 (MQ=255)                                                                                                                                                                                              cAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCg                                                 >  1:364102/1‑251 (MQ=255)                                                                                                                                                                                               aGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCAGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCGc                                                >  1:404160/1‑251 (MQ=255)                                                                                                                                                                                                                     ggctgACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCGCTTTCGGTTCCAGACGTTGaa                            <  1:455675/249‑1 (MQ=255)                                                                                                                                                                                                                      gctgACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTGCGCTCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGCTTCACAACCCGCACCTGCGATGTCGCTTTCGGTTCCAGACGTTGAATAg                         <  1:144836/251‑1 (MQ=255)                                                                                                                                                                                                                        tgACCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTAGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCGCTTTCGGTTCCAGACGTTGAATAggg                       >  2:156236/1‑251 (MQ=255)                                                                                                                                                                                                                          aCCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATg                                                    >  1:386133/1‑220 (MQ=255)                                                                                                                                                                                                                          aCCTGTAACTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATg                                                    <  2:386133/220‑1 (MQ=255)                                                                                                                                                                                                                                aaCTGCAGTTCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTGTCACAACCCGCACCTGCGATGTCGCTTTCGGTTCCAGACGTTGAATAGGGGGCAAAGc               <  1:388945/251‑1 (MQ=255)                                                                                                                                                                                                                                         tCGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCGCTTTCGGTTCCAGACGTTGAATAGGGGGCAAAGCGACCagaa       <  2:224474/250‑1 (MQ=255)                                                                                                                                                                                                                                          cGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCGCTTTCGGTTCCAGACGTTGAATAGGGGGCAAAGCGACCagaaga     <  1:310022/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGACCTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCGCTTTCGGCTCCAGACGTTGAATAGGGGGCAAAGCGACCaga        >  1:106091/1‑248 (MQ=255)                                                                                                                                                                                                                                             ccTTTTCTTCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCGCTTTCGGTTCCAGACGTTGAATAGGGGGCAAAGCGACCAGAAGAGCa  <  1:505495/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   AGCGTCGTCATGCGTAATGCGCCGTAATCATCCATGTAACCGAGATAA‑AACTGGCTGCCGCTATAACTTCCGGTGTTTGTGTTGACCGTACTAAAAGGTGCAACCATCACGGTTTTGAAGCCTGGGAGAACGCCTTTTTCATTTCTCCCCAGATAGGCAATGGTCAGATAATACGCGGTAGGATTTTTCAGTGTCAGTTGGTTGCCCTGCTGGCTGACCTGTAACTGCAGTTCGACCTTTTCTCCTGCTTTCTTGCGTAATGCAGCCGGTCGCCAGAACAATTTAATACGACTTTGAATAGCGACCTGAAGTATCGCATGGTCGCTGCTTTTATCCGGAGCGGGAGGAATTTCGCGCATATTGTAGAAAAATAGCGTTTCGCGATCGCCAGGTAATTGTGTGGTTGAGGCTTGTTTCACAACCCGCACCTGCGATGTCGCTTTCGGTTCCAGACGTTGAATAGGGGGCAAAGCGACCAGAAGAGCA  >  NC_000913/749025‑749510

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