BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,807,093	A→T	5.0%	intergenic (+19/‑39)	proW → / → proX	glycine betaine/proline ABC transporter permease/glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,807,093	0	A	T	5.0%	1182.7 / 35.5	332	intergenic (+19/‑39)	proW/proX	glycine betaine/proline ABC transporter permease/glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein
Reads supporting (aligned to +/- strand): ref base A (169/144); new base T (9/8); total (180/152)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

Only 100 of 334 total aligned reads displayed.

GGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTTCGCGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGG  >  NC_000913/2806831‑2807361

       tCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTTCGCGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCCCGCCATTATCCTCGATTGTCTGACGCACGCCGTTGGGCGCGACTCCCGCAGCCGCGGCAACCGTCGCCGGGACAACACTGGCCCATCTGGCCTGCCGACCCCCCCAGTCAATCAGTAAACTTGCCCNGTCCCNNGGGCGGCGGGCAATAACCCCCCGCCa                                                                                                                                                                                                                                                       >  3:23255/1‑279 (MQ=255)

                    tCGCCCTGATGATTGCCGGCGGCGGGTTGGGTCGGTTGGTCCGTCGCGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATaa                                                                                                                                                                                                                                          <  1:23255/279‑1 (MQ=255)

                        cTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTTCGCGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCTCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATaa                                                                                                                                                                                                                                          <  2:127353/275‑1 (MQ=255)

                                                            aCTTCGCGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGGACTTTTTGCGACAGCGtt                                                                                                                                                                                                        >  2:253534/1‑273 (MQ=255)

                                                                cgcgGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTCTCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCAc                                                                                                                                                                                                  >  2:113588/1‑275 (MQ=255)

                                                                                           tttGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTtgctg                                                                                                                                                                        <  2:410399/273‑1 (MQ=255)

                                                                                                              tCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGc                                                                                                                                                                                                                                                                           >  2:344841/1‑156 (MQ=255)

                                                                                                                        gATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCgg                                                                                                                                                                                                                                                                        >  8:1635330/1‑149 (MQ=255)

                                                                                                                        gATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCgg                                                                                                                                                                                                                                                                        >  5:280973/1‑149 (MQ=255)

                                                                                                                          tCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGc                                                                                                                                                                                                                                                                      >  8:557552/1‑149 (MQ=255)

                                                                                                                          tCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGc                                                                                                                                                                                                                                                                      <  6:815652/149‑1 (MQ=255)

                                                                                                                            cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaa                                                                                                                                                                                                                                                                    <  5:209390/149‑2 (MQ=255)

                                                                                                                              cGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagt                                                                                                                                                                                                                                                                  <  8:745109/149‑4 (MQ=255)

                                                                                                                                ccATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATc                                                                                                                                                                                                                                                                <  6:1222580/149‑1 (MQ=255)

                                                                                                                                 cATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCa                                                                                                                                                                                                                                                                     <  6:1075816/143‑1 (MQ=255)

                                                                                                                                 aattttcccgtgtCGTCTAGCGCAGTCCGTTGGGGGCGAGTCACGCAGTCGCGGCAACCGTGGCTGTTACACAACCGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCACGGTGACGCCGGGCATTATCACCCTGCCAAAAACAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTGTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGtt                                                                                                                                  <  4:164323/262‑1 (MQ=255)

                                                                                                                                     aTCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCt                                                                                                                                                                                                                                                           <  7:760387/149‑1 (MQ=255)

                                                                                                                                     aTCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCt                                                                                                                                                                                                                                                           >  5:817907/1‑149 (MQ=255)

                                                                                                                                         tCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCa                                                                                                                                                                                                                                                       >  6:849162/1‑149 (MQ=255)

                                                                                                                                           gATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAcc                                                                                                                                                                                                                                                             <  1:196949/141‑1 (MQ=255)

                                                                                                                                            aTCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaa                                                                                                                                                                                                                                                    >  6:1508150/1‑149 (MQ=255)

                                                                                                                                               gTCTGACGCAGGCCGTGGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaaaaa                                                                                                                                                                                                                                                 <  5:988842/149‑1 (MQ=255)

                                                                                                                                                tCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACttttt                                                                                                                                                                                                                   >  4:613182/1‑178 (MQ=255)

                                                                                                                                                tCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc                                                                                                                                                                                         >  4:381805/1‑204 (MQ=255)

                                                                                                                                                tCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGCCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTGTTCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAg                                                                                                                                                                                                                                                <  8:1037430/149‑1 (MQ=255)

                                                                                                                                                 cTGACGCGGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAgg                                                                                                                                                                                                                                               <  6:349332/149‑1 (MQ=255)

                                                                                                                                                      gCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATaa                                                                                                                                                                                                                                          <  6:1690528/149‑1 (MQ=255)

                                                                                                                                                        aGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACa                                                                                                                                                                                                                                        <  6:487447/149‑1 (MQ=255)

                                                                                                                                                          gCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAAt                                                                                                                                                                                                                                      <  5:517582/149‑1 (MQ=255)

                                                                                                                                                          gCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaa                                                                                                                                                                                                                                                                    <  6:1420509/119‑2 (MQ=255)

                                                                                                                                                          gCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaa                                                                                                                                                                                                                                                                    >  5:1420509/1‑118 (MQ=255)

                                                                                                                                                           ccGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAcc                                                                                                                                                                                                                                                             <  7:634311/125‑1 (MQ=255)

                                                                                                                                                           ccGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAcc                                                                                                                                                                                                                                                             >  8:634311/1‑125 (MQ=255)

                                                                                                                                                            cGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGc                                                                                                                                                                                                                                    <  7:689310/149‑1 (MQ=255)

                                                                                                                                                               tGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGa                                                                                                   <  2:577253/275‑1 (MQ=255)

                                                                                                                                                                 ggCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaaaaa                                                                                                                                                                                                                                                 <  6:675295/131‑1 (MQ=255)

                                                                                                                                                                 ggCGCGACTCACGCAGTCGCCGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaaaaa                                                                                                                                                                                                                                                 >  5:675295/1‑131 (MQ=255)

                                                                                                                                                                    gcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAAt                                                                                                                                                                                                                                            >  1:120712/1‑133 (MQ=255)

                                                                                                                                                                       aCTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTa                                                                                                                                                                                                                         <  8:157339/149‑1 (MQ=255)

                                                                                                                                                                           aCGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACttt                                                                                                                                                                                                                     >  8:316933/1‑149 (MQ=255)

                                                                                                                                                                           aCGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACttt                                                                                                                                                                                                                     >  5:839580/1‑149 (MQ=255)

                                                                                                                                                                             gCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACttttt                                                                                                                                                                                                                   >  8:383641/1‑149 (MQ=255)

                                                                                                                                                                                    cgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGaa                                                                                                                                                                                                                                             >  6:753069/1‑116 (MQ=255)

                                                                                                                                                                                             tCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATaa                                                                                                                                                                                                                                          >  7:26564/1‑110 (MQ=255)

                                                                                                                                                                                             tCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGaaa                                                                      <  2:547719/274‑1 (MQ=255)

                                                                                                                                                                                             tCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGaaa                                                                      <  2:164352/274‑1 (MQ=255)

                                                                                                                                                                                                         cacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAc                                                                                                                                                                                                                                                              <  8:1301649/78‑1 (MQ=255)

                                                                                                                                                                                                         cacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTa                                                                                                                                                                                       <  5:397241/149‑1 (MQ=255)

                                                                                                                                                                                                           cttgCCCTGTTGGTCTGCTGACCCGCCCAATCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGAAATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCg                                                         <  2:668903/270‑1 (MQ=255)

                                                                                                                                                                                                              gCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                                                                                                                                                                                             >  4:448193/1‑106 (MQ=255)

                                                                                                                                                                                                              gCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACgctgc                                                                                             >  4:339723/1‑234 (MQ=255)

                                                                                                                                                                                                               cccTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCAtt                                                                                                                                                 <  4:261219/181‑1 (MQ=255)

                                                                                                                                                                                                                 cTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACt                                                                                                                                                                               >  6:780693/1‑149 (MQ=255)

                                                                                                                                                                                                                  tGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGAc                                                                                                                                                                                                                                 >  5:1188656/1‑98 (MQ=255)

                                                                                                                                                                                                                  tGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTa                                                                                                                                                                                                                         <  2:448197/106‑1 (MQ=255)

                                                                                                                                                                                                                    ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                                                                                                                                            >  7:595489/1‑149 (MQ=255)

                                                                                                                                                                                                                    ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                                                                                                                                            >  5:1478504/1‑149 (MQ=255)

                                                                                                                                                                                                                    ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                                                                                                                                            <  8:255348/149‑1 (MQ=255)

                                                                                                                                                                                                                     tGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCTACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAGCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCa                                             <  2:97769/275‑1 (MQ=255)

                                                                                                                                                                                                                            ctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGt                                      >  2:395132/1‑275 (MQ=255)

                                                                                                                                                                                                                            ctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGAATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAg                                       <  2:30652/274‑1 (MQ=255)

                                                                                                                                                                                                                              gACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCTACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTACCTCTACACAAACTTTTGCTGCCGATCGTCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCACCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGACAAATTGGGTTATACCGTCAACAAAACAAACGAAGTAg                                    >  4:509203/1‑275 (MQ=255)

                                                                                                                                                                                                                               aCCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaaaaa                                                                                                                                                                                                                                                 <  7:518592/69‑1 (MQ=255)

                                                                                                                                                                                                                               aCCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACATAAACTTTTGCTGCCGATCt                                                                                                                                                                 >  5:1651990/1‑149 (MQ=255)

                                                                                                                                                                                                                               aCCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCt                                                                                                                                                                 >  6:278637/1‑149 (MQ=255)

                                                                                                                                                                                                                               aCCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCt                                                                                                                                                                 >  8:682673/1‑149 (MQ=255)

                                                                                                                                                                                                                                cccGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTg                                                                                                                                                                <  6:1599711/149‑1 (MQ=255)

                                                                                                                                                                                                                                 ccGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGc                                                                                                                                                                                                           >  7:93155/1‑105 (MQ=255)

                                                                                                                                                                                                                                 ccGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGc                                                                                                                                                                                                           <  8:93155/105‑1 (MQ=255)

                                                                                                                                                                                                                                 ccGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGc                                                                                                                                                               >  8:439780/1‑149 (MQ=255)

                                                                                                                                                                                                                                 ccGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCAATCTGc                                                                                                                                                               >  6:1282926/1‑149 (MQ=255)

                                                                                                                                                                                                                                  cGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGcc                                                                                                                                                              <  7:1635330/149‑1 (MQ=255)

                                                                                                                                                                                                                                  cGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGcc                                                                                                                                                              >  6:609199/1‑149 (MQ=255)

                                                                                                                                                                                                                                    cccATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGtt                                                                                                                                                                                                        >  2:721279/1‑105 (MQ=255)

                                                                                                                                                                                                                                    cccATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCgg                                                                                                                                                            <  5:761955/149‑1 (MQ=255)

                                                                                                                                                                                                                                          cattAAGTGACTCTGCACTTGCCCGGTGACGCTGGGCATTATCACCCTGCCAAAAAAAGGACTAACAATGCGACATAGCGTACTTTTTGCTACAGCGTTTGCCATGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAg                                                                                                                                                      <  7:1051926/149‑1 (MQ=255)

                                                                                                                                                                                                                                           attAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAgg                                                                                                                                                     >  5:1011593/1‑149 (MQ=255)

                                                                                                                                                                                                                                           attAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAgg                                                                                                                                                     >  5:159686/1‑149 (MQ=255)

                                                                                                                                                                                                                                               aGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTa                                                                                                                                                >  4:75784/1‑150 (MQ=255)

                                                                                                                                                                                                                                                gTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACttt                                                                                                                                                                                                                     <  5:636257/80‑1 (MQ=255)

                                                                                                                                                                                                                                                gTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTa                                                                                                                                                >  7:1568380/1‑149 (MQ=255)

                                                                                                                                                                                                                                                gTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTa                                                                                                                                                >  5:127443/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTAc                                                                                                                                               >  6:42579/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTAc                                                                                                                                               >  5:43314/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACt                                                                                                                                              >  8:204360/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    ctctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTCTTGCTGCCGATCTGCCGGGCAAANGCATTACTGTTAATCCACTTCAGCGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGCGCGCTGGAGAAANTAGGTTATACCGGAAAAAAACCCATNAAAGAAGATTACAACGTTGGcacggtctttta           >  3:170839/1‑267 (MQ=255)

                                                                                                                                                                                                                                                     tctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTtgc                                                                                                                                                                          >  7:995595/1‑118 (MQ=255)

                                                                                                                                                                                                                                                     tctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGtt                                                                                                                                           >  7:1189619/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       tGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAgg                                                                                                                                                     >  2:758952/1‑137 (MQ=255)

                                                                                                                                                                                                                                                       tGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTaa                                                                                                                                         >  7:339253/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        gCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAAt                                                                                                                                        <  8:595489/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         cACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATc                                                                                                                                       <  7:47847/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          aCTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAGCTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATcc                                                                                                                                      >  8:1519769/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          aCTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATcc                                                                                                                                      >  8:86002/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           aatGCCCGGCGTCACCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc                                                                                                                                                                                         >  7:1512199/3‑97 (MQ=255)

                                                                                                                                                                                                                                                              gCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAAAAACCTTCAAGACGCTGCTGGTCAGTCGGGCGCTGGAAAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGCTTACAACTGTGGATAACCCTCGCTTGTTTcc    >  4:560650/1‑275 (MQ=255)

                                                                                                                                                                                                                                                              gCCCGGCGTCACCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGtt                                                                                                                                  >  7:1242981/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                ccGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCgg  <  2:509207/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                  ggTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCaga                                                                                                                              <  8:517388/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   gTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCa                                                                                                                                                                                                   >  4:705273/1‑79 (MQ=255)

GGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTTCGCGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGG  >  NC_000913/2806831‑2807361

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