BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	314,258	G→A	7.6%	intergenic (‑453/‑99)	ykgP ← / → eaeH	pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	314,258	0	G	A	7.6%	489.8 / 21.9	157	intergenic (‑453/‑99)	ykgP/eaeH	pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand): ref base G (77/68); new base A (6/6); total (83/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

Only 100 of 159 total aligned reads displayed.

CTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGAATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCAT  >  NC_000913/314006‑314526

cTATCTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCCCTATGAATAGTGATTTATTCACTGTCCGCAGAGTATGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAATGCTAAATTTTCCCCCATAATTAAAAATATAATCCCGCGCAAAACCACCTGATGAGTAGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTAACTTAAGGTATTTCCTTTTCTGGGGAAAGGGAACTAGGGAAATgcc                                                                                                                                                                                                                                                        >  4:292178/1‑275 (MQ=255)

                                                                                       aaGAAATATAACTTAGGTATCTATTTAATGACTTTCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCTCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCATTTTCTTACATAAAGGTATTTCCTTTTCTGCGGACAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTTTTCAGGAAATTATCTCATTGTCACATta                                                                                                                                                                 >  2:711364/1‑275 (MQ=255)

                                                                                                         aTCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCa                                                                                                                                                                                                                               <  4:398382/195‑1 (MQ=255)

                                                                                                             aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                                                                                                                                         <  7:257853/149‑1 (MQ=255)

                                                                                                             aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaa                                                                                                                                                                                                                                                                          <  6:805871/148‑1 (MQ=255)

                                                                                                             aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGaaa                                                                                                                                                                                                                                                                          <  5:448143/148‑1 (MQ=255)

                                                                                                                tAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTTCATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                                                                                                                                                                                                      <  7:592366/149‑1 (MQ=255)

                                                                                                                tAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                                                                                                                                                                                                      <  7:6546/149‑1 (MQ=255)

                                                                                                                       ttGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                                                                                                                                                >  7:225329/1‑148 (MQ=255)

                                                                                                                            caAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACCTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATc                                                                                                                                                                                                                                                          >  8:286195/1‑149 (MQ=255)

                                                                                                                                aaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcccc                                                                                                                                                                                                                                                       >  5:876359/1‑148 (MQ=255)

                                                                                                                                aaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCg                                                                                                                                                                                                                                                      >  8:805325/1‑149 (MQ=255)

                                                                                                                                    tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACTTAAAGGTATTTCCTTTTCCGCAGAAAAGGaa                                                                                                                                                                                                                                                                     <  5:447244/130‑1 (MQ=255)

                                                                                                                                    tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaa                                                                                                                                                                                                                                                                     >  6:447244/1‑130 (MQ=255)

                                                                                                                                     aaaTTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGttt                                                                                                                                                                                                                                                  >  5:303926/1‑148 (MQ=255)

                                                                                                                                           tCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGAc                                                                                                                                                                                                                                   >  2:661106/1‑157 (MQ=255)

                                                                                                                                             ccccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACaa                                                                                                                                                                                                                                         >  8:121575/1‑149 (MQ=255)

                                                                                                                                               ccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTAACAAGc                                                                                                                                                                                                                                       >  8:644206/1‑149 (MQ=255)

                                                                                                                                                     ataaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACaa                                                                                                                                                                                                                                         >  8:346349/1‑141 (MQ=255)

                                                                                                                                                     ataaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACaa                                                                                                                                                                                                                                         <  7:346349/141‑1 (MQ=255)

                                                                                                                                                         aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGAc                                                                                                                                                                                                                                           >  8:260163/1‑135 (MQ=255)

                                                                                                                                                         aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGAc                                                                                                                                                                                                                                           <  7:260163/135‑1 (MQ=255)

                                                                                                                                                         aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatt                                                                                                                                                                                                                             >  7:130897/1‑149 (MQ=255)

                                                                                                                                                         aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatt                                                                                                                                                                                                                             >  8:189348/1‑149 (MQ=255)

                                                                                                                                                         aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCACTGCGTTGCGTGGGCaaa                                                                                               >  1:179753/1‑275 (MQ=255)

                                                                                                                                                          gataacaaaccGCGCCCAACCACCTGAGTGGTGGCTATAGGCACGGGGTATATCAGGGGGCGAGACACTTACTTACAGAAATGTATTTCCTTTTCTGCGGAAAAGTAAATCGGTAAATCCCCTGTTTTTCTTACAAGCAGACGCAATTATTTGTGTCTGCCTATGTTCATTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCGGTTCTGGCCCGCTGCGTGGCGTGGGCaaa                                                                                               <  4:99415/265‑1 (MQ=255)

                                                                                                                                                               aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgt                                                                                                                                                                                                                       <  8:890577/149‑1 (MQ=255)

                                                                                                                                                                aTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                                                                                                                                                <  7:901242/107‑1 (MQ=255)

                                                                                                                                                                 tCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAAGCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCGCAATGTCACATTATAAAACAGTTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCGTCGTGGCTTGGGCAAATATctctt                                                                                       >  4:536292/1‑274 (MQ=255)

                                                                                                                                                                     gcgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgt                                                                                                                                                                                                                       <  2:726550/143‑1 (MQ=255)

                                                                                                                                                                     gcgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGc                                                                                                                                                                                                                 <  8:869052/149‑1 (MQ=255)

                                                                                                                                                                      cgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGc                                                                                                                                                                                                                                       >  4:802564/1‑126 (MQ=255)

                                                                                                                                                                      cgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTgcgtg                                                                                                          >  2:233505/1‑251 (MQ=255)

                                                                                                                                                                      cgcCCAACCACCTGATGAGTGGCTATAGGCACTGGACATATCAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTGTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAACTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATCCAGTTCGGGACCCATGCGTGCGCTGGCCAAattgccatgttca                                                                                   >  4:712025/1‑263 (MQ=255)

                                                                                                                                                                        cccAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  2:503119/1‑185 (MQ=255)

                                                                                                                                                                          cAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAg                                                                                                                                                                                                                                     <  2:24178/124‑1 (MQ=255)

                                                                                                                                                                          cAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGAc                                                                                                                                                                                                                                   <  2:802566/126‑1 (MQ=255)

                                                                                                                                                                          cAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTgcgtggcgt                                                                                                      <  4:233504/251‑1 (MQ=255)

                                                                                                                                                                            accaccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                                      <  4:503117/185‑1 (MQ=255)

                                                                                                                                                                               accTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGt                                                                                                                                                                                                       >  7:460582/1‑149 (MQ=255)

                                                                                                                                                                                ccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATTAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGtt                                                                                                                                                                                                      <  7:644206/149‑1 (MQ=255)

                                                                                                                                                                                   gatgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGAAAATAGCTCTGTTCAGGTTCTTTTTTCAc                                                                     >  1:97318/1‑275 (MQ=255)

                                                                                                                                                                                     tgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  3:187827/1‑172 (MQ=255)

                                                                                                                                                                                        gTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGtt                                                                                                                                                                                              <  8:460582/149‑1 (MQ=255)

                                                                                                                                                                                         tGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                                      <  1:187828/172‑1 (MQ=255)

                                                                                                                                                                                         tGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTCCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtc                                                                                                                                                                                             <  7:288656/149‑1 (MQ=255)

                                                                                                                                                                                           gCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCACGCTGCGTGGCGTGGGCAAATATCTCTGTTAAGGTTCTTTTTCCACTCGCTGta                                                             >  2:613722/1‑274 (MQ=255)

                                                                                                                                                                                            cTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatta                                                                                                                                                                                                                            >  7:690754/1‑115 (MQ=255)

                                                                                                                                                                                             tataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCAGCGGAAAAGGAAATCGGGTAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatca                                                                                                                                                                                         >  8:419129/1‑149 (MQ=255)

                                                                                                                                                                                                     cTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAAtt                                                                                                                                                                                  <  6:185666/148‑1 (MQ=255)

                                                                                                                                                                                                      tGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  4:395828/1‑155 (MQ=255)

                                                                                                                                                                                                        gATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtct                                                                                                                                                                              >  7:724977/1‑149 (MQ=255)

                                                                                                                                                                                                        gATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtct                                                                                                                                                                              >  7:169192/1‑149 (MQ=255)

                                                                                                                                                                                                          tatatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCa                                                                                                                                                                            >  7:86031/1‑149 (MQ=255)

                                                                                                                                                                                                          tatatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                                      <  2:395830/155‑1 (MQ=255)

                                                                                                                                                                                                           atatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCaa                                                                                                                                                                           >  8:257915/1‑149 (MQ=255)

                                                                                                                                                                                                              ttAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGt                                                                                                                                                                        >  7:629585/1‑149 (MQ=255)

                                                                                                                                                                                                               tAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGcc                                                                                                                                                                                                                >  6:806885/1‑108 (MQ=255)

                                                                                                                                                                                                               tAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGt                                                                                                                                                                        >  6:286320/1‑148 (MQ=255)

                                                                                                                                                                                                                 ggTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCgg                                                                                                                                                                                                                                                     >  8:635731/1‑69 (MQ=255)

                                                                                                                                                                                                                  gTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  3:19907/1‑143 (MQ=255)

                                                                                                                                                                                                                  gTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  3:157335/1‑143 (MQ=255)

                                                                                                                                                                                                                   tGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTcaca                                                                                                                                                                    >  5:510257/1‑148 (MQ=255)

                                                                                                                                                                                                                    ggCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACAtt                                                                                                                                                                  <  7:189348/149‑1 (MQ=255)

                                                                                                                                                                                                                      cGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                                      <  1:19907/143‑1 (MQ=255)

                                                                                                                                                                                                                       ggTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATtata                                                                                                                                                               <  7:875968/149‑1 (MQ=255)

                                                                                                                                                                                                                        gTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaa                                                                                                                                                              >  8:5272/1‑149 (MQ=255)

                                                                                                                                                                                                                        gTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTAtt                                                                                                                              <  1:113714/181‑1 (MQ=255)

                                                                                                                                                                                                                         tGCACTTTCTTACATAAAGGTATTTCCTTTTTTGCGTAAAAGGAGTCCGGGAATTGCCCGGTTTTTCTTACACGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGACACGTgcg                               <  4:613720/275‑1 (MQ=255)

                                                                                                                                                                                                                               ttCTTACATAAAGGTATTCCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGt                                                                                                                                                       <  8:156296/149‑1 (MQ=255)

                                                                                                                                                                                                                                    aCATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                                                                                                                                           >  5:831680/1‑140 (MQ=255)

                                                                                                                                                                                                                                     cATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATaaa                                                                                                                                                 <  8:849467/149‑1 (MQ=255)

                                                                                                                                                                                                                                        aaaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  4:670728/1‑121 (MQ=255)

                                                                                                                                                                                                                                        aaaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAaacaa                                                                                                                                              <  7:419129/149‑1 (MQ=255)

                                                                                                                                                                                                                                          aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  4:732826/1‑119 (MQ=255)

                                                                                                                                                                                                                                          aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  3:201336/1‑119 (MQ=255)

                                                                                                                                                                                                                                          aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATtata                                                                                                                                                               >  2:784319/1‑130 (MQ=255)

                                                                                                                                                                                                                                          aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGt                                                                                                                          >  2:369439/1‑167 (MQ=255)

                                                                                                                                                                                                                                           ccgaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCa                                                                                                                                           >  7:134457/4‑149 (MQ=255)

                                                                                                                                                                                                                                              aTTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                                      <  1:201337/119‑1 (MQ=255)

                                                                                                                                                                                                                                              aTTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTCTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACg                                                                                                                                         <  5:275368/148‑1 (MQ=255)

                                                                                                                                                                                                                                               tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  1:123891/1‑114 (MQ=255)

                                                                                                                                                                                                                                               tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                          >  2:416442/1‑114 (MQ=255)

                                                                                                                                                                                                                                               tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCACGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCCCATTATAAAANAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCACGCTGCTTGGCGTGGGCAAATATCTCTGTNCAGGTTCTTTTTCCATTCGTTGTCACCTTTTCTCCAGTAATTGCGNCNCTTTCCCACCCTGCGGTTCTGCCCCGGt         >  3:72089/1‑275 (MQ=255)

                                                                                                                                                                                                                                                 tCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAttt                                                                                                                                     >  7:8310/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 tCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAtt                                                                                                                                      >  6:883303/1‑148 (MQ=255)

                                                                                                                                                                                                                                                 tCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAtt                                                                                                                                      >  6:321871/1‑148 (MQ=255)

                                                                                                                                                                                                                                                   cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGt                                                                                                                                                                                               <  7:392810/89‑1 (MQ=255)

                                                                                                                                                                                                                                                   cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                                      <  3:123890/114‑1 (MQ=255)

                                                                                                                                                                                                                                                   cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                                                                                                                                           <  8:708995/125‑1 (MQ=255)

                                                                                                                                                                                                                                                   cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                                                                                                                                           >  7:708995/1‑125 (MQ=255)

                                                                                                                                                                                                                                                    ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAg                                                                                                                                                         <  7:885949/126‑1 (MQ=255)

                                                                                                                                                                                                                                                    ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCg                                                                                                                                   <  6:267394/148‑1 (MQ=255)

                                                                                                                                                                                                                                                     tttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGtt                                                                                                                                 <  7:805325/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     tttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCt                                                         >  2:495441/1‑221 (MQ=255)

                                                                                                                                                                                                                                                      ttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCAt  <  2:611041/275‑1 (MQ=255)

                                                                                                                                                                                                                                                         tGCGGAAAAGGAAATCGTGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTAc                                                     <  4:495439/221‑1 (MQ=255)

                                                                                                                                                                                                                                                         tGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTc                                                                                                                             >  8:730775/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          gCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                                      <  4:613977/107‑1 (MQ=255)

                                                                                                                                                                                                                                                            ggAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAtt                                                                                                                                      >  8:379889/1‑137 (MQ=255)

CTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGAATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCAT  >  NC_000913/314006‑314526

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