BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	940,751	T→C	7.7%	intergenic (+31/‑208)	serS → / → dmsA	seryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	940,751	0	T	C	7.7%	620.9 / 29.0	195	intergenic (+31/‑208)	serS/dmsA	seryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand): ref base T (91/89); new base C (9/6); total (100/95)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

Only 100 of 197 total aligned reads displayed.

CCTACCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGAGCAAGCTTGAGTGAGCCATTATGAAAACGAAAATCCCT  >  NC_000913/940482‑940976

ccTACCGTGGGATCCTTTCCTGCGCCAACGTTTGGGATTTCCAGGAACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCggg                                                                                                                                                                                                                              <  4:122893/275‑1 (MQ=255)

 cTACCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGAATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCATGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGGAGTACGAGAAGTTCTGCGTCAGTATATGAACGGACTGGAATATATTGGATAATACACAATTTTTCTGAATCTAAAAAACGCCTGCggg                                                                                                                                                                                                                              >  1:67848/1‑274 (MQ=255)

                 tCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTAAAGTACCAGACGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGCCGCTTTTTTTGTCTcct                                                                                                                                                                                                             >  2:49273/1‑274 (MQ=255)

                          aaCGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCACTTTTTCTGAATCTAAAAAGCGACTGCGGGCGCTTTTTTTGTCCCCCttttaata                                                                                                                                                                                                     >  4:80111/1‑269 (MQ=255)

                           aCGTGTGGGTTTTCCATGCACGTCGGATGAAGGCAAGTTGCCGAAGAAAGCCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGGACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCg                                                                                                                                                                                                   <  4:284231/275‑1 (MQ=255)

                                         cAGGAACGTCGTAGGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                                                                                                     <  1:50746/275‑1 (MQ=255)

                                           ggCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                                                                                                                                                   <  2:260710/275‑1 (MQ=255)

                                                                                 aCAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATCACCAATTTTTCTGAATCTAAAAAGCGCCTGGGGGCGCTTTTTTTGTATCCCTTTGATACCGAACAATAATTACTCCTCAATTACAAGTAATACAACTTTTGAGTGaaaaacaccc                                                                                                                                             >  4:205572/1‑269 (MQ=255)

                                                                                                            tGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAANTTTTCTGAATCTAAAAAGAGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGCAAAATAATTACTCCTCACTTACACGTACTCCTACTTTCGAGTGaaaaaccncn                                                                                                                                             >  3:63802/1‑242 (MQ=255)

                                                                                                              aaagttgttgggttggCGGTTGGACGTACGCTGGTTGCAGTAAAGGAAAACTATCAGCAGGCGGATTGTCGTATTGAAGTACCAGAAGTGCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTCATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTGTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATGCGAt                                                                                                                 <  4:133060/261‑1 (MQ=255)

                                                                                                                 ggTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                           <  2:203527/213‑1 (MQ=255)

                                                                                                                   ttCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAtttt                                                                                                                              >  4:19805/1‑256 (MQ=255)

                                                                                                                           tggCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGc                                                                                                                                                                                                                                 >  5:673730/1‑149 (MQ=255)

                                                                                                                           tggCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGc                                                                                                                                                                                                                                 <  5:51136/149‑1 (MQ=255)

                                                                                                                                 ttGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgcgc                                                                                                                                                                                                                           >  8:612837/1‑149 (MQ=255)

                                                                                                                                 ttGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                           <  2:358377/197‑1 (MQ=255)

                                                                                                                                      cGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttt                                                                                                                                                                                                                      <  6:1205638/149‑1 (MQ=255)

                                                                                                                                       gTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                                                                                                                                                     <  7:1277291/149‑2 (MQ=255)

                                                                                                                                         aCGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                                                                                                                                                   <  5:435810/149‑1 (MQ=255)

                                                                                                                                         aCGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttaga                                                                                                                                                                                                                   <  7:367646/149‑4 (MQ=255)

                                                                                                                                          cGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGt                                                                                                                                                                                                                  <  8:933749/149‑1 (MQ=255)

                                                                                                                                                ttGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCt                                                                                                                                                                                                            <  7:919664/149‑1 (MQ=255)

                                                                                                                                                 tGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCtt                                                                                                                                                                                                           <  5:324847/149‑1 (MQ=255)

                                                                                                                                                  gCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttt                                                                                                                                                                                                          <  5:2296/149‑1 (MQ=255)

                                                                                                                                                      tAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATa                                                                                                                                                                                                      >  5:846962/1‑149 (MQ=255)

                                                                                                                                                      tAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATa                                                                                                                                                                                                      >  8:866195/1‑149 (MQ=255)

                                                                                                                                                         tGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttt                                                                                                                                                                                                                      <  6:235883/130‑1 (MQ=255)

                                                                                                                                                           gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                 >  7:177458/1‑149 (MQ=255)

                                                                                                                                                           gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                 <  7:1090081/149‑1 (MQ=255)

                                                                                                                                                           gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                 >  8:288968/1‑149 (MQ=255)

                                                                                                                                                              aaCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGATTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaa                                                                                                                                                                                              >  5:575283/1‑149 (MQ=255)

                                                                                                                                                               accaTCAGCCTGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaat                                                                                                                                                                                             >  5:1315900/4‑149 (MQ=255)

                                                                                                                                                                 tATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                                                                                                                                                           >  7:890257/1‑149 (MQ=255)

                                                                                                                                                                 caTCAGCCTGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                                                                                                                                                           >  8:155932/2‑149 (MQ=255)

                                                                                                                                                                  aTCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataat                                                                                                                                                                                          <  8:1177049/149‑1 (MQ=255)

                                                                                                                                                                   tCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttaga                                                                                                                                                                                                                   <  8:224493/123‑4 (MQ=255)

                                                                                                                                                                     agcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAc                                                                                                                                                                                       >  8:63919/1‑149 (MQ=255)

                                                                                                                                                                     agcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttaga                                                                                                                                                                                                                   >  6:407829/1‑118 (MQ=255)

                                                                                                                                                                     agcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttaga                                                                                                                                                                                                                   <  5:407829/121‑4 (MQ=255)

                                                                                                                                                                     agcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCATGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAc                                                                                                                                                                                       <  8:890257/149‑1 (MQ=255)

                                                                                                                                                                      gcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAct                                                                                                                                                                                      >  6:776733/1‑149 (MQ=255)

                                                                                                                                                                      gcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAct                                                                                                                                                                                      >  6:186794/1‑149 (MQ=255)

                                                                                                                                                                       cagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                                                                                                     <  5:909798/149‑1 (MQ=255)

                                                                                                                                                                              aTGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                           <  4:31784/152‑1 (MQ=255)

                                                                                                                                                                                 gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                           <  5:658496/149‑1 (MQ=255)

                                                                                                                                                                                  tCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAc                                                                                                                                                                                       >  5:1218908/1‑136 (MQ=255)

                                                                                                                                                                                  tCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTacac                                                                                                                                                                          >  6:69853/1‑149 (MQ=255)

                                                                                                                                                                                    gTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTATTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGACATAATTACACAATACGTTTTTTTACTTGACTCAATCGt                                          >  2:81787/1‑275 (MQ=255)

                                                                                                                                                                                      aTTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTaa                                                                                                                                                                      <  6:1286314/149‑1 (MQ=255)

                                                                                                                                                                                      aTTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTaa                                                                                                                                                                      >  6:1335332/1‑149 (MQ=255)

                                                                                                                                                                                        tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTAGACGTAAta                                                                                                                                                                    <  7:1345434/149‑1 (MQ=255)

                                                                                                                                                                                         gAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                           >  7:642570/1‑141 (MQ=255)

                                                                                                                                                                                         gAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                           <  8:642570/141‑1 (MQ=255)

                                                                                                                                                                                               ccAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACttt                                                                                                                                                             <  6:512362/149‑1 (MQ=255)

                                                                                                                                                                                                cAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                                                                                                                                                           <  8:974099/118‑1 (MQ=255)

                                                                                                                                                                                                cAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                                                                                            >  6:1015271/1‑149 (MQ=255)

                                                                                                                                                                                                 aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCg                                                                                                                                                           >  7:700686/1‑149 (MQ=255)

                                                                                                                                                                                                 aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCg                                                                                                                                                           >  5:1471463/1‑149 (MQ=255)

                                                                                                                                                                                                        cTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcc                                                                                                                                                                                    >  2:128633/1‑117 (MQ=255)

                                                                                                                                                                                                         tGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaaa                                                                                                                                                   >  5:257392/1‑149 (MQ=255)

                                                                                                                                                                                                                      gAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTaa                                                                                                                                                                      <  8:1383565/117‑1 (MQ=255)

                                                                                                                                                                                                                      gAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATctct                                                                                                                                      <  5:1335332/149‑1 (MQ=255)

                                                                                                                                                                                                                          ggACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCAc                                                                                                                                                                                >  6:628466/1‑103 (MQ=255)

                                                                                                                                                                                                                              tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                                                                                              <  8:532361/133‑1 (MQ=255)

                                                                                                                                                                                                                              tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAtttt                                                                                                                              <  6:17775/149‑1 (MQ=255)

                                                                                                                                                                                                                                gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCa                                                                                                                            >  7:1258666/1‑149 (MQ=255)

                                                                                                                                                                                                                                gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCa                                                                                                                            >  6:29550/1‑149 (MQ=255)

                                                                                                                                                                                                                                    atatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAatt                                                                                                                        <  6:919681/149‑1 (MQ=255)

                                                                                                                                                                                                                                       ttGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAattatt                                                                                                                     >  8:1221768/1‑149 (MQ=255)

                                                                                                                                                                                                                                       ttGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACAAGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAattatt                                                                                                                     >  8:1420351/1‑149 (MQ=255)

                                                                                                                                                                                                                                          gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGa                                                                                                                  <  5:10258/149‑1 (MQ=255)

                                                                                                                                                                                                                                          gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGa                                                                                                                  >  6:556371/1‑149 (MQ=255)

                                                                                                                                                                                                                                          gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGa                                                                                                                  <  7:866195/149‑1 (MQ=255)

                                                                                                                                                                                                                                          gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGAt                                                                                                                 <  4:362338/150‑1 (MQ=255)

                                                                                                                                                                                                                                          gCTAATACCCAATTGTTCTGAAACTAAAAAGCGACTGCGGGCGCTTTTTTTGTCTCAATTTGATACCGAACAATAATTGCTTTGAAAATACACGTAATACTTCTTTCGAGTGCAATTCTACCTATGGCTTTCATTTTCAACTTCTCCGa                                                                                                                  <  6:1108710/149‑1 (MQ=255)

                                                                                                                                                                                                                                            tAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtactac                                                                                                                                                                >  5:946683/1‑101 (MQ=255)

                                                                                                                                                                                                                                             aaTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGt                                                                                                               >  7:1402945/1‑149 (MQ=255)

                                                                                                                                                                                                                                               tACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                                                                                            <  2:303284/102‑1 (MQ=255)

                                                                                                                                                                                                                                                  ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                                                                                                                                                   <  5:1037095/76‑1 (MQ=255)

                                                                                                                                                                                                                                                  ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTTGATGTATACa                                                                                                          >  5:1392542/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                                                                                                          >  8:1018183/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   cAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                                                                                                          <  2:216290/148‑1 (MQ=255)

                                                                                                                                                                                                                                                    aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAg                                                                                                        >  6:431976/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAg                                                                                                        >  8:1089475/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       ttttCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCt                                                                                                     <  7:1142174/149‑1 (MQ=255)

                                                                                                                                                                                                                                                           cTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTCCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtatat                                                                                                 >  7:1010750/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtatata                                                                                                >  7:956774/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             gAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAg                                                                                               <  7:1144871/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGa                                                                                            >  6:989980/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  tAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAc                                                                                          >  8:1342003/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  tAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAc                                                                                          >  5:124732/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  tAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAc                                                                                          >  6:1075227/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   aaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAat                                                                                                                         <  5:539046/117‑1 (MQ=255)

                                                                                                                                                                                                                                                                   aaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACt                                                                                         >  6:1027336/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    aaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCCAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTg                                                                                        >  5:1416511/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        gcgcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt                                                                                                                                    <  7:232515/101‑1 (MQ=255)

                                                                                                                                                                                                                                                                          gcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAg                                                                                  >  5:764462/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                           ccTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGa                                                                                 >  7:298342/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                            cTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGaa                                                                                <  6:37441/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                             tGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAAAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGaaa                                                                               <  6:575283/149‑1 (MQ=255)

CCTACCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGAGCAAGCTTGAGTGAGCCATTATGAAAACGAAAATCCCT  >  NC_000913/940482‑940976

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