BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,838,222	A→T	6.6%	intergenic (‑117/+32)	hydN ← / ← ascG	formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,838,222	0	A	T	6.6%	522.5 / 20.0	151	intergenic (‑117/+32)	hydN/ascG	formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor
Reads supporting (aligned to +/- strand): ref base A (66/74); new base T (5/5); total (72/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 155 total aligned reads displayed.

GAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGT  >  NC_000913/2837948‑2838488

gAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAAAACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGa                                                                                                                                                                                                                                                                            <  4:565162/275‑1 (MQ=255)

      ttcaacccTTTCAGGGCATGGGTCCTCGGGAAAAAATTTTCCTGGGTCAGCGATGCACCGTCCTGATTTTCCTTATGGGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTGATGCAGGAAAAACGATGTCATCGACACTAg                                                                                                                                                                                                                                                                      <  4:72464/269‑1 (MQ=255)

                       aTGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAACCGATGTCATCACCAATAGTGACGATGACATTtga                                                                                                                                                                                                                                                      >  4:48965/1‑274 (MQ=255)

                                                                            ccTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTcc                                                                                                                                                                                                >  2:628542/1‑275 (MQ=255)

                                                                                                                    aCCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtga                                                                                                                                                                                                                     >  2:357787/1‑214 (MQ=255)

                                                                                                                        aTACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTcg                                                                                                                                                                                                                 <  4:357787/214‑1 (MQ=255)

                                                                                                                                aCTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGaca                                                                                                                                                                                                                                                                          >  5:286057/1‑149 (MQ=255)

                                                                                                                                  tCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACt                                                                                                                                                                                                                                                                        >  6:850737/1‑149 (MQ=255)

                                                                                                                                  tCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACt                                                                                                                                                                                                                                                                        >  8:202651/1‑149 (MQ=255)

                                                                                                                                  tCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACt                                                                                                                                                                                                                                                                        >  5:904045/1‑149 (MQ=255)

                                                                                                                                    gcgTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAg                                                                                                                                                                                                                                                                      >  6:263814/1‑149 (MQ=255)

                                                                                                                                     cgTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGt                                                                                                                                                                                                                                                                     <  7:820066/149‑1 (MQ=255)

                                                                                                                                        cAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGAc                                                                                                                                                                                                                                                                  <  8:381408/149‑1 (MQ=255)

                                                                                                                                          gCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGa                                                                                                                                                                                                                                                                <  5:397861/149‑1 (MQ=255)

                                                                                                                                           cAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGAt                                                                                                                                                                                                                                                               <  8:419119/149‑1 (MQ=255)

                                                                                                                                           cAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGAt                                                                                                                                                                                                                                                               <  8:550721/149‑1 (MQ=255)

                                                                                                                                           cAATCATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTATTGCGCCGATACATTCACAAATCATGGCAGTTTTCAAATTACTGTTCTTTCAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACTAt                                                                                                                                                                                                                                                               <  8:672112/149‑1 (MQ=255)

                                                                                                                                            aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                                                                                                                                                                              <  5:983916/149‑1 (MQ=255)

                                                                                                                                            aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                                                                                                                                                                              >  7:391730/1‑149 (MQ=255)

                                                                                                                                              tgatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGAc                                                                                                                                                                                                                                                            <  7:754304/149‑1 (MQ=255)

                                                                                                                                               gatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTGATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACa                                                                                                                                                                                                                                                           <  7:670521/149‑1 (MQ=255)

                                                                                                                                               gatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACa                                                                                                                                                                                                                                                           <  7:864596/149‑1 (MQ=255)

                                                                                                                                                atgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAt                                                                                                                                                                                                                                                          <  6:71337/149‑1 (MQ=255)

                                                                                                                                                 tgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAtc                                                                                                                                                                                                                                                         <  7:884442/149‑2 (MQ=255)

                                                                                                                                                  gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgt                                                                                                                                                                                                                                                        >  8:174714/1‑149 (MQ=255)

                                                                                                                                                   aaaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTGTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                                                                                                                                       >  8:748703/1‑149 (MQ=255)

                                                                                                                                                    aaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                                                                                                                                                                                                      <  8:332095/149‑1 (MQ=255)

                                                                                                                                                       ggTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                                                                                                                                       <  5:638083/145‑1 (MQ=255)

                                                                                                                                                       ggTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                                                                                                                                       >  6:638083/1‑145 (MQ=255)

                                                                                                                                                       ggTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatga                                                                                                                                                                                                                                                   <  7:782912/149‑1 (MQ=255)

                                                                                                                                                       ggTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatga                                                                                                                                                                                                                                                   >  7:474610/1‑149 (MQ=255)

                                                                                                                                                            atcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATg                                                                                                                                                                                                                                              >  5:1019339/1‑149 (MQ=255)

                                                                                                                                                             tcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                                                                                                                                                                                                      >  6:795356/1‑140 (MQ=255)

                                                                                                                                                             tcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                                                                                                                                                                                                      <  5:795356/140‑1 (MQ=255)

                                                                                                                                                             tcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGt                                                                                                                                                                                                                                             >  5:705863/1‑149 (MQ=255)

                                                                                                                                                               aaaTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGttt                                                                                                                                                                                                                                           <  5:809872/149‑1 (MQ=255)

                                                                                                                                                                 aTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTAt                                                                                                                                                                                                                                         <  5:35732/149‑1 (MQ=255)

                                                                                                                                                                  ttCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATc                                                                                                                                                                                                                                        >  8:27028/1‑149 (MQ=255)

                                                                                                                                                                       cAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTTGTGACGATGACATGTGATGACAATGTTTATCGCGaa                                                                                                                                                                                                                                   <  8:1022702/149‑1 (MQ=255)

                                                                                                                                                                       cAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGaa                                                                                                                                                                                                                                   <  6:13454/149‑1 (MQ=255)

                                                                                                                                                                            aTGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTTTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGc                                                                                                                                                                                                                              <  6:526839/149‑1 (MQ=255)

                                                                                                                                                                              gACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCaa                                                                                                                                                                                                                            <  6:904045/149‑1 (MQ=255)

                                                                                                                                                                               aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAt                                                                                                                                                                                                                           <  7:950695/149‑1 (MQ=255)

                                                                                                                                                                                 aGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtga                                                                                                                                                                                                                         <  5:774962/149‑1 (MQ=255)

                                                                                                                                                                                 aGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtga                                                                                                                                                                                                                         <  5:850737/149‑1 (MQ=255)

                                                                                                                                                                                  gTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                                                                                                                                       <  5:944782/118‑1 (MQ=255)

                                                                                                                                                                                  gTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgag                                                                                                                                                                                                                        >  7:355737/1‑149 (MQ=255)

                                                                                                                                                                                   ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgat                                                                                                                                                                                                                                                     <  7:972706/119‑1 (MQ=255)

                                                                                                                                                                                   ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagt                                                                                                                                                                                                                       <  6:382984/149‑1 (MQ=255)

                                                                                                                                                                                     aaTGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtga                                                                                                                                                                                                                     <  5:36422/149‑1 (MQ=255)

                                                                                                                                                                                          gcCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCg                                                                                  <  2:8171/275‑1 (MQ=255)

                                                                                                                                                                                           ccGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAAGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCgg                                                                                 >  2:121750/1‑275 (MQ=255)

                                                                                                                                                                                            cGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatga                                                                                                                                                                                                                                                   <  7:1009000/112‑1 (MQ=255)

                                                                                                                                                                                             taTACATTCACAAATTATGCCAGTTTTGAATTTACTGTTATTTAAGGAAATTAATTTCTGGAATGGAGGAAAAACGATATAATCGACACTAGTGACGATGACATGTGATGAGAATGGTTATCGCGAAGGGGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGaa                                                                               <  4:628542/274‑1 (MQ=255)

                                                                                                                                                                                             cagaaatttaaaaaTCAGGCCATGTTTTAAGTTACGTTTATTTAATGAAATTAATTTTTGTAATGCAGGAAAAACGATGTGATCTACACTATGGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGGGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGaa                                                                               <  4:54117/264‑1 (MQ=255)

                                                                                                                                                                                                 cATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGAt                                                                                                                                                                                                         <  6:293678/149‑1 (MQ=255)

                                                                                                                                                                                                     cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGt                                                                                                                                                                                                     <  8:721087/149‑1 (MQ=255)

                                                                                                                                                                                                         aaTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTc                                                                                                                                                                                                 >  5:1026970/1‑149 (MQ=255)

                                                                                                                                                                                                         aaTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTc                                                                                                                                                                                                 <  5:891177/149‑1 (MQ=255)

                                                                                                                                                                                                          aTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTcc                                                                                                                                                                                                <  7:202651/149‑1 (MQ=255)

                                                                                                                                                                                                              tGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGgaga                                                                                        >  2:232199/1‑249 (MQ=255)

                                                                                                                                                                                                                 cAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                                                                                                                                                                                         >  8:209473/1‑149 (MQ=255)

                                                                                                                                                                                                                  aGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAgcgc                                                                                    <  4:232199/249‑1 (MQ=255)

                                                                                                                                                                                                                    tttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCa                                                                                                                                                                                                       >  5:177358/1‑132 (MQ=255)

                                                                                                                                                                                                                    tttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCa                                                                                                                                                                                                       <  6:177358/132‑1 (MQ=255)

                                                                                                                                                                                                                      tttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGaa                                                                                                                                                                                                                                   <  7:365641/102‑1 (MQ=255)

                                                                                                                                                                                                                      tttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGtttt                                                                                                                                                                                    <  5:567540/149‑1 (MQ=255)

                                                                                                                                                                                                                              aCTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGTGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGaa                                                >  2:220498/1‑273 (MQ=255)

                                                                                                                                                                                                                              aCTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTg                                                                                                                                                                            <  6:839143/149‑1 (MQ=255)

                                                                                                                                                                                                                              aCTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGGCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTCCGCTGGAAAGGGCCGGAACGGTGTAGGGGGAAAGAGCGATATCATCGaag                                               >  4:216414/1‑273 (MQ=255)

                                                                                                                                                                                                                               cTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTga                                                                                                                                                                           <  8:116996/149‑1 (MQ=255)

                                                                                                                                                                                                                                tGTTATTTAAGGAAATTAATTTCTCTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGttt                                                                                                                                                                                                                                           <  8:910568/84‑1 (MQ=255)

                                                                                                                                                                                                                                 gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTgaga                                                                                                                                                                         >  5:116833/1‑149 (MQ=255)

                                                                                                                                                                                                                                 gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCCGCGGTgaga                                                                                                                                                                         >  5:202918/1‑149 (MQ=255)

                                                                                                                                                                                                                                  ttATTTAAGGAAATGAATTTATGTAATGCAGGAAAAACGATGTCATCGACACTAGGGACGATGAGATGTGATGACAAAGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGaa                                                                               <  1:148298/238‑1 (MQ=255)

                                                                                                                                                                                                                                       tAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATccc                                                                                                                                                                   <  5:518032/149‑1 (MQ=255)

                                                                                                                                                                                                                                         aGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCg                                                                                                                                                                 >  8:692677/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 aattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAg                                                                                                                                                         <  8:381849/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       tGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATTACTTACACCTTCTCTGGCa                      >  2:355524/1‑274 (MQ=255)

                                                                                                                                                                                                                                                         tAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATaaaaa                                                                                                                                                 <  8:355737/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             gCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAg                                                                                                                                             >  6:376692/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             gCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAg                                                                                                                                             <  5:74711/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGt                                                                                                                                         >  6:757175/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGt                                                                                                                                         >  8:52836/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGt                                                                                                                                         >  8:844118/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    aaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCa                                                                                                                                      >  5:589254/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     tcacaTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaa                                                                                                                                     >  8:153621/5‑149 (MQ=255)

                                                                                                                                                                                                                                                                     tcacaTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaa                                                                                                                                     >  7:727783/5‑149 (MQ=255)

                                                                                                                                                                                                                                                                         aTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaataat                                                                                                                                 >  8:816363/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                          tGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGa                                                                                                                                                          <  5:522328/123‑1 (MQ=255)

                                                                                                                                                                                                                                                                          tGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGa                                                                                                                                                          >  6:522328/1‑123 (MQ=255)

                                                                                                                                                                                                                                                                            tCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTc                                                                                                                              <  6:421472/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            tCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTc                                                                                                                              <  7:1062416/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                              aTCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATc                                                                                                                                                                     <  8:1027212/108‑1 (MQ=255)

                                                                                                                                                                                                                                                                              aTCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCt                                                                                                                            >  5:409391/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                              aTCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGCGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCt                                                                                                                            >  8:270721/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                              aTCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCt                                                                                                                            >  5:597400/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                               tCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTAcgcc            >  4:46570/1‑260 (MQ=255)

                                                                                                                                                                                                                                                                                  acacTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAat                                                                                                                        <  6:949213/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                                  acacTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAat                                                                                                                        <  5:263923/149‑1 (MQ=255)

GAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGT  >  NC_000913/2837948‑2838488

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