Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	858,892	C→T	G57D (GGC→GAC)	ybiU ←	DUF1479 domain‑containing protein YbiU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,892	0	C	T	97.7%	134.8 / NA	43	G57D (GGC→GAC)	ybiU	DUF1479 domain‑containing protein YbiU
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (18/24);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTATAAATA  >  NC_000913/858656‑859128                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             cGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGa                                                                                                                                                                                                                                  <  1:302490/249‑1 (MQ=255)    gACGTGAAGCGCTGAGCGACCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAGCCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCg                                                                                                                                                                                                                             <  2:27006/251‑1 (MQ=255)        tGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTg                                                                                                                                                                                                                          >  2:114798/1‑250 (MQ=255)             cgcTGAGACTCCAGAAAAAATTAACGCCGGCGCCATTGTAGAACTCGAAGGAGCAGTTGCGGTCCAGATAATCCAGCATCGACTGAGCCCAGCCTAGGGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAAACGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCAGGACGGCATCGTCCTGTGCATTGAGTGCGTTGATTTCa                                                                                                                                                                                                                    <  2:399326/251‑1 (MQ=255)                   gagTCCCGAAGAAATGATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGTTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGAGGGGTAAATGGCCTTTTATCACATAACAACCGCGACGTTTAATCGTTTCGCGCTGATCTGCAGGAACAAGACCTGTTTTGATATCGGAATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCGTTGAGTGCGTTGATTTCAGCCACt                                                                                                                                                                                                              <  2:318307/249‑1 (MQ=255)                                  tATCGCCGGGGCCGTTGTAGATCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAAGGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCGCTGCAGGAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGAGTTCAGCCACTCGCGTGGCGATGTc                                                                                                                                                                                                <  2:303062/250‑1 (MQ=255)                                       ccGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCt                                                                                                                                                                                          <  2:216546/251‑1 (MQ=255)                                         ggggCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCACACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGAGTGGCAATGTCATCGCtt                                                                                                                                                                                         >  1:144217/1‑250 (MQ=255)                                                ttGTAGCCCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGAt                                                                                                                                                                                  <  1:17545/250‑1 (MQ=255)                                                ttGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTCAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGAtt                                                                                                                                                                                 >  2:133434/1‑251 (MQ=255)                                                          cGTCAAAGCGGTTGCGGTCCAGATAGTCCAGTATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctg                                                                                                                                                                        <  1:54531/250‑1 (MQ=255)                                                                aGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACg                                                                                                                                                                 <  2:223628/251‑1 (MQ=255)                                                                        cGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAACTGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAACGATCTGCTGGCCGTCGCCa                                                                                                                                                           >  2:170593/1‑249 (MQ=255)                                                                        cGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCa                                                                                                                                                           >  2:168383/1‑249 (MQ=255)                                                                            ccAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTgc                                                                                                                                                     <  2:202691/251‑1 (MQ=255)                                                                                aTAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGTGTGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGccc                                                                                                                                                 <  1:435689/251‑1 (MQ=255)                                                                                                        gCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACgg                                                                                                                          >  1:412258/1‑250 (MQ=255)                                                                                                         ccTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGAt                                                                                                                        <  2:176936/251‑1 (MQ=255)                                                                                                                       gcgGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGAt                                                                                                          <  1:114798/251‑1 (MQ=255)                                                                                                                        cgGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTACTCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCACGTGCTTCATCTGACGGATAGCTGCTTTGTGATc                                                                                                         >  1:155275/1‑251 (MQ=255)                                                                                                                                 tGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggcaggca                                                                                                 >  1:326984/1‑250 (MQ=255)                                                                                                                                        ttATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTc                                                                                          <  1:442034/250‑1 (MQ=255)                                                                                                                                                    aaCCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGt                                                                              >  2:392364/1‑250 (MQ=255)                                                                                                                                                        gcgACGTTTAAGCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGTCGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAg                                                                         <  1:9178/251‑1 (MQ=255)                                                                                                                                                         cgACGTTTAATCGGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCGTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGc                                                                        <  2:182763/251‑1 (MQ=255)                                                                                                                                                            cGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCAt                                                                     <  1:428386/251‑1 (MQ=255)                                                                                                                                                                tAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac                                                                  <  1:11389/250‑1 (MQ=255)                                                                                                                                                                           cgcTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGAGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCACGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCAtt                                                       >  2:278946/1‑250 (MQ=255)                                                                                                                                                                               gctcggaaGTAACATGACCTGCTTTGATATCGGTAGAAGACAGCCCCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTGGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAAc                                                    <  2:307966/242‑1 (MQ=255)                                                                                                                                                                               gctCTGCAGTAACATGCCCTGCTTTGATATCGGCATAAGACAGCCCCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACtt                                                  <  1:146587/251‑1 (MQ=255)                                                                                                                                                                                 tctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACAGGCAAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGCTCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAGCTTTc                                                >  1:200622/1‑251 (MQ=255)                                                                                                                                                                                     cAGTAACATGCCCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTACAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTc                                            >  2:400447/1‑251 (MQ=255)                                                                                                                                                                                               cccTGCTTTGGTATCGGCATAAGACAGCCCCCGCCCGACGGCATCGTCCTGTGCTTTGAATGGGTTTATTTCCGCCACCCGTGTGGCAATGGCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTGTCATTCGAAACCATTc                                  <  1:327039/250‑1 (MQ=255)                                                                                                                                                                                                 cTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTg                                >  1:84903/1‑251 (MQ=255)                                                                                                                                                                                                                 aTAAGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGGGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCAGTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAat                <  1:170593/251‑1 (MQ=255)                                                                                                                                                                                                                    aGACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTACCTTTCATTCGAAACCATTCTGATTTGCATGTGCCAatat              >  2:400916/1‑250 (MQ=255)                                                                                                                                                                                                                     gACAGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAAtata             >  1:113718/1‑250 (MQ=255)                                                                                                                                                                                                                        aGCACCGGCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGAAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGtt         <  1:413740/251‑1 (MQ=255)                                                                                                                                                                                                                              ggCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTataaat   >  1:151129/1‑251 (MQ=255)                                                                                                                                                                                                                              ggCCAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCACCCACACGCGGGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAGGCCATAACCACCTCTCATTAAATTTCATTCGAAACCATTCTGATTTACATGTGACAATATAAGTTataaat   >  1:204764/1‑251 (MQ=255)                                                                                                                                                                                                                                ccAGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTataaata  >  1:14644/1‑250 (MQ=255)                                                                                                                                                                                                                                  aGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCttt                                                                                                               >  1:325988/1‑139 (MQ=255)                                                                                                                                                                                                                                  aGACGGCATCGTCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCttt                                                                                                               <  2:325988/139‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             CGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATATAAGTTATAAATA  >  NC_000913/858656‑859128

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