BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	129,811	T→G	14.4%	T484P (ACT→CCT)	yacH ←	DUF3300 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	129,811	0	T	G	14.4%	211.5 / 3.5	124	T484P (ACT→CCT)	yacH	DUF3300 family protein
Reads supporting (aligned to +/- strand): ref base T (60/42); new base G (5/17); total (65/59)
Fisher's exact test for biased strand distribution p-value = 3.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAGTGTTCAGCTTCATTA  >  NC_000913/129663‑129959

gACTTATTACTCTCACGCGCGGCGTCTCTTCGCTGCTGGTTCTGCGGGGTAGGCTGGTATTTCTCGGGAACCGCCTGGCGCTGCTCAGGCGAGGCGTACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGt                                                                                                                                                      <  2:1548979/149‑1 (MQ=25)

gACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTGTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGt                                                                                                                                                      <  1:819290/149‑1 (MQ=255)

gACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGt                                                                                                                                                      <  1:826049/149‑1 (MQ=255)

  atGAATACGCTCTCGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCt                                                                                                                                                    <  1:473293/148‑1 (MQ=255)

     aaTACGCCCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCa                                                                                                                                                 <  2:787044/149‑1 (MQ=255)

          gCTCTCGCGCTGTGTCTCTTCGCTGCTGGTTCTGCGGGATAGTCTGCAATTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATtctc                                                                                                                                            <  2:1259787/149‑1 (MQ=255)

          gCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATTGCGCTCACGCGCTGCCTCAGTTAGGTGCTGGGTCTGAGGGTTAGTCGGCCTtctc                                                                                                                                            >  2:148008/1‑149 (MQ=18)

          gCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAAAGCGCTCACGCGCTGCCTCGCTTAGCTGCTGGATCGTCGGGTTTCTCTGCATgctc                                                                                                                                            >  2:1272824/1‑149 (MQ=21)

          gCTCACGCGCAGCGTCTCATCGCTGCTGGTTCAGCGCGTTAGTCTGCCTTTTCTCGCGAACAGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCATCGCTGCTGCATCTTGGGGGTACTCAGcctacta                                                                                                                                            >  2:200301/1‑143 (MQ=18)

            tCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATTGCGCTCACGCGCTGCCTCGGTTAGCTGCTGGTTCATAGGGTGTGGCTGGATtctgtc                                                                                                                                          >  2:257341/1‑149 (MQ=17)

            tCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATCCGCTCACGCGCTGCGTCTATTAGATGCTGGTTCTTACGTTTAGTCTGCATtctgtc                                                                                                                                          >  2:952812/1‑149 (MQ=18)

             cACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGAGGGATCTTCGGGATAGTCTGCATTCTCTcg                                                                                                                                         >  1:1157008/1‑149 (MQ=255)

             cACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCGTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTcg                                                                                                                                         >  1:1197619/1‑149 (MQ=255)

              aCGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCACTTCGCTGCTGGATCTGCGGGTTAGTCAGCCTTCTCTcgc                                                                                                                                        >  1:901180/1‑149 (MQ=255)

              aCGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTAGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATTGCGCTCACGCGCTGCCTCTGTTAGCTGCTGGATCTTAGGGTTAGTCTGCATTCTCTcgc                                                                                                                                        >  1:467756/1‑149 (MQ=25)

                    tGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTg                                                                                                                                  >  1:1295973/1‑149 (MQ=255)

                    tGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTAGGTGCTGGATCTGCGGGTTAGGCTGCATTCACTCGCGAACTg                                                                                                                                  >  1:1540094/1‑149 (MQ=25)

                         ctctTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTGACGCGCTGCCTCTCTTCGGTGCTGGATCAGAGGGTTAGTCTGCATGCTCTCGAGAACTGCCAgg                                                                                                                             >  1:1057441/1‑149 (MQ=21)

                               gctgctGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTAGCTTCTGGATCTTAGGGTTAGTCTGCATTCTCTCACGAAATGACTGGAtctgc                                                                                                                       >  1:311441/1‑149 (MQ=25)

                                 tgctgTTTCTGTGGCCGAGTCTGCACTTTTTCGCGAACCCCCTGGCGCTGCTCAGGCGAGGCGGACTGAATACGCTCACGCGCTGCCTCTCTTCGCTGCTGTATCTTCGGGCTAGGCTGCATTCTCTCCCGAACTGCCTGGCGCTGCTc                                                                                                                     <  2:514148/149‑1 (MQ=18)

                                 tgctgGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAAATCGCTCACGAGCTGCCTCTCTTAGCTGCTGGATCTTCGGGTTAGTAAGCATTCTCTCGAGAACTGCCTGGAGCAGCTc                                                                                                                     >  2:322269/1‑149 (MQ=25)

                                    tgGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACAGCCTGGCGCTGCTCAGGCGAGGCGGCCTGAATGAGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCGTCGGGTGAGTCTGCATTCTCTCGCGAACTGCCTGACGCTGCTCagg                                                                                                                  >  1:1156572/1‑149 (MQ=38)

                                    tgGTTCTGCGGGGTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAAGGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCagg                                                                                                                  <  1:481889/149‑1 (MQ=255)

                                     ggTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCGGTTCGCTGCTGGATCTGCGGGTTAGTCGGCATTCTCTCGCGGACTGCCTGGCGCTGCTCaggg                                                                                                                 >  2:141235/1‑148 (MQ=38)

                                      gTTCTGCGGGGGAGTCTGCAGTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCGGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCaggcg                                                                                                                <  1:305139/149‑1 (MQ=255)

                                          tGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAAAGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGTTCTTCGGGTTAGTCTGCCTTCTCTCGCGAACTGCCTGCCGCTGCTGcggcgaggc                                                                                                            >  2:1414718/1‑149 (MQ=37)

                                           gCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCaggcgaggcg                                                                                                           >  1:236297/1‑149 (MQ=255)

                                           gCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCaggcgaggcg                                                                                                           >  1:990703/1‑149 (MQ=255)

                                             gggCTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTGCGGGTTAGGCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGa                                                                                                         <  1:1265701/149‑1 (MQ=255)

                                              gggtcggcTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGGTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGAc                                                                                                        <  1:777231/142‑1 (MQ=255)

                                              ggTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCGCAGGCGAGGCGGAc                                                                                                        >  1:427897/1‑149 (MQ=255)

                                               gatAGTCTGCAGTTTCTCCCTAACCTCCTGGCGCTGCTCAGGCGATGCGGACTGAATTCGCTCACGCGCTGCCTCTCTTCTCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACt                                                                                                       <  2:808165/147‑1 (MQ=17)

                                                   ctCTGTAGTTTCTCGTTCAACGCTTGCCGATGCTCAGGGTTTGCGGACTGAATTCGCTCACGCGCTGACTCTCTTCGCTGCTGGATCTTTGGGTTAGTCTGCATTCTCTCCCTAACAGCCTGGCGCTGCTCAGGCGAGGCGGACTGATa                                                                                                   <  2:922579/148‑1 (MQ=9)

                                                      tGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATTCGCTCACGCGCTGCCTCTAATGGCTGCTGGATCTGAGGGGAAGTCAGCATTCTAGCGCGATCGTCATGGCGCTGCTCAGGCGAGGCGGGCTCATAACg                                                                                                >  2:1059446/1‑149 (MQ=18)

                                                      tGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGAATCACTTCGCTGATGGATCTTCGGGTTAGTCTGCTTTCTCTAGCGAACTGCCTGGCGCTGCTCAGCCGAGGCGGACTGATAACg                                                                                                >  2:1430652/1‑149 (MQ=38)

                                                       gCATTTTCTCGCGAACAGCCTGGCGCTGCTCAGGCGAGGCGGTCTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGATGGATATTCGGGGTAGTCTGCAGTCTCTCGCGAACTGCGTGGCGCTGCTCAGGCGTGGCGGAGTGATAACGc                                                                                               >  2:1404060/1‑149 (MQ=38)

                                                        cATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCt                                                                                              >  2:460284/1‑149 (MQ=255)

                                                        cATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTAGGCTTAGTCTGCCTTCACTCGAGGACGGCCTGGGGGAGCGCAGGCGAGGCGGACTGATAACGCt                                                                                              >  2:385033/1‑149 (MQ=37)

                                                           tttCTCCCGAACCGCCTGGCGCTGCGCAGGCGCGGCGGACTGACTGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGAc                                                                                           <  2:1197619/149‑1 (MQ=255)

                                                           ttactcGCGAACCGCCTGGCGCTGCGCAGCCGAGGCGGTCTGCCTGGGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGAc                                                                                           <  2:1157008/146‑1 (MQ=21)

                                                                 gcgAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAgcgg                                                                                     >  1:485008/1‑149 (MQ=255)

                                                                 gcgAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGGCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAgcgg                                                                                     >  1:966786/1‑149 (MQ=255)

                                                                   gAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAgcggcg                                                                                   >  2:565533/1‑149 (MQ=255)

                                                                    aaCCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGAATGATAACGCTGAGGAgcggcgg                                                                                  >  1:885468/1‑148 (MQ=255)

                                                                       cGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGCGGCGGACTGATAACGCTGCCGAGCGGCGTCCt                                                                               >  2:1268779/1‑149 (MQ=255)

                                                                         ccTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTGCGCTGCTGGCTCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCGGCGCAGGCGAGGCGGACTGATAACGCGGACGAGCGGAGTCCtgt                                                                             >  1:707599/1‑147 (MQ=255)

                                                                          cTGGCGCTGCTCAGGCGATGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTCGGCTGAATTTTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCtttt                                                                            <  2:631884/149‑1 (MQ=18)

                                                                           tGGCCCTGCTCTGGCGAGGCGGACTGAATTCGCTCTCGCGCTGCCTCTCTTCGCTGCTGGATCTGCGGGTTAGGCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTttg                                                                           <  2:716855/149‑1 (MQ=18)

                                                                           gggCGCTGCTCAGGCGAGGCGGACTGAATCCGCTCACCCGCTGCGTCTCTTCGCTGCGCGAGCTTCGGGCTAGGCAGCATTCTCTCGCTAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTttg                                                                           <  1:1255310/148‑1 (MQ=18)

                                                                                  gctctcgcgaggcgGACTTAATGCGCGCACGCGCTGCCTCGCTTCGCTGCTCGCTCTTCGGGGTAGTCTGCAGTCTCTCGCGAACACCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTgg                                                                    <  1:191458/143‑1 (MQ=21)

                                                                                   gtCAGGCGAGGCGGAGTGAATTAGCTGACGGGCTGTTTCTCTTTGCTGATGTTTCTTTGGGTTAGTCTGCATTTTATCGTTAAATTCCTGGCGCTGCTCAGGCGAGGAGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTggg                                                                   <  2:908667/148‑1 (MQ=16)

                                                                                   ctCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTggg                                                                   >  1:262426/1‑149 (MQ=255)

                                                                                   ctCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCGCAGGGGGGGCGGCGTGATAACGCTGAGGGGGGGCGTGGGTTTGtgggtggg                                                                   >  1:1571969/1‑149 (MQ=255)

                                                                                    tCAGGCGTGGCGGAGTGTATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGGTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGt                                                                  <  1:1059446/149‑1 (MQ=255)

                                                                                    tCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGt                                                                  <  1:669965/149‑1 (MQ=255)

                                                                                        gcgaggcgGACTGAATGCGCTCACGCGCTGCCTCTCTTCTCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGt                                                              <  2:120783/149‑1 (MQ=255)

                                                                                           aggcgGACTGAATCCGCTCCCCCGCGGACTCTCTTCGCGGCTGGCTCTGCGGGGGAGGCAGCAGTCGCTCGCTCACCGCCTGGCGCGGCGCCGGCGAGGCGGACTGATAACGCTGACGCGCGGCGTCCTTTTGTTGCTGGGGCAGTGGt                                                           <  2:940185/149‑1 (MQ=255)

                                                                                            ggcgGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCGGGCTCTTCGGGGTAGGCTGCATTCTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGtt                                                          <  1:323431/149‑1 (MQ=21)

                                                                                              cgGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTgg                                                        <  1:1146713/149‑1 (MQ=255)

                                                                                               ggACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCCGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGc                                                       >  2:868017/1‑149 (MQ=255)

                                                                                               gcaCTGAATGCGCTCCCGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGGTCGGCTGCATTCTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGc                                                       <  2:1002802/147‑1 (MQ=25)

                                                                                                gACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGCTCTGATAACGCTGACGAGCGGCGTCCGTGAGTTGCTGGGTCAGTGGTTGGCg                                                      >  2:591562/1‑149 (MQ=21)

                                                                                                  cTGAATCCGCTCCCGCTCTGCCTCTCTTCGCTGCTGCCTCTTCGGGTTAGTCTGCATTCTCTCCAGAACTTCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGATGCTGGGTCAGTGGTTGGCGAc                                                    <  1:1040604/149‑1 (MQ=25)

                                                                                                   tGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACg                                                   >  1:35684/1‑149 (MQ=255)

                                                                                                    gAATTCGCTCTCGGGCGGCTTCTCTTGGCTGCTGGAGCTGCGGGCTCTTGTGTATTTTCTCGTTAACTCCCTGGCTCTGCTCAGGCGAGTAGGACTTATAACTCTGACTAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTTGCGACgg                                                  <  2:1015565/149‑1 (MQ=255)

                                                                                                     aaTGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAAATGCCTGGCGATGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCATTTTGTTGCTGGGTCAGTGGTTGGCGACGGc                                                 >  2:322941/1‑149 (MQ=255)

                                                                                                     aaTGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGGTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCGCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGGCCGTGTGTTGCTGGGGCAGTGGGTGGCGACGGc                                                 >  1:523214/1‑149 (MQ=255)

                                                                                                      ctGCGCTCACGAGCTGCCTCTCTTCGCTGCGGGATCTCCGGGGTAGTCTGCATTGTCTCGCGACCTGCCTGGAGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCt                                                <  1:805716/148‑1 (MQ=25)

                                                                                                       tGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAg                                                               >  2:1076389/1‑133 (MQ=255)

                                                                                                       tGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAg                                                               <  1:1076389/133‑1 (MQ=255)

                                                                                                        cagcTCTCGCGCTGAATCTCTTCGCTGCTGCATCAGCGGGTGAGTCTGCAATATCTCGATAAATTCCTGGCGCTGCTCAGGCGAGGCGGACTTATAACGCTGACGATCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGa                                              <  2:1451425/147‑1 (MQ=21)

                                                                                                            tCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTc                                          >  1:669818/1‑149 (MQ=255)

                                                                                                              aCGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGCCGGCTGAAGTCGt                                        >  1:684239/1‑149 (MQ=255)

                                                                                                              aCGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGt                                        >  2:436500/1‑149 (MQ=255)

                                                                                                              aCGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTGTTGTTGCTGGGTCAGTGGTTGGCGACTGCTGAAGTCGt                                        >  1:1400483/1‑149 (MQ=255)

                                                                                                               gtcgcgGCCTCTCTGCGCTGCGGGATCTGCGGGTGAGGCTGCATTATCTCGCGAACAGCCGGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTg                                       <  1:717947/147‑1 (MQ=25)

                                                                                                               cgcgCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTTTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTg                                       <  2:1077143/149‑1 (MQ=255)

                                                                                                               cgcgCTGCCTCTCTTCGCTGCTGGAGCTTCGGGTTCTGCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTg                                       <  2:523214/149‑1 (MQ=255)

                                                                                                               cgcgCTGCCTCTCTTCGCTGCTGGAGCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTg                                       <  1:1373563/149‑1 (MQ=255)

                                                                                                                gcgcTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGGCTGCATTCTCTCGCGAACTGCCTGGCGCTGCGCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGGCAGTGGTTGGCGACGGCTGAAGTCGTgg                                      <  1:129901/149‑1 (MQ=255)

                                                                                                                gcgcTGCCTCTCATCGCTGCTGGATCTGCGGGTTAGTCTGCATTATCTCGCGAACGGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTgg                                      <  2:843454/149‑1 (MQ=255)

                                                                                                                    tGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAgt                                  <  1:688980/149‑1 (MQ=255)

                                                                                                                      ccTCTCTTCGCTGCTGGATCTTCGCGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTCATAACGCTGACGAGCGGCGTCCTTTAGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAgtcg                                >  1:35757/1‑149 (MQ=255)

                                                                                                                       ctctctTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAgtcgt                               <  2:577161/149‑1 (MQ=255)

                                                                                                                         ctcgTCGCTGCTGGCTCTGCGGGTTAGTCTGGATGTTCTCGCTACCAGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCa                             <  1:1356939/149‑1 (MQ=25)

                                                                                                                          tctTCGCTGCTTGATCTTCGGGGTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCAt                            <  1:620290/149‑1 (MQ=255)

                                                                                                                              cGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGGCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCt                        <  1:141235/149‑1 (MQ=255)

                                                                                                                              cGCTGCGGGATCTTCGGGGTCGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCt                        <  1:447402/149‑1 (MQ=255)

                                                                                                                                 tgctgGATCTGCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCTGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTccc                     <  1:952812/149‑1 (MQ=255)

                                                                                                                                   ctgGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCCGGCGAGGCGGACTGATAACGCTGACGCGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCAt                   >  2:88165/1‑149 (MQ=255)

                                                                                                                                   ctgGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTGAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTGATAGCTCACAt                   >  2:1540891/1‑149 (MQ=255)

                                                                                                                                   ctgGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCAt                   >  2:790649/1‑149 (MQ=255)

                                                                                                                                    tgcttCTTCGGGTTCGGCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATa                  <  2:820936/145‑1 (MQ=255)

                                                                                                                                      gATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAgt                <  2:1540094/149‑1 (MQ=255)

                                                                                                                                      gATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGAGCTGCTCAGGCGAGGCGGACTGATAACGCTGCCGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGCCGCCTGAAGTCGTGGAAGTCGTCATAGCTCCCCTAgt                >  2:657836/1‑149 (MQ=38)

                                                                                                                                            cGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAGTGTTCAg          <  1:148008/149‑1 (MQ=255)

                                                                                                                                                 tAGTATGAATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGGTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAGTGTTCAGCTTCa     <  1:789214/149‑1 (MQ=255)

                                                                                                                                                    tCTGCATTCTCTCGCGCACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATCACGCTGACGAGCGGCGTCCTGTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGGCGTGGAAGTCGTCATAGCTCCCATAGTGTTCAGCTTCATTa  >  2:690334/1‑149 (MQ=255)

                                                                                                                                                    tCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAGTGTTCAGCTTCATTa  >  2:1340473/1‑149 (MQ=255)

                                                                                                                                                    tCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGAGTCCTTTTGTTGCTCGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTGATAGATCCCATTGTGTTCAGCtgagtta  >  2:1030199/1‑143 (MQ=38)

GACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAGTGTTCAGCTTCATTA  >  NC_000913/129663‑129959

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