Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,689,746	T→C	16.5%	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.5%	158.8 / 3.2	68	E172G (GAG→GGG)	bcsZ	endo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand):  ref base T (21/34);  new base C (5/7);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCC  >  NC_000913/3689509‑3689967                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                aCCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCAACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGCAGCCGGGGTTAAAATGCCAGCTGTTATCC‑Tc                                                                                                                                                                                                                                >  1:160218/1‑239 (MQ=255)       aTGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGGGGTCCCCGGCGCGCCAAAGCGGGGGAAATACTGGGCCCGCGTCCGCCGGAGGGAGCCGGGGGTTAAAAGCCAACTTGTATCC‑CCA‑‑Caaaa                                                                                                                                                                                                                        >  2:312026/1‑239 (MQ=255)            aaCGCTGATAGCGTCGTAGCTGCTGATCACTGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTTATTGGTTTCGCGCCGGCGTGG‑CCACGGGGCGCCCAAGCGGGTGAAAAACTTGGCCAGCGTCCGCGGCACGGAGCCGGGGGTAAAACACCAGCTGCTCTCC‑CCAC‑AAAAA‑acc                                                                                                                                                                                                                    >  2:35694/1‑238 (MQ=255)            aaCGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATAATGCGCCAGCGTCGGCGGCCGGGAGCCGGGGTTAAAAACCCAAGTGTTATCC‑CCAG‑CAAAA‑ccc                                                                                                                                                                                                                    >  1:592472/1‑238 (MQ=255)                 ggATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTTAAACGCCAGCTGGTATCC‑TCAG‑ACAAA‑CCCCCttt                                                                                                                                                                                                               >  1:298776/1‑238 (MQ=255)                 ggATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGGTTCCCG‑CAGCGGGGTCCACGGCGCGCCAAAGCGGGTGCAATACTGCGCCAGCGTCGGCGGCCGGTAGCCGGGGGTAAAAACCCAGCTGTTATAC‑CCAG‑CAAAA‑CACACTTTg                                                                                                                                                                                                              >  2:431369/1‑239 (MQ=255)                          cttttgttcttaacatttttttttCGGCTTTTAGCTCCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTaa                                                                                                                                                                                                      <  2:629131/220‑1 (MQ=255)                          cgtAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCGTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTaa                                                                                                                                                                                                      <  1:1027978/238‑1 (MQ=255)                           gttgcttattaTCATTTTTTTTTCGGCTTTGAACTGCCAGCCTGTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATTCTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTaaca                                                                                                                                                                                                    <  2:929367/229‑1 (MQ=255)                           gtggctgcggaacatgtgtttttCGCGTTTTAGCGTCCAACCTTTGGCTTTCCCATACCGCACCCATGCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTGTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTaaca                                                                                                                                                                                                    <  2:277299/221‑1 (MQ=255)                               cctttgttCAATTTTTTTTCGGATTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTTCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTATCAACAt                                                                                                                                                                                                <  2:592472/232‑1 (MQ=255)                                tgctgATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCTCCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACAt                                                                                                                                                                                                <  1:503630/238‑1 (MQ=255)                                             tttttCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAAcccagcccag                                                                                                                                                                                  <  1:42160/239‑1 (MQ=255)                                             tttttCGGCTTTTACGTGCCCGCCTTTGTCTTTCTCAAAACGCACCCAGTCGGGCGAAAAACCGTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAAcccagcccag                                                                                                                                                                                  <  2:298776/239‑1 (MQ=255)                                              ttttCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAgg                                                                                                                                                                                 <  1:10236/239‑1 (MQ=255)                                               tttCGGCTGTTACCCGCCAGCCTTTGTCTTTCTCAGAGCGCACCCAGGCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGc                                                                                                                                                                                <  1:312026/239‑1 (MQ=255)                                                    gCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTATCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑t                                                                                                                                                                                                                                 <  2:653/187‑1 (MQ=255)                                                    gCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTATCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑t                                                                                                                                                                                                                                 >  1:653/1‑187 (MQ=255)                                                      tttttgcgccccgcttttGTTTTTCCCATAGCGCCCCCAGTCTGGCGAAAACCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTTGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTca                                                                                                                                                                         <  2:1021308/224‑1 (MQ=255)                                                      ttttAGCTTCCCGCCTTTTTTTTTCTCACAGCCCACCCGGTCTGGCGAAAAACCTTTCGGGGCCGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGGAGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTca                                                                                                                                                                         <  2:160218/239‑1 (MQ=255)                                                      ttttAGCCGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTc                                                                                                                                                                          <  1:547168/238‑1 (MQ=255)                                                      ttttACCTGCCTGCCTTTGTCTTTCCCATAGCCCCCCCAGGCTGGCGAAAACCCCTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTTGTCCACGGCGCACCAAAGCGGGTGAAATACTGCGCCAGCGGCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑CCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTca                                                                                                                                                                         <  2:1066879/239‑1 (MQ=255)                                                                cAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAACCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTtcctc                                                                                                                                                               <  1:671354/239‑1 (MQ=255)                                                                   aagggaacgtttccaaggggcACCCAGTTTGGCTAAAGCCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCcg                                                                                                                                                            <  2:93656/220‑1 (MQ=255)                                                                     ttttctttttccATAGCCCACCCCGCTTGGCGAAAACCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGGGGTCCACGGCGCGCCAAAGCGGGTTAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCcgcg                                                                                                                                                          <  2:793328/228‑1 (MQ=255)                                                                       tGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGGGAAAAACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTTGTATCC‑TCCGCCAAAA‑CCCCCTTTCCCCGGAACAAAAAGGGAACACAACCAAGG‑ACCCCCACCCCCTCCCCCcacgcg                                                                                                                                                        >  2:103261/1‑239 (MQ=255)                                                                        gTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACCTCCTCCCCCCCGa                                                                                                                                                       >  1:589425/1‑239 (MQ=255)                                                                              ctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTTGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTACCCGGGGTTAAAACGCCAACTGGTATCC‑CCAG‑CAAAA‑CCCACTTTGCCCGGTAACAAAATGGGACCCACCCCAGGCACCCTCCCCAATTCCCCCCCCCCCAAACGtt                                                                                                                                                 >  2:784065/1‑239 (MQ=255)                                                                              ctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGAAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAAAAACATGGAACCCAGCCCAGCACCCGGCACAACTTCCCCCCCCCCGAAACGtt                                                                                                                                                 >  2:359957/1‑239 (MQ=255)                                                                                       cACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCCCTTCCTCCCGCGCGCTACGttctgcccctt                                                                                                                                        >  1:93656/1‑230 (MQ=255)                                                                                       cACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGACCCGTCACCCCTTCCCCCCCCCCGATCCGTTTTGCCAAc                                                                                                                                         >  1:579349/1‑238 (MQ=255)                                                                                             gTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGCTACGTTTTAGCAAACCGCTg                                                                                                                                   >  1:686365/1‑238 (MQ=255)                                                                                               cTGGCGAAAACCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTg                                                                                                                                   <  1:768490/236‑1 (MQ=255)                                                                                               cTGGCGAAAACCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑AAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGcc                                                                                                                                 <  1:65624/238‑1 (MQ=255)                                                                                                 ggCGAAACCCCTTTCGGGGGGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGAt                                                                                                                              <  2:589425/239‑1 (MQ=255)                                                                                                  gCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGGTAAAACGCCAGCTGTTATCC‑TCCG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGAACCGCCACCCCTTCCCCCCCCCCGCTACGTTTTAGCaaaccgcaccccgt                                                                                                                              >  2:65624/1‑226 (MQ=255)                                                                                                    gAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGt                                                                                                                            <  1:431369/238‑1 (MQ=255)                                                                                                    caaaaaCCCTTGGGGGCGGGTTCCCGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTTAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGt                                                                                                                            <  1:1061277/237‑1 (MQ=255)                                                                                                        accccttCGGGGCGGGTTCCAGGAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGta                                                                                                                       <  1:170953/234‑1 (MQ=255)                                                                                                             ttCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACg                                                                                                                  <  1:716472/239‑1 (MQ=255)                                                                                                                           ccgCAATAAACGTTGATTGGTTTCGCG‑CGGGGTGGTCCACGGCGCGCAAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCcaca                                                                                                    <  2:559829/237‑1 (MQ=255)                                                                                                                                 tAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGGCGGTCTAACGCTGCTCTTTCCACAAACGc                                                                                               >  1:21731/1‑238 (MQ=255)                                                                                                                                         tttttgttttccgcgaggttGGTCCCGGGCGCGCCAAGGCGGGTGAAATCCTGCGCCAGCGGGGGGGGCAGGTAGCTGGGGTTAAACCGCCAGCTGTTACCC‑TCAG‑CAAAC‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCGGGCCCCGTCACCACTTCCTCCCGCGCGATCCGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTc                                                                                       <  2:587007/220‑1 (MQ=255)                                                                                                                                         gATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGa                                                                                                                                                       >  1:280125/1‑174 (MQ=255)                                                                                                                                         gATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGa                                                                                                                                                       <  2:280125/174‑1 (MQ=255)                                                                                                                                         gATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCCCTTTCCACAAACGCCCCGCCTc                                                                                       >  1:688069/1‑238 (MQ=255)                                                                                                                                            tGGGTTCGCG‑CCGGGTGTGCCACCGCGGCCAAAAGGGGGGGAAATACTGCCCCAGGGGCGGCGGCAGGTAGCTGGGGGTAAAACCCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCCCTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCACCGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGc                                                                                   <  2:686365/239‑1 (MQ=255)                                                                                                                                              tgtttcggcgcgCGTGGTCCAGGGCCCCCCAAGGGGGGGGAAATACTGGGCCAGCGTCGGCGGCAGGTAGCGGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCCCTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCCCCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCa                                                                                  <  2:579349/228‑1 (MQ=255)                                                                                                                                               tttCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATACCGCTGCCCTTTCCAACAACGCCCCGCCCCCAGCCAc                                                                                >  1:1059178/1‑239 (MQ=255)                                                                                                                                                    cg‑cAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTCCTCTTTCCAAAAACGCCCCGCCCCCCGCAACGCCCa                                                                           >  1:1026575/1‑239 (MQ=255)                                                                                                                                                       cggggttttCCCCGGCCCGCCAAAGGGCGCAAAATACTGCCCCCGGCGCGGCGGCCGGGAGGTTGGGTTAAAACCCAAGTTTGTCTCC‑CCAG‑AAAAC‑CCCCCTTCGCCCGGTAACAATATGTAGCCCAGGCCAAGCACCGGCACCTCGTCTTCCCGCGCATTGCGTTTTAGCAACGCGCTGCCGATGTGGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAgg                                                                         <  2:394753/231‑1 (MQ=255)                                                                                                                                                          cGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGCCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGGATAACGCTGCTCTTTCCACAAACTCCCCGCCTCCAGCAACGACCAGGcc                                                                       >  1:647458/1‑238 (MQ=255)                                                                                                                                                             tggccccGGCGCGCAAAGGTGGGGGAAATGCTGCCCCGCGGTCGGGGGCGGGGAGCTGGTGTTAAAACCCCAGCGTTTAACC‑ACAG‑CAAAA‑CACACTTTGCCCGGTAACAACATGGAACCCCGCCCAGGCCACGCCACCACTTCCCCCCGCGGGATACGTTTTAGCAACGCGCTCCCGATGTCGGTATAACGCTGCTCTTTCCCCAAACGCCCCGCCTCCAGCAACGACCAGGCCATc                                                                    <  2:999700/232‑1 (MQ=255)                                                                                                                                                             ggTCCAGGGCGCGCAAAGGCGGGTGAAATCCTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTACCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATc                                                                    <  1:758006/238‑1 (MQ=255)                                                                                                                                                                    ggCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCACTGCCGATTTCGGTATAACGCTGCTCTTTCCACAAACCCCCCCCCCCCAGCAACGCCCAGGCCCCCCCGACAta                                                            >  2:700807/1‑238 (MQ=255)                                                                                                                                                                         gcCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACAtaccac                                                        >  1:156798/1‑233 (MQ=255)                                                                                                                                                                         gcCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGAGGTCGGTATAACGCTGCCCTTTCCACAAACCCCCCGCCCCCAACAAACACCAAGCCATCCaaaaacaaaaaa                                                       >  2:452906/1‑228 (MQ=255)                                                                                                                                                                                       aatccgcggcCCGCGTGGGGGGCAGGTAGCGGGGGTTAAAACGCCAGTTGTTATCC‑TCAG‑CAAAC‑CCCCCTTTGCCCGGTAACAACAGGGAACCCAGCGCAGGCACCGTCCCCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTg                                         <  2:53953/231‑1 (MQ=255)                                                                                                                                                                                       aaTACGGCACCAGCGGCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTTTTATCC‑TCAG‑CAAAC‑CCCACTTTGCCCGGAAACAACATGGACCCCAGCCCAGGCCCCGTCCCCACTTCCTCCCGCGCGATACGTTTTAGCAACGCCCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTg                                         <  1:380774/239‑1 (MQ=255)                                                                                                                                                                                            tGCGCCAGCGGCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTAACC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGACCCCAGCCCAGGCACCGTAACCACTGCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTc                                    <  1:784065/239‑1 (MQ=255)                                                                                                                                                                                             gcgcCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACACCACCACCGGAGGCCGAATTCCTGTc                                    >  1:622850/1‑238 (MQ=255)                                                                                                                                                                                             gcgcCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACCCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGCCCAGGCCATCCACACATCACCACCCGACGCCGAAtacctctc                                    >  2:236016/1‑231 (MQ=255)                                                                                                                                                                                                    ggccggcggCTGGTTGCGGGGGTTAAACCGCCA‑‑TGTTTTCCACCAA‑CAAAACCCCACTTACCCGGTTAAGAACGAGGAACCCAGCCCGGGAATCGTCACCCATTCCCCCCGCGCGATACGTTTTAGCAACGCCCTTCCGGTGGCGGTATAACGCTGCTCTTCCCCCACACCTCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACtt                            <  2:629100/236‑1 (MQ=255)                                                                                                                                                                                                          cggcAGGTAGCTGGGGTTAAAACCCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCGGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCAc                       <  1:452906/238‑1 (MQ=255)                                                                                                                                                                                                               ggTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTa                    <  1:50620/236‑1 (MQ=255)                                                                                                                                                                                                                    cTGGGGTTAAACCGCCAGCTGTTATCC‑CCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCCCCGTCCCCACTTCCCCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTtct             <  1:359957/238‑1 (MQ=255)                                                                                                                                                                                                                      ggggGTAAAAACCCAGCTGTTCACC‑CCAC‑CAAAA‑CCCACTTTGCCCGGTACAAACCTGGACCCAAGCCCAGGCACCGTCACCACTTCCCCCCGCGGGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCCCAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCtctt          <  2:21731/239‑1 (MQ=255)                                                                                                                                                                                                                         gTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCCCCATCGGAGGCCGAATCTCGGTCCAGCACTTCCCACTTCCTGTTCtctt          >  2:50620/1‑236 (MQ=255)                                                                                                                                                                                                                          taaaaaCCCCAGTTGTTACCC‑CCCG‑CAAAC‑CCCCCCTGGCCCGGTAACAAAAGGGACCCCAGCCCAGGCCCCGTCACCACTCCCCCCCGCGCGATCCGTTTTACCAACGCGCTGCCGATGCCGGTATAACGCTGCTCTTCCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCtctttct       <  2:1059178/236‑1 (MQ=255)                                                                                                                                                                                                                              aaCCCCGGCTTTTACCC‑CCAG‑CAACA‑CCCGTTTTCCCCGGGAACAACATGGAACCACCCCCGGGCCCCGTCACCCCTTCCCCCCGCGCGAGACGTGTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCCCTTTCCCCAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGccc  <  2:622850/239‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                ACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGCCTCCAGCAACGACCAGGCCATCCAGACATCACCATCGGAGGCCGAATTGCTGTCCAGCACTTCCCACTTACTGTTCTCTTTCTTGCCC  >  NC_000913/3689509‑3689967

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: ‑