BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	940,754	G→A	6.0%	intergenic (+34/‑205)	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,754	0	G	A	6.0%	979.0 / 35.5	302	intergenic (+34/‑205)	serS/dmsA	seryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand): ref base G (144/140); new base A (7/11); total (151/151)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-04

Only 100 of 304 total aligned reads displayed.

GTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGAGCAAGCTTGAGTGAGCCATTATGAAAACGAAAATCCCTGATGCGGTATTGGCTGCTGAGGTGAGTCGCCGTGGTTTGGTAAAAAC  >  NC_000913/940488‑941023

        tctTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGtct                                                                                                                                                                                                                                                               >  2:639119/1‑275 (MQ=255)

                   cAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATa                                                                                                                                                                                                                                                     <  2:171643/274‑1 (MQ=255)

                                   aaGGCACGTGGTATGCAGGAAAGTTGCGGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                                                                                                                                                    <  2:523476/274‑1 (MQ=255)

                                                                cagcaAGTCGGACAAGAAAACCCGTCTGGGTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTg                                                                                                                                                                                                      <  1:111490/276‑1 (MQ=255)

                                                                                                      tGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTATCCTTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCCATTTAACATCTACTTATCTTTTTGATTTTCAAattattgta                                                                                                                                                                 >  4:437046/1‑272 (MQ=255)

                                                                                                            gTTCTGGGTTGGGGGTTGGTCGTTCGCTGGTTGCAGTAATGGAAAACTATCGGCATGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTTCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGtata                                                                                                                                                           <  4:579425/275‑1 (MQ=255)

                                                                                                              tctggtctggCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                                                                                                                                                         <  2:437050/275‑1 (MQ=255)

                                                                                                                       gCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑ttttt                                                                                                                                                                                                                                                                     >  2:774045/1‑158 (MQ=255)

                                                                                                                       gCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATg                                                                                                                                                               >  2:579429/1‑260 (MQ=255)

                                                                                                                           ttGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGgcgc                                                                                                                                                                                                                                                                           >  6:1019498/1‑149 (MQ=255)

                                                                                                                                cGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGC‑ttttt                                                                                                                                                                                                                                                                     >  7:38411/1‑149 (MQ=255)

                                                                                                                                 gTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgcgc                                                                                                                                                                                                                                                                           <  8:477810/143‑1 (MQ=255)

                                                                                                                                 gTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgcgc                                                                                                                                                                                                                                                                           >  7:477810/1‑143 (MQ=255)

                                                                                                                                 gTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGC‑ttttta                                                                                                                                                                                                                                                                    <  6:995615/149‑2 (MQ=255)

                                                                                                                                         gTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTa                                                                                                                                                                                                                                       >  4:191598/1‑170 (MQ=255)

                                                                                                                                            gCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCttt                                                                                                                                                                                                                                                         <  7:1479623/149‑1 (MQ=255)

                                                                                                                                            gCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCttt                                                                                                                                                                                                                                                         <  8:1409878/149‑1 (MQ=255)

                                                                                                                                              aGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGa                                                                                                                                                                                                                                                       <  8:1060469/149‑1 (MQ=255)

                                                                                                                                                   tGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTaca                                                                                                                      >  2:349196/1‑273 (MQ=255)

                                                                                                                                                    ggAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGa                                                                                                                                                                                                                                                 >  6:1422501/1‑149 (MQ=255)

                                                                                                                                                    ggAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGa                                                                                                                                                                                                                                                 >  5:427502/1‑149 (MQ=255)

                                                                                                                                                     gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                                                                >  6:904846/1‑149 (MQ=255)

                                                                                                                                                     gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                                                                >  6:424619/1‑149 (MQ=255)

                                                                                                                                                        aaCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACaa                                                                                                                                                                                                                                             >  5:1654951/1‑149 (MQ=255)

                                                                                                                                                        aaCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCGGGCGC‑TTTTTTTTTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATACCGAACGGCATTAGAAATAATTacacctgttt                                                                                                               >  4:418420/1‑269 (MQ=255)

                                                                                                                                                           tATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGC‑tttt                                                                                                                                                                                                                                                                      >  7:1438012/1‑121 (MQ=255)

                                                                                                                                                            aTCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGa                                                                                                                                <  2:418424/254‑1 (MQ=255)

                                                                                                                                                                gcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTAct                                                                                                                                                                                                                                     >  7:30968/1‑149 (MQ=255)

                                                                                                                                                                 cagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCCAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                                                                                                                                                    <  7:1223820/149‑1 (MQ=255)

                                                                                                                                                                  agGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcc                                                                                                                                                                                                                                   >  6:1442696/1‑149 (MQ=255)

                                                                                                                                                                   ggCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                                                                                                                                                                                                  <  6:171967/149‑1 (MQ=255)

                                                                                                                                                                    gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAGCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                                                                                                                                                    <  1:51258/146‑1 (MQ=255)

                                                                                                                                                                    gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGCCTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCTTTTGATANCGAACAATNATTActc                                                                                                                                                                                                                                    >  3:51258/1‑146 (MQ=255)

                                                                                                                                                                    gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATAcc                                                                                                                                                                                                                                                   >  4:702510/1‑131 (MQ=255)

                                                                                                                                                                     cTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                                                                                                                                                                                                  <  3:252216/147‑1 (MQ=255)

                                                                                                                                                                     cTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                                                                                                                                                                                                  >  1:252216/1‑147 (MQ=255)

                                                                                                                                                                       gATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGaat                                                                                                <  2:340065/275‑1 (MQ=255)

                                                                                                                                                                            tCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTacac                                                                                                                                                                                                                         <  6:1354070/149‑1 (MQ=255)

                                                                                                                                                                                 aagaaacagcAGAAGTACTGCTTTCGTATATGAACGGACTGGAATATATTTGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACCATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGa                                                                                       <  4:241050/265‑1 (MQ=255)

                                                                                                                                                                                  tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGtata                                                                                                                                                           <  2:21091/205‑1 (MQ=255)

                                                                                                                                                                                    aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCttt                                                                                                                                                                                                                                                         >  8:631858/1‑109 (MQ=255)

                                                                                                                                                                                    aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTCGATACCGAACAATAATTACTCCTCACTTACACGTAAtact                                                                                                                                                                                                                 >  6:723611/1‑149 (MQ=255)

                                                                                                                                                                                      gTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGta                                                                                                                                                             >  4:505588/1‑199 (MQ=255)

                                                                                                                                                                                            gAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                                                                                                                                                                         <  7:186127/149‑1 (MQ=255)

                                                                                                                                                                                              aGTTCTGCGTTCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGt                                                                                                                                                                                                       <  5:261516/149‑1 (MQ=255)

                                                                                                                                                                                                   tGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaaa                                                                                                                                                                                                  >  7:710297/1‑149 (MQ=255)

                                                                                                                                                                                                   tGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAattatt                                                                                                                                                                    >  4:613633/1‑179 (MQ=255)

                                                                                                                                                                                                                gAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATctct                                                                                                                                                                                     >  7:1316445/1‑149 (MQ=255)

                                                                                                                                                                                                                  aCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt                                                                                                                                                                                   <  8:1316445/149‑1 (MQ=255)

                                                                                                                                                                                                                  aCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGAGCAAGCTTGAGTGAGCCATTATGAAAACGAAAAt                                                     >  2:241922/1‑275 (MQ=255)

                                                                                                                                                                                                                    ggACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGa                                                                                                                                                                                 >  6:292718/1‑149 (MQ=255)

                                                                                                                                                                                                                      aCTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAAt                                                                                                                                                                                                 <  1:215994/131‑1 (MQ=255)

                                                                                                                                                                                                                        tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAtttt                                                                                                                                                                             <  8:306578/149‑1 (MQ=255)

                                                                                                                                                                                                                        tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAtttt                                                                                                                                                                             >  8:751559/1‑149 (MQ=255)

                                                                                                                                                                                                                         ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGt                                                                                                                                                                                                                                                                 <  5:126823/64‑1 (MQ=255)

                                                                                                                                                                                                                          gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                                                                                                                                                    <  4:344495/92‑1 (MQ=255)

                                                                                                                                                                                                                          gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCa                                                                                                                                                                           >  6:1515899/1‑149 (MQ=255)

                                                                                                                                                                                                                          gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTGCGAGTGAAAATCTACCTATCTCTTTGATTTTCa                                                                                                                                                                           >  8:1483802/1‑149 (MQ=255)

                                                                                                                                                                                                                          gAATATATTGCCTAATATCGAATTTTTGTGAGTCTAAAAAGCGCCTTCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCa                                                                                                                                                                           >  6:1124235/1‑149 (MQ=255)

                                                                                                                                                                                                                           aaTATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaa                                                                                                                                                                          >  5:780786/1‑149 (MQ=255)

                                                                                                                                                                                                                            atatatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaaa                                                                                                                                                                         <  5:1614099/149‑1 (MQ=255)

                                                                                                                                                                                                                            atatatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaaa                                                                                                                                                                         <  5:1359803/149‑1 (MQ=255)

                                                                                                                                                                                                                                atTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAattat                                                                                                                                                                     <  8:1525991/149‑1 (MQ=255)

                                                                                                                                                                                                                                  tGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                                                                                                                                                                   <  5:292718/149‑1 (MQ=255)

                                                                                                                                                                                                                                  tGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGAAACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTGGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                                                                                                                                                                   <  5:498545/149‑1 (MQ=255)

                                                                                                                                                                                                                                   ggCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt                                                                                                                                                                                   <  6:643/132‑1 (MQ=255)

                                                                                                                                                                                                                                    gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGa                                                                                                                                >  3:144521/1‑182 (MQ=255)

                                                                                                                                                                                                                                      tAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCtt                                                                                                                                                                                                                                                          >  6:1647588/1‑58 (MQ=255)

                                                                                                                                                                                                                                      tAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                                                                                                                                                                         <  6:1173827/107‑1 (MQ=255)

                                                                                                                                                                                                                                      tAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                                                                                                                                                                         >  5:1173827/1‑107 (MQ=255)

                                                                                                                                                                                                                                         tACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAat                                                                                                                                                                        >  2:72435/1‑137 (MQ=255)

                                                                                                                                                                                                                                         tACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGtat                                                                                                                                                            <  8:349834/149‑1 (MQ=255)

                                                                                                                                                                                                                                            ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGTCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTTAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                                                                                                                                                         >  5:924289/1‑149 (MQ=255)

                                                                                                                                                                                                                                            ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                                                                <  5:845668/62‑1 (MQ=255)

                                                                                                                                                                                                                                            ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                                                                >  6:845668/1‑62 (MQ=255)

                                                                                                                                                                                                                                            ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                                                                                                                                             <  6:1375166/113‑1 (MQ=255)

                                                                                                                                                                                                                                            ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaaa                                                                                                                                                                         >  5:607643/1‑133 (MQ=255)

                                                                                                                                                                                                                                            ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                                                                                                                                                         >  5:1056549/1‑149 (MQ=255)

                                                                                                                                                                                                                                            ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                                                                                                                                                         >  8:393867/1‑149 (MQ=255)

                                                                                                                                                                                                                                            ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                                                                                                                                                         >  6:234830/1‑149 (MQ=255)

                                                                                                                                                                                                                                             cAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTc                                                                                                                                                                            >  5:1528185/1‑130 (MQ=255)

                                                                                                                                                                                                                                                  tttCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCta                                                                                                                                                   <  5:522422/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   ttCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtat                                                                                                                                                  <  5:1442696/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   ttCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtat                                                                                                                                                  >  7:265326/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      tGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtatata                                                                                                                                               >  5:1413733/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      tGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCCATGCATACAAGCCtatgat                                                                                                                                               <  8:1468035/149‑4 (MQ=255)

                                                                                                                                                                                                                                                       gAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAg                                                                                                                                              <  8:1254260/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       gAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTCGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAg                                                                                                                                              <  5:723611/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        aaTCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                                                                                                                                                         <  4:758628/137‑1 (MQ=255)

                                                                                                                                                                                                                                                          tCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGa                                                                                                                                           >  8:936671/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          tCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGa                                                                                                                                           >  8:655801/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            tAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                                                                                                                                             >  6:1609294/1‑97 (MQ=255)

                                                                                                                                                                                                                                                            tAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAc                                                                                                                                         >  7:944808/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            tAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCCTTGATTTTCAAATTATTCGATGTATACAAGcc                                                                                                                                                     >  6:521830/1‑137 (MQ=255)

                                                                                                                                                                                                                                                             aaaaaGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACt                                                                                                                                        >  7:63317/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             aaaaaGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACt                                                                                                                                        >  7:1123181/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               aaaGCGCCCGCAGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGc                                                                                                                                      >  8:1027286/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                aaGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCt                                                                                                                                     >  8:193190/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                aaGCGCCCGCAGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCt                                                                                                                                     >  6:406050/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 aGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCta                                                                                                                                    >  5:371115/1‑149 (MQ=255)

GTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGAGCAAGCTTGAGTGAGCCATTATGAAAACGAAAATCCCTGATGCGGTATTGGCTGCTGAGGTGAGTCGCCGTGGTTTGGTAAAAAC  >  NC_000913/940488‑941023

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