BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,517,484	C→G	9.0%	I32M (ATC→ATG)	ydcX →	DUF2566 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,517,484	0	C	G	9.0%	338.9 / 18.5	99	I32M (ATC→ATG)	ydcX	DUF2566 family protein
Reads supporting (aligned to +/- strand): ref base C (47/43); new base G (5/4); total (52/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 102 total aligned reads displayed.

AAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGT  >  NC_000913/1517336‑1517632

aaaaaaTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATc                                                                                                                                                      >  4:472405/1‑149 (MQ=255)

 aaaaaTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATcc                                                                                                                                                     <  3:389626/149‑1 (MQ=255)

        aCGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCtt                                                                                                                                              <  4:196311/149‑1 (MQ=255)

          gCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTaa                                                                                                                                            >  4:185886/1‑149 (MQ=255)

           cagatTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTCTGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatc                                                                                                                                           <  4:127368/146‑3 (MQ=255)

           cATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatc                                                                                                                                           <  3:352854/149‑3 (MQ=255)

            atatTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgc                                                                                                                                          <  1:682393/149‑1 (MQ=255)

                        gTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggtg                                                                                                                              >  1:94888/1‑149 (MQ=255)

                        gTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggtg                                                                                                                              <  1:197100/149‑1 (MQ=255)

                        gTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggtg                                                                                                                              <  1:745937/149‑1 (MQ=255)

                         ttttATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggtgg                                                                                                                             >  2:777827/1‑149 (MQ=255)

                          tttATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTggg                                                                                                                            <  1:721387/149‑1 (MQ=255)

                          tttATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTggg                                                                                                                            <  3:183168/149‑1 (MQ=255)

                           ttATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTgggg                                                                                                                           <  2:443716/149‑1 (MQ=255)

                             aTTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGc                                                                                                                         >  3:693740/1‑149 (MQ=255)

                              ttCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatc                                                                                                                                           >  4:681596/1‑128 (MQ=255)

                              ttCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatc                                                                                                                                           <  3:681596/130‑3 (MQ=255)

                                 aTCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACa                                                                                                                     >  3:818795/1‑149 (MQ=255)

                                   ctggctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATg                                                                                                                   >  3:609010/1‑149 (MQ=255)

                                     ggctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGc                                                                                                                 >  3:285266/1‑149 (MQ=255)

                                        tggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGa                                                                                                              >  3:767639/1‑149 (MQ=255)

                                        tggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGa                                                                                                              <  3:40383/149‑1 (MQ=255)

                                         ggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGAt                                                                                                             <  2:163199/149‑1 (MQ=255)

                                         ggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGAt                                                                                                             <  4:621401/149‑1 (MQ=255)

                                             gCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGt                                                                                                         <  3:157702/149‑1 (MQ=255)

                                             gCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGt                                                                                                         >  2:884622/1‑149 (MQ=255)

                                               ccGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGttt                                                                                                       >  3:804836/1‑149 (MQ=255)

                                                 gTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTcc                                                                                                     <  3:844363/149‑1 (MQ=255)

                                                 gTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTcc                                                                                                     <  2:673679/149‑1 (MQ=255)

                                                     aTGTCTCTCTTTCCACACAGTCTTGTTTTTTATGCGGTTACGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGAGAAGAAACGCCTCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATAAGTTTCCCGGt                                                                                                 <  2:221937/149‑1 (MQ=255)

                                                      tGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTg                                                                                                <  1:884622/149‑1 (MQ=255)

                                                      tGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTg                                                                                                <  1:174970/149‑1 (MQ=255)

                                                          tgtctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCCCCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGgcgc                                                                                            <  4:758311/147‑1 (MQ=255)

                                                          tctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTgtgg                                                                                                                           >  4:776200/1‑118 (MQ=255)

                                                          tctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTgggg                                                                                                                           <  3:776200/118‑1 (MQ=255)

                                                            tctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtg                                                                                          >  1:884878/1‑149 (MQ=255)

                                                            tctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtg                                                                                          >  1:106433/1‑149 (MQ=255)

                                                              tATCAACACATGCTTGTTTTTTATGCGGTTATGTCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtgtt                                                                                        <  4:584572/149‑1 (MQ=255)

                                                               aTCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtgttg                                                                                       >  4:314078/1‑149 (MQ=255)

                                                                tCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAGGATAAGAAACGCATCGGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtgttgt                                                                                      >  4:700748/1‑149 (MQ=255)

                                                                  aaCACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGttt                                                                                    >  2:100336/1‑149 (MQ=255)

                                                                   acacATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGtttt                                                                                   <  3:713816/149‑1 (MQ=255)

                                                                     acaTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCa                                                                                 >  2:12529/1‑149 (MQ=255)

                                                                     acaTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCa                                                                                 >  4:796157/1‑149 (MQ=255)

                                                                      caTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCAc                                                                                >  1:159331/1‑149 (MQ=255)

                                                                      caTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCAc                                                                                >  2:749482/1‑149 (MQ=255)

                                                                           ttGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGttt                                                                           >  4:616779/1‑149 (MQ=255)

                                                                              ttttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaa                                                                        <  3:96378/149‑1 (MQ=255)

                                                                               tttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                                                                       <  1:434844/149‑1 (MQ=255)

                                                                               tttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                                                                       <  2:522889/149‑1 (MQ=255)

                                                                               tttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                                                                       <  4:285266/149‑1 (MQ=255)

                                                                               tttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                                                                       <  3:339615/149‑1 (MQ=255)

                                                                               tttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                                                                       <  3:554673/149‑1 (MQ=255)

                                                                                ttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAAt                                                                      <  1:777827/149‑1 (MQ=255)

                                                                                    aTGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTCCACTGTTTTAAATCTTa                                                                  >  2:750910/1‑149 (MQ=255)

                                                                                     tGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTAc                                                                 <  3:616779/149‑1 (MQ=255)

                                                                                          ttATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTCCACTGTTTTAAATCTTACTCACt                                                            <  1:750910/149‑1 (MQ=255)

                                                                                              ggCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCt                                                        >  4:735902/1‑149 (MQ=255)

                                                                                              ggCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCt                                                        >  4:551938/1‑149 (MQ=255)

                                                                                              ggCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCt                                                        <  4:417431/149‑1 (MQ=255)

                                                                                              ggCAGCAACCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCt                                                        >  3:814365/1‑149 (MQ=255)

                                                                                               gcagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCtt                                                       >  2:239979/1‑149 (MQ=255)

                                                                                                 agcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGa                                                     >  1:808308/1‑149 (MQ=255)

                                                                                                  gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGAc                                                    <  3:472405/149‑1 (MQ=255)

                                                                                                  gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGAc                                                    <  4:609010/149‑1 (MQ=255)

                                                                                                   caATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACg                                                   >  2:375345/1‑149 (MQ=255)

                                                                                                         cATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGcacca                                             <  2:884878/149‑1 (MQ=255)

                                                                                                           tttCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACg                                           >  3:223409/1‑149 (MQ=255)

                                                                                                            ttCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGc                                          >  1:467920/1‑149 (MQ=255)

                                                                                                            ttCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGc                                          <  1:848596/149‑1 (MQ=255)

                                                                                                                  tCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtgg                                                                                                                                <  1:230731/57‑1 (MQ=255)

                                                                                                                  tCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtgg                                                                                                                                >  2:230731/1‑57 (MQ=255)

                                                                                                                    acttGGTTTCTCTCTCACGATAAGAAACGCATCCGTGTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCTCTGGCTTGTCGGCACCACGCCGACTtat                                  >  2:31363/4‑149 (MQ=255)

                                                                                                                    aCCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTtat                                  >  4:863358/1‑149 (MQ=255)

                                                                                                                    aCCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTtat                                  <  1:375345/149‑1 (MQ=255)

                                                                                                                     ccTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTtata                                 >  3:189243/1‑149 (MQ=255)

                                                                                                                      cTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTatat                                >  1:811114/1‑149 (MQ=255)

                                                                                                                       tGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATc                               >  1:181619/1‑149 (MQ=255)

                                                                                                                              tttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGa                        >  4:79747/1‑149 (MQ=255)

                                                                                                                               ttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGAt                       <  2:885631/149‑1 (MQ=255)

                                                                                                                                  tcACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGaa                    >  2:340720/1‑149 (MQ=255)

                                                                                                                                   cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAAt                   <  4:818147/149‑1 (MQ=255)

                                                                                                                                   cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAAt                   <  1:239979/149‑1 (MQ=255)

                                                                                                                                    aCGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATc                  >  1:556766/1‑149 (MQ=255)

                                                                                                                                    aCGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATc                  >  2:736222/1‑149 (MQ=255)

                                                                                                                                      gATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAAt                   <  4:267982/146‑1 (MQ=255)

                                                                                                                                      gATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAAt                   >  3:267982/1‑146 (MQ=255)

                                                                                                                                      gATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGa                >  1:177944/1‑149 (MQ=255)

                                                                                                                                       ttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTa                                                                  >  4:485154/2‑98 (MQ=255)

                                                                                                                                       ttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTa                                                                  <  3:485154/97‑1 (MQ=255)

                                                                                                                                       ttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGAt               >  3:477275/2‑149 (MQ=255)

                                                                                                                                       ttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGAt               >  4:849902/2‑149 (MQ=255)

                                                                                                                                       ttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGAt               >  2:522862/2‑149 (MQ=255)

                                                                                                                                           gAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAAt           >  2:142146/1‑149 (MQ=255)

                                                                                                                                           gAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAAt           <  2:467920/149‑1 (MQ=255)

                                                                                                                                                gCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaa                                                                        <  1:220772/83‑1 (MQ=255)

                                                                                                                                                gCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaa                                                                        >  2:220772/1‑83 (MQ=255)

                                                                                                                                                  aTCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACaa    >  2:261762/1‑149 (MQ=255)

                                                                                                                                                  aTCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACaa    <  1:855917/149‑1 (MQ=255)

                                                                                                                                                    ccGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGt  <  2:748421/149‑1 (MQ=255)

AAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGT  >  NC_000913/1517336‑1517632

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