Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,689,746	T→C	16.4%	E172G (GAG→GGG)	bcsZ ←	endo‑1,4‑D‑glucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,689,746	0	T	C	16.4%	130.1 / 2.5	56	E172G (GAG→GGG)	bcsZ	endo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (4/6);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACA  >  NC_000913/3689514‑3689971                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  caTGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCCATAAACGTTGATTGGTTTCGCCCACCGGGGTCCCCGGCGCGCCCAACCGGGTGGAATACTGCGCCCGCGTCGGCGGCAGGTAGCTGGGGGTAAAACACCCACTGTATACCTCCGCaa                                                                                                                                                                                                                             >  2:532557/1‑239 (MQ=255)  aTGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCCAACCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGAAGCTGGGGTTAAAAAGCCAGCCGTTATCCTCAGCaaa                                                                                                                                                                                                                            >  1:707621/1‑239 (MQ=255)  aTGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCCCGCCAAAGCGGGGGAAATACTGCGCCAGCGTCGGCGGCAGGTAACTGGGGTTAAAACGCCAGCCGTTATCCTCAGCaaa                                                                                                                                                                                                                            >  1:568791/1‑239 (MQ=255)        aCGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCCCAGCGTGGTCCACGGCGCGCCAAACCGGGGGAAATACTGCGCCAGCGTGGGCGGCAGGAAGCTGCGGTTAAAACCCCCGCTGTTATCCCCACCAAAAccc                                                                                                                                                                                                                        >  2:663911/1‑237 (MQ=255)            ggATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGATGATTGGTTTCGCGCAGCGTGGTCCCCGGAGCGCCAACACGGGGGAAATACTGCCCCCGCTGCGTCGGCGGATAGACGGGGTAAAAAAGCCAGCGGCAACCCTACACAAAACCCACCTTg                                                                                                                                                                                                                  >  2:288234/1‑239 (MQ=255)               tAGCGTCGTAGCTGGTGATCAATGTTTTTGCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGcc                                                                                                                                                                                                                <  1:347969/238‑1 (MQ=255)                  cgtcgtAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTACCTGGGGTTAAAAAGCCAGCGGTTATCCTCCGCAAAACCCCATTTGccccgg                                                                                                                                                                                                            >  1:36297/1‑236 (MQ=255)                      gtAGCTTCTGTTAAATGTTTTTGCGGCTTTTAGCTGCCCGCCTGTGTCTTTCTCATAGCGCACCCAGTCGGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaaca                                                                                                                                                                                                        <  2:707621/239‑1 (MQ=255)                      gtAACTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATCGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaac                                                                                                                                                                                                         <  1:593986/238‑1 (MQ=255)                            gctgATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAGAACCTTTCGGGGCGGTTTCCCGCAATAAACGTTGATTGGTTTCGCGCACCGTGGTCCACGGCGCGCCAAACGGGGGGAAATACTGGCCCAGCGTCGGCGGCAGGTACCCGGGGGTAAAACGCCAACCGTTTTCCCCCCCAAAACCCACCTTACCCgggcaaaaaaggg                                                                                                                                                                                                  >  2:357305/1‑228 (MQ=255)                            gctgATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGGGGTCCCCGGGCCGCCACAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCCAGTAGCTGGGGGTAAAACGCCAACTGTTTTTCTCCACAAAACCCCCCTTTCCCGGGaaaaaaaaag                                                                                                                                                                                                  >  2:891881/1‑234 (MQ=255)                                 tctattttttttGGTCTTTGTGCTCCCGTCCTTTGTCTTTTTCAATGCGCCCCCGGTCTGGCGAACAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAccc                                                                                                                                                                                             <  2:872730/233‑1 (MQ=255)                                     tGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACCGCGCCAGCGTCGGCGGCAGGGAGCCTGGGGTAAAACGGCCCCTGTTTTCCCCAGCCAAACCCACATTCCCCGGGAAAAACAGGGAAcccagc                                                                                                                                                                                          >  2:411542/1‑238 (MQ=255)                                         tttgctgCTTTTCGTTGCCAGCCTTTGTCTTTCTCAAAGCGCCCCTAGTCTGGGGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGGTTCGCGCAGCGTGGACCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAgg                                                                                                                                                                                     <  2:111819/233‑1 (MQ=255)                                          tttCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATTGGACCCAACCCAGGc                                                                                                                                                                                    >  1:824950/1‑239 (MQ=255)                                             cGGCTTTTGGTTTCCCGCCTTTGTCTTTTCCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCAcc                                                                                                                                                                                 <  2:568791/239‑1 (MQ=255)                                               gCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCCCTTTTCCCGGTAACAAAATGGAAcccagcccaagacacgc                                                                                                                                                                               >  2:882131/1‑231 (MQ=255)                                                 ttttAGCTGCCTGCCTGTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCTGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTTGTCCACGGCGCGCCAAAGCGGGTGAAAGACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGTTAACAACATGGAACCCAGCCCAGGCACCGTca                                                                                                                                                                             <  1:288234/239‑1 (MQ=255)                                                 tttgtgCTGCCCGTTTTTTTTTTTCTCAAGGCGCACCCAGGCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAAAAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTca                                                                                                                                                                             <  2:875465/234‑1 (MQ=255)                                                  tttagcgcacAGCCTGTGTCTTGCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGGTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTca                                                                                                                                                                             <  1:787494/229‑1 (MQ=255)                                                                   gtccttcctaTAGCGCCACCAGTCTGGCGAAAAACCGTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGa                                                                                                                                                           <  1:393690/230‑1 (MQ=255)                                                                   gTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTGCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGa                                                                                                                                                           <  1:637555/239‑1 (MQ=255)                                                                        tctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCCCAGCGTGGGCCACGGCGGGCCAAAACGGGTGAAATACTGCGCCCGCGTCGGCGGCAGGTAGGCGGGGGTCAACAGCCACACGGTTTCCCTACAAAAAACCACCTTTCCCGGGAACAACATGGAACCCAAGCCAGGCACCCGcaccccccacacccccccgacacaa                                                                                                                                                      >  2:20174/1‑218 (MQ=255)                                                                                  cACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGa                                                                                                                                                           <  2:6335/224‑1 (MQ=255)                                                                                  cACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGa                                                                                                                                                           >  1:6335/1‑224 (MQ=255)                                                                                   aCCCAGGCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGGGGGCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAAcgc                                                                                                                                           <  1:905383/239‑1 (MQ=255)                                                                                         tCTGGCGAAAAACCTTTGCGGGCGTTTTCCAGAAATAAGCGTTGATGGTGTGCGCGCATGGTTGGCCCCGGCGGGCCAAAGCGGGTGAAATACTGGGCCAGCTTCGGCGGCAGGGAGCTGGGGTTAAAACGCCAGCTGTTAGCCTCAGCAAAACCCACTTTGGCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGcc                                                                                                                                     <  2:272886/239‑1 (MQ=255)                                                                                         tCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGGAACACCATGGAACCCAACCCAGGCACCCTCACCACTTCCTCCCGCCCCATACGATTTAGCAACCCGCTGcc                                                                                                                                     >  1:240575/1‑239 (MQ=255)                                                                                          cTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGGTTCGCGCCGCGTGGTCCACGGCGCGCCGAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAACAGCCCGCTGTTATCGTCAGCAAAAACCACTTTTCCCGCTAACAACATGGAACCCCACCCAGGCACCGGCCCCACCTCCTCCCCCCCGaaaaggtttgacaacacgacgcca                                                                                                                                    >  2:151103/1‑216 (MQ=255)                                                                                            ggCGAAAATCCTTTCGGGGCGGTTTCAAGCAATAAACGTTGATTGGTTTGGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGAt                                                                                                                                  <  1:20174/239‑1 (MQ=255)                                                                                               gAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGGAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGt                                                                                                                                <  1:882131/238‑1 (MQ=255)                                                                                                   aaCCTTTCGGGGCGGTTTCCCGCAATAAGCGGTGATTTGTTTCGCGCAGCGTGGTCCACGGGGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGta                                                                                                                           <  2:36297/239‑1 (MQ=255)                                                                                                   aaCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAACGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTCGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCCCTTTGCCCGGGAACAACATGGAACCCCGCCCCGGCACCCGCACCCCTTCCCCCCCCGCGAAACATTTTACCAAcccccccccgacgccggga                                                                                                                           >  2:350837/1‑221 (MQ=255)                                                                                                                       aGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACCCGCTGCCGAGTTCGGTATAACCCTGCCCTTTCCCCaa                                                                                                       >  1:443701/1‑239 (MQ=255)                                                                                                                            taaccttttgtgggTTTCGCGCGGGGGGGGCCCCGGGCCTCTAGAGCGGGTAAAATCCTGCGCCAGGGTCGGGTGCCGGTAGCGGGGGTTAAAACCCCAGTGTTTACCCCCAGCAAACCCCACTTTGCCCGGTAAAAACATGGACCCCAGCCAAGGCCCCTTCCCCTCTTCCTCCCGCGCGATCCGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAACCGcc                                                                                                  <  2:8229/227‑1 (MQ=255)                                                                                                                              aaCGTTGATTGGTTCCGCGCGTCGTGGCCCACGGCGCGCCGAAGCGGGTGAAATCCTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGcccc                                                                                                <  1:151103/239‑1 (MQ=255)                                                                                                                                               cgcAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCa                                                                               >  1:487216/1‑239 (MQ=255)                                                                                                                                                    cGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCCCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCCGCAACGACCAGccc                                                                           >  1:439287/1‑238 (MQ=255)                                                                                                                                                        gTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcc                                                                                                                                                                                              <  2:394343/119‑1 (MQ=255)                                                                                                                                                        gTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcc                                                                                                                                                                                              >  1:394343/1‑119 (MQ=255)                                                                                                                                                              ggCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCCGACAt                                                                 >  1:43066/1‑238 (MQ=255)                                                                                                                                                                  cgcCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATCCGTTTTAGCAACGCCCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCCCCAGCCACACCCAAGCCATCCCGAcaccacaa                                                            >  2:202182/1‑237 (MQ=255)                                                                                                                                                                  cgcCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGCCCAGGCCATCCAGACATAAcc                                                             >  1:713364/1‑238 (MQ=255)                                                                                                                                                                          gggggGAAATACCGCCCCGGGTTGGGCGTCGGGGAGGTGGGTTTAAAGCGGCGACGGTTAGCCCTAGAAAACCCCACTGTGCCAGGTACCAAAATGGAACCCAGCCGCGGCACGGTCACCACTTCCTCCCGGCCGATTCGTTTTAGCACCCCGCTGCCGATGGCGGTATAACGCTGCTCTTTCAACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGa                                                       <  2:130982/235‑1 (MQ=255)                                                                                                                                                                               gAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCGCTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCCCAAACGCCCCGCCTCCACCAACGACCAGCCCATCCAAACATACCCACCAGAGGCCgcaa                                               >  2:358175/1‑236 (MQ=255)                                                                                                                                                                                 aagaCTGCGCCAGCGCCGGCGGAAGGTAGCGGGGGTTAAAACCCCAGCTGTTAACCTCAGCAAAACCCACTTTGCCCGGTAACAACAGGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAACCGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTg                                             <  1:357305/236‑1 (MQ=255)                                                                                                                                                                                       gcgcCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACt                                                                                                                                                                        <  2:412779/110‑1 (MQ=255)                                                                                                                                                                                       gcgcCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACt                                                                                                                                                                        >  1:412779/1‑110 (MQ=255)                                                                                                                                                                                            acggTCGGCGGCAGGTACCTGGGGTTAAAACGCCAGCTGTTAACCCAAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCAc                                  <  1:202182/236‑1 (MQ=255)                                                                                                                                                                                               ggcggcggCAGGGAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGACCCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACtt                                <  1:663911/236‑1 (MQ=255)                                                                                                                                                                                               gTCGGCGGCAGTGAGCTGGGGTTAAAACGCCAGCTGTTACCCCCAGCAAAACCCCCTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCAc                                  <  1:411542/236‑1 (MQ=255)                                                                                                                                                                                                 gggcggtgggttgtcggggTTAAACCCCCTGGTTTTACCCCCAGCAAACCCCCCTTTGCCGGGTAACAACAGGGAACCCACCCCAGGCACCGTCACCACTTCCGCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTccc                             <  2:443701/224‑1 (MQ=255)                                                                                                                                                                                                        gggttgttGGGGTTAAAACCCCAGCTGTTACCCCCAGCAAAACCCCCTTTCCCCGGAAACAACATGGAACCCAGCCCAGGCCCCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATACCGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTAc                       <  2:439287/231‑1 (MQ=255)                                                                                                                                                                                                                             cAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCCCTTACTTTTCCCTTTCTTTCCCcacc  >  1:309480/1‑237 (MQ=255)                                                                                                                                                                                                                              aGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTc                                        >  2:851543/1‑199 (MQ=255)                                                                                                                                                                                                                              aGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTc                                        <  1:851543/199‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  CATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCCCACA  >  NC_000913/3689514‑3689971

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