BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,918,020	A→C	5.6%	intergenic (+207/+25)	barA → / ← gudD	hybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,918,020	0	A	C	5.6%	784.2 / 27.6	224	intergenic (+207/+25)	barA/gudD	hybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand): ref base A (103/107); new base C (5/8); total (108/116)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-02

Only 100 of 225 total aligned reads displayed.

TGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAA  >  NC_000913/2917752‑2918276

                                                             aTGGGTTTCTTTTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGGCAAAGTGCCTGTGGTTATGCCGGATGCGGCGTTAACGCCTTAGCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCGGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcc                                                                                                                                                                                                <  4:118556/274‑1 (MQ=255)

                                                                           gggTAGAATTCGGCTTTTTATGAGCCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATTCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTg                                                                                                                                                                                 <  4:570207/273‑1 (MQ=255)

                                                                                                                 ggTGACAAAGTGCC‑‑TGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                         >  2:326468/1‑275 (MQ=255)

                                                                                                                          gTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATa                                                                                                                                                                                                                                                                <  5:750576/149‑1 (MQ=255)

                                                                                                                            gCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATaga                                                                                                                                                                                                                                                              >  7:36812/1‑149 (MQ=255)

                                                                                                                            gCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATaga                                                                                                                                                                                                                                                              <  7:449745/149‑1 (MQ=255)

                                                                                                                               tgtgTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGcc                                                                                                                                                                                                                                                           >  5:722742/1‑149 (MQ=255)

                                                                                                                                gtgtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCg                                                                                                                                                                                                                                                          >  7:1050490/1‑149 (MQ=255)

                                                                                                                                 tgtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCgg                                                                                                                                                                                                                                                         >  6:5495/1‑149 (MQ=255)

                                                                                                                                   tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGGTGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtt                                                                                                                                                                                                                                                       >  8:331108/1‑149 (MQ=38)

                                                                                                                                   tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtt                                                                                                                                                                                                                                                       <  8:1066288/149‑1 (MQ=255)

                                                                                                                                   tttATGCCGGAGGCGGCGTTAGCGCCTTATCCAGCCTACAAAACCAATTAGATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtt                                                                                                                                                                                                                                                       >  8:926038/1‑149 (MQ=255)

                                                                                                                                      aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttga                                                                                                                                                                                                                                                    <  5:70181/149‑3 (MQ=39)

                                                                                                                                      aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttg                                                                                                                                                                                                                                                    <  8:439363/149‑1 (MQ=255)

                                                                                                                                          cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggttt                                                                                                                                                                                                                                                >  7:194663/1‑149 (MQ=255)

                                                                                                                                             aTGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt                                                                                                                                                                                                                                             <  7:860452/149‑1 (MQ=255)

                                                                                                                                              tGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt                                                                                                                                                                                                                                             >  8:848781/1‑148 (MQ=255)

                                                                                                                                              tGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt                                                                                                                                                                                                                                             <  7:848781/148‑1 (MQ=255)

                                                                                                                                                 ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCt                                                                                                                                                                                                                                         >  6:526131/1‑149 (MQ=255)

                                                                                                                                                  gcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggttt                                                                                                                                                                                                                                                >  4:716789/1‑141 (MQ=255)

                                                                                                                                                  gcgTTAAAGCCTTATCCGGCCTACAAAACCAATTAAATTAAATGAATTAAAAAATGATTTAGACCGAGTTCATAGAACAACAGGCAATTTTGGGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGTGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACGGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGTTACAGCTCATGGGCTTTCATCACTTGAt                                                                                                          <  4:447932/275‑1 (MQ=255)

                                                                                                                                                    gTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTaa                                                                                                                                                                                                                                      <  5:955943/149‑1 (MQ=255)

                                                                                                                                                    gTTAACGCCCTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGCTAGAGCCGGTTTTGGTTTTCTGTCTTat                                                                                                                                                                                                                                      >  8:694439/1‑148 (MQ=255)

                                                                                                                                                      tAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg                                                                                                                                                                                                                                            <  2:716789/141‑1 (MQ=255)

                                                                                                                                                         cGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAc                                                                                                                                                                                                                                 <  5:584246/149‑1 (MQ=255)

                                                                                                                                                         cGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAc                                                                                                                                                                                                                                 <  6:6235/149‑1 (MQ=255)

                                                                                                                                                          gCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAcc                                                                                                                                                                                                                                >  5:535253/1‑149 (MQ=255)

                                                                                                                                                            cTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtttt                                                                                                                                                                                                                                                     <  8:288483/126‑1 (MQ=255)

                                                                                                                                                            cTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtttt                                                                                                                                                                                                                                                     >  7:288483/1‑126 (MQ=255)

                                                                                                                                                             ttATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCtt                                                                                                                                                                                                                                        >  1:45266/1‑138 (MQ=255)

                                                                                                                                                                   ggCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCAAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                                                                                                                                        >  8:1055521/1‑148 (MQ=255)

                                                                                                                                                                      cTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgcgc                                                                                                                                                                                                                    <  5:782763/149‑1 (MQ=255)

                                                                                                                                                                      cTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgcgc                                                                                                                                                                                                                    <  6:200239/149‑1 (MQ=255)

                                                                                                                                                                        aCAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCtt                                                                                                                                                                                                                  >  7:136246/1‑149 (MQ=255)

                                                                                                                                                                          aaaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGt                                                                                                                                                                                                                >  8:793497/1‑149 (MQ=255)

                                                                                                                                                                            aaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                                                                                                                                              <  7:645320/149‑1 (MQ=255)

                                                                                                                                                                             aCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAt                                                                                                                                                                                                             <  6:667086/149‑1 (MQ=255)

                                                                                                                                                                              ccAATTAAATTCAATGAATTAAAAAATGATGTAGTCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATc                                                                                                                                                                                                            <  8:235359/149‑1 (MQ=255)

                                                                                                                                                                              ccAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATc                                                                                                                                                                                                            <  6:736300/149‑1 (MQ=255)

                                                                                                                                                                                aaTTAAATTCAATTAATTAAAAAAGGATGTAGACCTGACACATAGCACAACAGGCTTTTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGTTAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGtt                                                                            <  4:326468/275‑1 (MQ=255)

                                                                                                                                                                                 aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGaa                                                                                                                                                                                                         <  7:1062485/149‑1 (MQ=255)

                                                                                                                                                                                 aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                                                                                                                                  >  2:446508/1‑204 (MQ=255)

                                                                                                                                                                                       ttCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTcca                                                                                                                                                                                                   >  7:143475/1‑149 (MQ=255)

                                                                                                                                                                                       ttCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTcca                                                                                                                                                                                                   >  6:1048171/1‑149 (MQ=255)

                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGt                                                                                                                                                                                                                                           >  4:197570/1‑106 (MQ=255)

                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgcgc                                                                                                                                                                                                                    <  6:635113/129‑1 (MQ=255)

                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAt                                                                                                                                                                                                             >  5:109203/1‑136 (MQ=255)

                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        <  8:454906/141‑1 (MQ=255)

                                                                                                                                                                                           gtttatcagagactgattTAGACCTGTNANATATCACAACAGGCTTTTTTGCGTTTCGCCTCAGTCTCAAACCGGNCCCAGATAGCGCCTGTTTTGGTTTTCTGTCGTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCCNCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCAc                                                                  <  3:60845/261‑1 (MQ=255)

                                                                                                                                                                                           aTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcca                                                                                                                                                                                               >  5:507735/1‑149 (MQ=255)

                                                                                                                                                                                           aTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcca                                                                                                                                                                                               >  5:297399/1‑149 (MQ=255)

                                                                                                                                                                                              aaTTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTa                                                                                                                                                                                                                                       <  2:197570/106‑1 (MQ=255)

                                                                                                                                                                                               aTTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGaa                                                                                                                                                                                           <  5:905248/149‑1 (MQ=255)

                                                                                                                                                                                                ttAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCAGCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAAt                                                                                                                                                                                          >  6:946112/1‑149 (MQ=255)

                                                                                                                                                                                                ttAAAAAATGATGTAGACCCGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAAt                                                                                                                                                                                          <  8:400650/149‑1 (MQ=255)

                                                                                                                                                                                                 tAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATc                                                                                                                                                                                         <  8:655802/149‑1 (MQ=255)

                                                                                                                                                                                                 tAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATc                                                                                                                                                                                         <  6:91091/149‑1 (MQ=255)

                                                                                                                                                                                                 tAAAAAACGATGTAAACCTGATACATAGCACAACAGGCTATTTTGCGTTTATCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTAAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGaaac                                                                                                                                                                                         <  6:574793/149‑3 (MQ=255)

                                                                                                                                                                                                               gACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTcc                                                                                                                                                                           >  7:471408/1‑149 (MQ=255)

                                                                                                                                                                                                                  cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcat                                                                                                                                                                        >  6:702283/1‑149 (MQ=255)

                                                                                                                                                                                                                  cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcat                                                                                                                                                                        >  5:552384/1‑149 (MQ=255)

                                                                                                                                                                                                                   tGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatc                                                                                                                                                                       <  6:409561/149‑1 (MQ=255)

                                                                                                                                                                                                                    gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCCCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGCTCCATATCGCTTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCaa                                        >  2:307218/1‑275 (MQ=255)

                                                                                                                                                                                                                     aTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTcc                                                                                                                                                                                                    >  4:687255/1‑118 (MQ=255)

                                                                                                                                                                                                                     aTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgc                                                                                                                                                                     <  6:818036/149‑1 (MQ=255)

                                                                                                                                                                                                                     aTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCCCTTGATc                                                                                                         >  1:192009/1‑209 (MQ=255)

                                                                                                                                                                                                                       aCATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGc                                                                                                                                                                                                                            <  5:586070/92‑1 (MQ=255)

                                                                                                                                                                                                                        cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgcatc                                                                                                                                                                  <  5:357544/149‑1 (MQ=255)

                                                                                                                                                                                                                        cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgcatc                                                                                                                                                                  <  5:838945/149‑1 (MQ=255)

                                                                                                                                                                                                                         aTAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcc                                                                                                                                                                                                <  2:687255/118‑1 (MQ=255)

                                                                                                                                                                                                                          tAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGt                                                                                                                                                                <  5:418805/149‑1 (MQ=255)

                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAcgcg                                                                                                                                                          >  7:408000/1‑149 (MQ=255)

                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAcgcg                                                                                                                                                          >  5:1016099/1‑149 (MQ=255)

                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                                                                                                                                        >  1:152238/1‑167 (MQ=255)

                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGGGTTTCt                                                                                                                                        >  1:76027/1‑167 (MQ=255)

                                                                                                                                                                                                                                 acaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCt                                                                       <  2:100627/231‑1 (MQ=255)

                                                                                                                                                                                                                                  caGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTaa                                                                                                                                                                                                                                      >  1:108848/1‑71 (MQ=255)

                                                                                                                                                                                                                                  caGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTcc                                                                                                                                                                                                    >  4:23075/1‑105 (MQ=255)

                                                                                                                                                                                                                                  caGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTgg                            >  2:624616/1‑273 (MQ=255)

                                                                                                                                                                                                                                     gCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCa                                                                                                                                                     >  5:251815/1‑149 (MQ=255)

                                                                                                                                                                                                                                      cTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACg                                                                                                                                                                                                                                    <  3:41914/69‑1 (MQ=255)

                                                                                                                                                                                                                                      cTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                                                                                                                                  <  3:108848/71‑1 (MQ=255)

                                                                                                                                                                                                                                        aTTTTCCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGGGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGgcgc                    <  2:530188/275‑1 (MQ=255)

                                                                                                                                                                                                                                         ttttGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGcc                                                                                                                                                 >  6:809436/1‑149 (MQ=255)

                                                                                                                                                                                                                                               gTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGtt                                                                                                                                           >  6:159648/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 ttAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                         <  5:406888/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  tAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                                                                                                                                        <  6:974156/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   aGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTg                                                                                                                                       >  7:316227/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    gCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGa                                                                                                                                      <  8:178773/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       tCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATAc                                                                                                                                   <  5:942736/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             tcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTccc                                                                                                                                                                          >  6:776899/1‑104 (MQ=255)

                                                                                                                                                                                                                                                             tcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCa                                                                                                                             >  7:769630/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             tcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCa                                                                                                                             >  8:148108/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             caaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCa                                                                                                                             >  8:729819/3‑149 (MQ=255)

                                                                                                                                                                                                                                                                aaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATgg                                                                                                                          >  8:881149/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtc                                                                                                                    <  5:677546/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtc                                                                                                                    >  6:382354/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       tCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca                                                                                                                   >  8:619481/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        ccAGATGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtcat                                                                                                                  <  6:473763/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                         cAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca                                                                                                                   <  2:206460/147‑1 (MQ=255)

TGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAA  >  NC_000913/2917752‑2918276

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