BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,689,747	T→C	18.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,747	0	T	C	18.4%	138.6 / 3.9	62	E172G (GAG→GGG)	bcsZ	endo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand): ref base T (20/29); new base C (5/7); total (25/37)
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/3689510‑3689968

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)

                         cgtAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCGTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTaa                                                                                                                                                                                                      <  1:1027978/238‑1 (MQ=255)

                          gttgcttattaTCATTTTTTTTTCGGCTTTGAACTGCCAGCCTGTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCG‑CAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATTCTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCC‑TCAG‑CAAAA‑CCCACTTTGCCCGGTaaca                                                                                                                                                                                                    <  2:929367/229‑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)

                                                     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)

                                                                    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)

                                                                                                   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)

                                                                                                                                        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/3689510‑3689968

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