BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	129,802	T→G	12.2%	K487Q (AAG→CAG)	yacH ←	DUF3300 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	129,802	0	T	G	12.2%	237.3 / 5.3	131	K487Q (AAG→CAG)	yacH	DUF3300 family protein
Reads supporting (aligned to +/- strand): ref base T (61/50); new base G (7/12); total (68/63)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

GTGATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAGT  >  NC_000913/129655‑129945

ttGATGCGGACTGAATACGCTCTCGCGCGGCGGCTCTTCGCTGCTGTTTCTGCGGGTTAGGCTGCATTTTATCGCGAAACTCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTc                                                                                                                                                <  1:120866/148‑1 (MQ=25)

gtgtTGCGGACTGACGACGCGCCCGCGCGGCGGCTCTTCGCTGAGGGTTCTGCGGGATCGGCGGCCTTTTCTCGCGCCCCGCCTGGCGCGGCTCAGGCGAGGCGGCCTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTc                                                                                                                                                <  2:748192/149‑1 (MQ=255)

 gtaggCGGACTGAATACGCTCCCGCTCTGCGTCTCGTCGCGGCTGCTTCTGCGGGCTCGGCTGCACTTTCTCCCGAACCTCCTGGCGCTGCTCAGGCGAGGCGGACTGAATTCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCg                                                                                                                                               <  2:295229/145‑1 (MQ=21)

     gCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCTCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACAGAATACGCTCACGAGCTGCCTCACTTCGATGCTGGATCTTCGGGtt                                                                                                                                           >  2:620290/1‑149 (MQ=18)

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

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

        gACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCGTCGATTAGGTGCTGGATCTTAGGGTTAgg                                                                                                                                        >  2:207014/1‑148 (MQ=18)

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

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

              aTACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAAAGCGCTCACGCGCGGCGTCGGTTCGGTGCTGGATCTTCGGTTTAGTcgggct                                                                                                                                  >  2:500694/1‑144 (MQ=18)

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

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

                  gCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATCCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATtctc                                                                                                                              <  1:1293252/149‑1 (MQ=255)

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

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

                    tCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGAGCTGCCTCACTTCGCTGCTGGATCTTCGGGTGAGTCTGCAttaagtc                                                                                                                            >  2:789214/1‑144 (MQ=255)

                    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)

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

                                            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)

                                                      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)

                                                              tGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACg                                                                                  >  2:669965/1‑149 (MQ=255)

                                                              tGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACg                                                                                  >  1:817654/1‑149 (MQ=255)

                                                              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)

                                                                 aTTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCCGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTg                                                                               >  2:578084/1‑149 (MQ=255)

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

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

                                                                        cgcgAACCTCCTGGCGCTGCTCAGGCGATGCGGACTGAATCCGCTCACGCGCTGCCTCTCTTCGCTGCTGTATCTTCGGGTTAGTCTGCATTATCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAgcg                                                                        <  2:253823/149‑1 (MQ=18)

                                                                         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)

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

                                                                                         tgccccggcgatgcgGCCTGAATTCGCTCACGCGCTGCCTCTCTTCGCTGCTGTATCTGCGGTTTAGGCTGCATGATCTCGCTAACTGCCTGGCGCTGCTCAGGCGAGTCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTg                                                       <  1:1399528/143‑1 (MQ=255)

                                                                                          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)

                                                                                                       ggACTGAATTCGCTCCCGCGCTGCATCTCTTTGCTGTTGCTTCTTCGGGTTAGTCTGTATTTTTTCGCTAACTTCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGc                                         <  2:355721/149‑1 (MQ=21)

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

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

                                                                                                          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                                                 <  1:1076389/133‑1 (MQ=255)

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

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

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

                                                                                                                      aCGCGCTGCCTCTCTTCGCTGCTGGCTGTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGCCTGATGACGCTGCCGAGCGGCGTCCTTTTGTGGCTGGGTGAGTGGTGGGCGACGGCGGAAGTCgg                          >  2:1212267/1‑148 (MQ=21)

                                                                                                                      aCGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGCCGGCTGAAGTCGt                          >  1:684239/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)

                                                                                                                       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)

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

                                                                                                                                    ttCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAg            >  1:472482/1‑149 (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)

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

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

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

GTGATGCGGACTGAATACGCTCACGCGCTGCGTCTCTTCGCTGCTGGTTCTGCGGGTTAGTCTGCATTTTCTCGCGAACCGCCTGGCGCTGCTCAGGCGAGGCGGACTGAATGCGCTCACGCGCTGCCTCTCTTCGCTGCTGGATCTTCGGGTTAGTCTGCATTCTCTCGCGAACTGCCTGGCGCTGCTCAGGCGAGGCGGACTGATAACGCTGACGAGCGGCGTCCTTTTGTTGCTGGGTCAGTGGTTGGCGACGGCTGAAGTCGTGGAAGTCGTCATAGCTCCCATAGT  >  NC_000913/129655‑129945

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: `‑`