Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	803,662	C→A	L54I (CTT→ATT)	ybhJ →	putative hydratase YbhJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	803,662	0	C	A	97.4%	145.1 / NA	39	L54I (CTT→ATT)	ybhJ	putative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (18/20);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAAT  >  NC_000913/803422‑803899                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              ccTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATAGGGATAAaattaaaattaa                                                                                                                                                                                                                                     <  2:576404/251‑1 (MQ=255)    gCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGtt                                                                                                                                                                                                                                  <  2:231019/251‑1 (MQ=255)         tCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATaaaattaaaatt                                                                                                                                                                                                                                       <  1:478039/241‑1 (MQ=255)          cGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                                                                                                                                                                                                            <  2:374997/251‑1 (MQ=255)            aaTTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTATCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAt                                                                                                                                                                                                                          >  2:626787/1‑251 (MQ=255)                attTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                                                                                                                                                                                                      <  1:184360/251‑1 (MQ=255)                  ttAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGcc                                                                                                                                                                                                                     <  2:99650/250‑1 (MQ=255)                    agagTATATCCATGTATTATCTTTCTGCGCACTTCCCGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCTGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                                                                                                                                                                                                   <  2:47915/250‑1 (MQ=255)                     gagTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATGCATTAGCCtctc                                                                                                                                                                                                                 <  2:35113/251‑1 (MQ=255)                                                                        ggAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCGCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATgg                                                                                                                                                              <  2:630193/251‑1 (MQ=255)                                                                                                 taaaaaGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACTTCCGGAAATATGGATAAAATTAGAAGTAAGTGTGAGTCATTAGCCGCTCACGATATTGCCTGTGTTGGTATGGTTCAGTCCGCGAAAGCGTGCGGTATGGAACGTTTCCCGCTGCCGTAtgt                                                                                                                                        <  1:88305/247‑1 (MQ=255)                                                                                                    aaaGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGCAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGAc                                                                                                                                  >  2:480948/1‑251 (MQ=255)                                                                                                                   gCCAGTAATAACTAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCTTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAactc                                                                                                                   <  1:576485/251‑1 (MQ=255)                                                                                                                   gCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGTAATTAAAAAAGAAGAAGCGAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAACGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAactc                                                                                                                   >  1:287638/1‑251 (MQ=255)                                                                                                                          ataaCGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAACAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGACTATTCTCTCTTCGCATAATACGTCCGGAACTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGATAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTgc                                                                                                            >  1:734761/1‑251 (MQ=255)                                                                                                                                aaaTAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCGCTCACGATATTGCCTTTGTTTGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAAGGTGTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTc                                                                                                      <  1:348063/251‑1 (MQ=255)                                                                                                                                    aatTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTcggcg                                                                                                  >  1:1831/1‑251 (MQ=255)                                                                                                                                        gCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCTTCGGCGGCa                                                                                               <  1:236637/250‑1 (MQ=255)                                                                                                                                                  aTTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGt                                                                                    <  2:616481/251‑1 (MQ=255)                                                                                                                                                          ggCGAAATTACAAAAGAAGAAGCCAAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGCAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTCCCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGACCGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTCTTAACTGtgatgccc                                                                             >  2:130599/1‑248 (MQ=255)                                                                                                                                                          ggCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGAccc                                                                            >  2:189185/1‑250 (MQ=255)                                                                                                                                                            cGAAATTAACAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATg                                                                          >  1:215081/1‑251 (MQ=255)                                                                                                                                                             gAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGCAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCCCTATTAACGGGTATGACCATGt                                                                         >  2:244507/1‑251 (MQ=255)                                                                                                                                                             gAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGt                                                                         >  2:439091/1‑251 (MQ=255)                                                                                                                                                                             aagCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGa                                                                                                                                                             >  2:70171/1‑151 (MQ=255)                                                                                                                                                                             aagCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATGACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGa                                                                                                                                                             <  1:70171/151‑1 (MQ=255)                                                                                                                                                                                 cAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGCTTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGccc                                                     >  1:627461/1‑251 (MQ=255)                                                                                                                                                                                      aa‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGt                                                >  1:185567/1‑251 (MQ=255)                                                                                                                                                                                       a‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGtt                                               <  1:189185/251‑1 (MQ=255)                                                                                                                                                                                             cTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGTGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATggg                                          >  1:198815/1‑250 (MQ=255)                                                                                                                                                                                                 tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCgggg                                      >  2:212981/1‑249 (MQ=255)                                                                                                                                                                                                  gCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGttt                                     >  1:504230/1‑249 (MQ=255)                                                                                                                                                                                                        tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCGGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTAGTAACGGTGATGGCCAGGTTTTTGGTTTATcggcg                                                         <  2:228605/225‑1 (MQ=255)                                                                                                                                                                                                           aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCCGCGTTATGGCGGTATTTTTTTGcc                            >  1:703921/1‑251 (MQ=255)                                                                                                                                                                                                           aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTtgt                               <  1:168886/248‑1 (MQ=255)                                                                                                                                                                                                           aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTtgt                               <  1:480948/248‑1 (MQ=255)                                                                                                                                                                                                                ctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCCCGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCGGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTAGTTTTGTGCCTcc                         <  2:379857/249‑1 (MQ=255)                                                                                                                                                                                                                 tctTCGCATAATACGTACGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTtcctccccc                      >  1:511680/1‑248 (MQ=255)                                                                                                                                                                                                                    tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGGGCCTCCGCata                    <  2:177409/250‑1 (MQ=255)                                                                                                                                                                                                                                  cGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGGCGCCGCATATTgtggggatcactc     >  1:603404/1‑239 (MQ=255)                                                                                                                                                                                                                                     aaaTATGGATAAAATTAAAATTAAGTTTGATTAATTAGCCTCTCACGATATTACCTTTGTTTGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGGCGGCGGCACTATTAACGGTGATGACCATGTTTGTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAat  <  2:703921/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              CCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAA‑GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAAT  >  NC_000913/803422‑803899

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