BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,177,337	G→T	5.4%	intergenic (‑133/+175)	gatY ← / ← fbaB	D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,177,337	0	G	T	5.4%	588.0 / 17.0	166	intergenic (‑133/+175)	gatY/fbaB	D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand): ref base G (77/80); new base T (6/3); total (83/83)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

Only 100 of 167 total aligned reads displayed.

CTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCG  >  NC_000913/2177108‑2177557

cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                                                                                                                                                   <  4:306595/273‑1 (MQ=255)

                                                                                 tCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCg                                                                                                                                                                                                                              >  5:1374304/1‑149 (MQ=255)

                                                                                        cACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACGAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGa                                                                                                                                                                                                                       <  5:793895/149‑1 (MQ=255)

                                                                                           gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGc                                                                                                                                                                                                                    <  8:545331/149‑1 (MQ=255)

                                                                                           gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGc                                                                                                                                                                                                                    >  5:1207343/1‑149 (MQ=255)

                                                                                           gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGc                                                                                                                                                                                                                    >  8:208790/1‑149 (MQ=255)

                                                                                               aTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                                                                                                                                                <  8:13705/149‑1 (MQ=255)

                                                                                                ttttCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGttt                                                                                                                                                                                                               <  7:848528/149‑3 (MQ=255)

                                                                                                 tttCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGt                                                                                >  4:168708/1‑275 (MQ=255)

                                                                                                    cATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaaa                                                                                                                                                                                                           <  7:208790/149‑1 (MQ=255)

                                                                                                    cATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaaa                                                                                                                                                                                                           <  7:202460/149‑1 (MQ=255)

                                                                                                      tatCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGt                                                                                                                                                                                                         >  5:1200893/1‑149 (MQ=255)

                                                                                                        tCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgc                                                                                                                                                                                                       >  5:1090065/1‑149 (MQ=255)

                                                                                                          cTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                                                                                                                                                                                                     <  8:244490/149‑1 (MQ=255)

                                                                                                                  tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt                                                                                                                                                                                             <  5:788648/149‑1 (MQ=255)

                                                                                                                     tttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                                                                                                                                                          >  7:880842/1‑149 (MQ=255)

                                                                                                                     tttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                                                                                                                                                          >  8:1192680/1‑149 (MQ=255)

                                                                                                                          aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAACAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTAATGCCgg                                                                                                                                                                                      <  8:1235843/148‑1 (MQ=255)

                                                                                                                          aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg                                                                                                                                                                                      >  7:1235843/1‑148 (MQ=255)

                                                                                                                          aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa                                                                                                                                                                                     >  6:301464/1‑149 (MQ=255)

                                                                                                                           tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                                                                                                                                                    >  5:150068/1‑149 (MQ=255)

                                                                                                                           tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                                                                                                                                                    <  6:3374/149‑1 (MQ=255)

                                                                                                                            ttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAg                                                     <  2:168708/275‑1 (MQ=255)

                                                                                                                               aaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                                                                                                                                                >  8:1032765/1‑149 (MQ=255)

                                                                                                                                aaaTATAATTAAATTTTTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgt                                                                                                                                                                               <  6:1109135/149‑1 (MQ=255)

                                                                                                                                  atatAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtgt                                                                                                                                                                             >  7:497416/1‑149 (MQ=255)

                                                                                                                                     taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGGGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                                                                                                     >  4:382540/1‑154 (MQ=255)

                                                                                                                                     taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                                                                                                                                                          <  8:497416/149‑1 (MQ=255)

                                                                                                                                     taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                                                                                                                                                          <  8:650061/149‑1 (MQ=255)

                                                                                                                                        tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGc                                                                                                                                                                                                                    <  8:52316/104‑1 (MQ=255)

                                                                                                                                        tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACAt                                                                                                                                                                       <  7:298900/149‑1 (MQ=255)

                                                                                                                                        tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGCGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGc                                                                                                                                                                                                                    >  7:52316/1‑104 (MQ=255)

                                                                                                                                         gAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGGGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTgg                                                                                                                                                                 <  2:382541/154‑1 (MQ=255)

                                                                                                                                          aaaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTACGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                                                                                                     >  8:949479/1‑149 (MQ=255)

                                                                                                                                             ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                                                                                                                                                  <  6:1374304/149‑1 (MQ=255)

                                                                                                                                                 ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAg                                                                                                                                                              >  7:8653/1‑149 (MQ=255)

                                                                                                                                                    tttttAATTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGc                                                                                                                                                           <  8:519504/149‑1 (MQ=255)

                                                                                                                                                       ttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTAt                                                                                                                                                        <  6:1455669/149‑1 (MQ=255)

                                                                                                                                                       ttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTAt                                                                                                                                                        <  7:1063073/149‑1 (MQ=255)

                                                                                                                                                         aaaTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGt                                                                                                                                                      <  5:301464/149‑1 (MQ=255)

                                                                                                                                                          aaTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                                                                                                                                                     <  8:89448/149‑1 (MQ=255)

                                                                                                                                                            tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                                                                                                                                                          <  6:467086/126‑1 (MQ=255)

                                                                                                                                                            tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                                                                                                                                                          >  5:467086/1‑126 (MQ=255)

                                                                                                                                                            tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAgg                                                                                                                                                   <  6:1200893/149‑1 (MQ=255)

                                                                                                                                                                 agaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                                                                                                                                                          >  7:481105/1‑121 (MQ=255)

                                                                                                                                                                 agaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGa                                                                                                                                              >  7:11692/1‑149 (MQ=255)

                                                                                                                                                                    tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATaa                                                                                                                                           <  7:420005/149‑1 (MQ=255)

                                                                                                                                                                      aCGATCACAAAAACGACAATATGAAAATTAGTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGa                                                                                                                                         <  5:136026/149‑1 (MQ=255)

                                                                                                                                                                       cGATCACAAAAACGACAATATGAAAAGTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                                                                                                                                        <  7:768284/149‑1 (MQ=255)

                                                                                                                                                                        gaccacaAAAACGACAATATGAAAACGATAAGACGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGAGCGATGTACGCATGATAAGACg                                                                                                                                       >  6:133776/4‑149 (MQ=255)

                                                                                                                                                                        gATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                                                                                                                                                >  6:1225715/1‑76 (MQ=255)

                                                                                                                                                                        gATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                                                                                                                                         >  7:142238/1‑115 (MQ=255)

                                                                                                                                                                        gATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                                                                                                                                       >  5:310915/1‑149 (MQ=255)

                                                                                                                                                                        gATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTACCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGTTGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATCAGACg                                                                                                                                       >  6:901074/1‑149 (MQ=255)

                                                                                                                                                                         aTCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg                                                                                                                                                                                 >  7:669814/1‑106 (MQ=255)

                                                                                                                                                                         aTCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGc                                                                                                                                      <  8:786129/149‑1 (MQ=255)

                                                                                                                                                                            acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTc                                                                                                                                   >  7:1239176/1‑149 (MQ=255)

                                                                                                                                                                             caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGCTAAGCCGTAGCAAAGTGCGCGTTTTTTTATGCCGGATGCGGTGTAAACATTCGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCa                                                                                                                                  <  8:1061595/149‑1 (MQ=255)

                                                                                                                                                                                aaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCTTATCAGGCGATAGTAATTTTGCTATCGAGTTAAACGTCCTGCACGGCg  >  2:157312/1‑274 (MQ=255)

                                                                                                                                                                                  aCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTc                                                                                                                             <  5:22356/149‑1 (MQ=255)

                                                                                                                                                                                    gACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGc                                                                                                                           <  6:1325469/149‑1 (MQ=255)

                                                                                                                                                                                        atatGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCa                                                                                                                       >  7:1439621/1‑149 (MQ=255)

                                                                                                                                                                                         tatGAAAATTATTCGAGGAGTGAACGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACTCTTCAGCGTCGCATCAg                                                                                                                      >  5:1018474/1‑149 (MQ=255)

                                                                                                                                                                                         tatGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                                                                                                                      <  7:666123/149‑1 (MQ=255)

                                                                                                                                                                                             aaaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGGTTCAGCGTCGCATCAGGCAt                                                                                                                  <  5:1180371/149‑1 (MQ=255)

                                                                                                                                                                                             aaaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCAt                                                                                                                  >  8:59218/1‑149 (MQ=255)

                                                                                                                                                                                              aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTGAAGCGTct                                                                                                                            <  5:478583/138‑2 (MQ=255)

                                                                                                                                                                                              aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTGAAGCGTct                                                                                                                            >  6:478583/1‑137 (MQ=255)

                                                                                                                                                                                               aaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg                                                                                                                >  6:344959/1‑149 (MQ=255)

                                                                                                                                                                                                    ttCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATa                                                                                                                 <  5:369577/143‑1 (MQ=255)

                                                                                                                                                                                                       gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAgg                                                                                                               <  2:47015/142‑1 (MQ=255)

                                                                                                                                                                                                       gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                                                                                                        >  8:139371/1‑149 (MQ=255)

                                                                                                                                                                                                       gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGGGTaa                                                                                                                                                                           >  6:295623/1‑82 (MQ=255)

                                                                                                                                                                                                         ggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                                                                                                      <  8:1310293/149‑1 (MQ=255)

                                                                                                                                                                                                            gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATgta                                                                                                   <  6:1090065/149‑3 (MQ=255)

                                                                                                                                                                                                            gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGAt                   >  4:333410/1‑229 (MQ=255)

                                                                                                                                                                                                              gAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGc                                                                                                 <  7:1096088/149‑1 (MQ=255)

                                                                                                                                                                                                               aaaGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCg                                                                                                >  5:734345/1‑149 (MQ=255)

                                                                                                                                                                                                                aaGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAAc               <  2:333410/229‑1 (MQ=255)

                                                                                                                                                                                                                    cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc                                                                                           <  8:1164011/149‑1 (MQ=255)

                                                                                                                                                                                                                     aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg                                                                                          >  8:742900/1‑149 (MQ=255)

                                                                                                                                                                                                                     aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGACGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg                                                                                          >  8:875623/1‑149 (MQ=255)

                                                                                                                                                                                                                      aaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                                                                                                      >  8:227222/1‑136 (MQ=255)

                                                                                                                                                                                                                      aaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                                                                                                      <  7:227222/136‑1 (MQ=255)

                                                                                                                                                                                                                      aaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGc                                                                                         <  8:1239176/149‑1 (MQ=255)

                                                                                                                                                                                                                       aaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                                                                                                      >  7:98273/1‑135 (MQ=255)

                                                                                                                                                                                                                       aaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGcc                                                                                        >  5:142810/1‑149 (MQ=255)

                                                                                                                                                                                                                          taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc                                                                                           >  8:159936/2‑143 (MQ=39)

                                                                                                                                                                                                                          taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc                                                                                           <  7:159936/142‑1 (MQ=39)

                                                                                                                                                                                                                          taCGGCTTCCCTATCGGGAGGCCATAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTa                                                                                     >  8:91857/2‑149 (MQ=39)

                                                                                                                                                                                                                          aaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTa                                                                                     <  8:434502/149‑1 (MQ=255)

                                                                                                                                                                                                                            cGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATc                                                                                   >  5:1031677/1‑149 (MQ=255)

                                                                                                                                                                                                                             ggCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATcc                                                                                  >  6:1029544/1‑149 (MQ=39)

                                                                                                                                                                                                                              gCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCAt                                                                                                                                                <  7:1326949/86‑1 (MQ=255)

                                                                                                                                                                                                                                cTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTc                                                                               <  8:777159/149‑1 (MQ=255)

                                                                                                                                                                                                                                cTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACCTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTc                                                                               >  8:1108970/1‑149 (MQ=255)

                                                                                                                                                                                                                                 tCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTcc                                                                              <  8:880842/149‑1 (MQ=255)

                                                                                                                                                                                                                                 tCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTcc                                                                              >  5:1481751/1‑149 (MQ=255)

                                                                                                                                                                                                                                   ccGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTa                                                                            >  8:370169/1‑149 (MQ=255)

                                                                                                                                                                                                                                   ccGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTa                                                                            >  7:271800/1‑149 (MQ=255)

CTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCG  >  NC_000913/2177108‑2177557

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`