BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,109,551	C→T	P240S (CCA→TCA)	sbp →	sulfate transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,109,551	0	C	T	100.0%	440.8 / NA	111	P240S (CCA→TCA)	sbp	sulfate transporter subunit
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (55/56); total (55/56)

Only 100 of 111 total aligned reads displayed.

GGCGCGCGTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAA  >  NC_000913/4109404‑4109696

ggCGCGCGTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCa                                                                                                                                                 <  3:1282938/150‑1 (MQ=255)

ggCGCGCGTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCa                                                                                                                                                 <  2:994171/150‑1 (MQ=255)

    cgcgTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCgt                                                                                                                                            >  3:1176181/1‑151 (MQ=255)

     gcgTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCgtg                                                                                                                                           <  1:339522/151‑1 (MQ=255)

      cgTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCgtgt                                                                                                                                          <  2:870585/151‑1 (MQ=255)

           cTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCggtg                                                                                                                                     >  3:77691/1‑151 (MQ=255)

              cACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTc                                                                                                                                  >  3:760897/1‑151 (MQ=255)

               aCTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCg                                                                                                                                 <  3:449110/151‑1 (MQ=255)

                cTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCg                                                                                                                                 >  2:1331468/1‑150 (MQ=255)

                  aaCACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGa                                                                                                                                <  3:1159221/149‑1 (MQ=255)

                  aaCACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGAt                                                                                                                               <  1:429763/150‑1 (MQ=255)

                   acacTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATaa                                                                                                                             <  3:1189066/151‑1 (MQ=255)

                   acacTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATaa                                                                                                                             <  2:522584/151‑1 (MQ=255)

                   acacTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATa                                                                                                                              >  2:327828/1‑150 (MQ=255)

                   acacTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCAGTGGTCGATa                                                                                                                              >  4:183528/1‑150 (MQ=255)

                    cacTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATaa                                                                                                                             >  2:690503/1‑150 (MQ=255)

                        tttGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTg                                                                                                                         <  2:191165/150‑1 (MQ=255)

                          tGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGt                                                                                                                       <  3:645699/150‑1 (MQ=255)

                           gTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTc                                                                                                                      <  3:807158/150‑1 (MQ=255)

                              gAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCg                                                                                                                     >  3:1103179/1‑148 (MQ=255)

                               aGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaa                                                                                                                 <  2:347361/151‑1 (MQ=255)

                               aGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaa                                                                                                                 <  2:635133/151‑1 (MQ=255)

                                 cgcgGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaaa                                                                                                                >  3:363090/1‑150 (MQ=255)

                                  gcgGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaaaa                                                                                                               >  3:619078/1‑150 (MQ=255)

                                  gcgGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaaaa                                                                                                               >  1:448020/1‑150 (MQ=255)

                                      aaTTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAgg                                                                                                           >  3:1364870/1‑150 (MQ=255)

                                      aaTTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGt                                                                                                          <  2:245736/151‑1 (MQ=255)

                                      aaTTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGt                                                                                                          <  1:303557/151‑1 (MQ=255)

                                       aTTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGt                                                                                                          >  4:249730/1‑150 (MQ=255)

                                         tGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTAc                                                                                                        <  3:1051024/150‑1 (MQ=255)

                                          ggCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTa                                                                                                      <  3:10156/151‑1 (MQ=255)

                                          ggCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTa                                                                                                      <  1:145518/151‑1 (MQ=255)

                                              aTGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAaga                                                                                                  <  4:141207/151‑1 (MQ=255)

                                               tGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAag                                                                                                   >  2:140516/1‑149 (MQ=255)

                                                      aTTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCgg                                                                                          >  3:1119840/1‑151 (MQ=255)

                                                       ttGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCgg                                                                                          >  3:261847/1‑150 (MQ=255)

                                                        tGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGa                                                                                         >  3:428963/1‑150 (MQ=255)

                                                         gCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGaa                                                                                        <  2:162997/150‑1 (MQ=255)

                                                         gCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGa                                                                                         <  2:289492/149‑1 (MQ=255)

                                                         gCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGa                                                                                         >  1:336636/1‑149 (MQ=255)

                                                          ccTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGc                                                                                      >  2:1401252/1‑151 (MQ=255)

                                                             gggAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCTGAAGCCt                                                                                    <  3:35965/150‑1 (MQ=255)

                                                             gggAAAACGAAGCCCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGTGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCt                                                                                    <  2:908051/150‑1 (MQ=255)

                                                              ggAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTAc                                                                                  >  2:178157/1‑151 (MQ=255)

                                                                aaaaCGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCt                                                                                >  2:392326/1‑151 (MQ=255)

                                                                  aaCGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTg                                                                               >  3:1222067/1‑150 (MQ=255)

                                                                  aaCGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTg                                                                               >  2:1355082/1‑150 (MQ=255)

                                                                   aCGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGaa                                                                             <  2:619051/151‑1 (MQ=255)

                                                                   aCGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGa                                                                              >  2:388876/1‑150 (MQ=255)

                                                                     gAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAat                                                                           <  3:434141/151‑1 (MQ=255)

                                                                      aaGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAat                                                                           >  3:110111/1‑150 (MQ=255)

                                                                            tgctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATctct                                                                     >  2:745697/1‑150 (MQ=255)

                                                                            tgctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATctct                                                                     >  2:1363659/1‑150 (MQ=255)

                                                                             gctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATAtctc                                                                      <  4:21851/148‑1 (MQ=255)

                                                                             gctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATctct                                                                     <  3:431551/149‑1 (MQ=255)

                                                                              ctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACt                                                                  >  3:1157171/1‑151 (MQ=255)

                                                                                    gCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAg                                                            <  2:668534/151‑1 (MQ=255)

                                                                                      gAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGaa                                                          <  3:185114/151‑1 (MQ=255)

                                                                                      gAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGa                                                           <  2:540859/150‑1 (MQ=255)

                                                                                        aTGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAgg                                                        <  3:367409/151‑1 (MQ=255)

                                                                                        aTGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAgg                                                        <  2:1201190/151‑1 (MQ=255)

                                                                                          gAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGt                                                       <  2:155240/150‑1 (MQ=255)

                                                                                          gAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGGGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGGGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTc                                                      >  2:437284/1‑151 (MQ=255)

                                                                                           aaCTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTc                                                      >  2:171464/1‑150 (MQ=255)

                                                                                           aaCTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCa                                                     <  2:841626/151‑1 (MQ=255)

                                                                                            aCTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAg                                                    <  3:815640/151‑1 (MQ=255)

                                                                                               ggggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAgg                                                   <  2:1176125/149‑1 (MQ=255)

                                                                                               ggggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGaa                                                 <  3:1157007/151‑1 (MQ=255)

                                                                                                   aaaGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTg                                             <  2:1364801/151‑1 (MQ=255)

                                                                                                    aaGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGc                                            >  3:1387802/1‑151 (MQ=255)

                                                                                                    aaGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGc                                            >  3:746229/1‑151 (MQ=255)

                                                                                                      gATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCcg                                          <  3:327857/151‑1 (MQ=255)

                                                                                                        tAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCcgcg                                        >  2:1174317/1‑151 (MQ=255)

                                                                                                          aaTTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGa                                       >  3:954636/1‑150 (MQ=255)

                                                                                                              cGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGaaaaa                                   <  3:107663/150‑1 (MQ=255)

                                                                                                              cGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAac                                  >  2:1257888/1‑151 (MQ=255)

                                                                                                                    cGTCACGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACtactac                            >  2:185303/1‑151 (MQ=255)

                                                                                                                        aCGCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTc                        >  4:371042/1‑151 (MQ=255)

                                                                                                                          gCCGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCcg                      <  2:753626/151‑1 (MQ=255)

                                                                                                                           ccGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCcgc                     <  2:1174476/151‑1 (MQ=255)

                                                                                                                           ccGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCcgc                     >  2:1403622/1‑151 (MQ=255)

                                                                                                                           ccGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCcgc                     >  2:646240/1‑151 (MQ=255)

                                                                                                                            cGAGTGAGTCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCcgcg                    >  3:12260/1‑151 (MQ=255)

                                                                                                                                 gagtCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGAc                <  3:1331541/150‑1 (MQ=255)

                                                                                                                                  agtCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGc              >  1:417911/1‑151 (MQ=255)

                                                                                                                                  agtCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGc              >  3:921968/1‑151 (MQ=255)

                                                                                                                                    tCTATCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTg            >  3:678268/1‑151 (MQ=255)

                                                                                                                                       aTCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGa           <  3:690528/149‑1 (MQ=255)

                                                                                                                                       aTCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGa           <  3:1401323/149‑1 (MQ=255)

                                                                                                                                       aTCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAgg         >  1:47327/1‑151 (MQ=255)

                                                                                                                                       aTCCTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAgg         <  4:61122/151‑1 (MQ=255)

                                                                                                                                          cTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGt        <  4:224393/149‑1 (MQ=255)

                                                                                                                                          cTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTgg      >  2:265784/1‑151 (MQ=255)

                                                                                                                                          cTCGCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTAAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTgg      >  3:177133/1‑151 (MQ=255)

                                                                                                                                             gCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCcgcgc                   >  2:1102147/1‑135 (MQ=255)

                                                                                                                                             gCAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCcgcgc                   <  3:1102197/135‑1 (MQ=255)

                                                                                                                                              cAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGaa  <  2:428932/151‑1 (MQ=255)

                                                                                                                                              cAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGa   >  2:109519/1‑150 (MQ=255)

                                                                                                                                              cAGAGTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGa   >  2:655079/1‑150 (MQ=255)

                                                                                                                                               agagTCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGaa  >  2:208675/1‑150 (MQ=255)

GGCGCGCGTGGCTCCACTAACACTTTTGTCGAGCGCGGAATTGGCGATGTACTGATTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAA  >  NC_000913/4109404‑4109696

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