BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,859,354	T→C	28.8%	N255D (AAC→GAC)	yidK ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,859,354	0	T	C	28.8%	1080.0 / inf	680	N255D (AAC→GAC)	yidK	putative transporter
Reads supporting (aligned to +/- strand): ref base T (270/214); new base C (90/106); total (360/320)
Fisher's exact test for biased strand distribution p-value = 2.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 704 total aligned reads displayed.

AAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATACCG  >  NC_000913/3859205‑3859503

aaaCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGtt                                                                                                                                                       >  1:2928841/1‑150 (MQ=255)

aaaCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGtt                                                                                                                                                       >  2:3052164/1‑150 (MQ=255)

aaaCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGtt                                                                                                                                                       >  2:370888/1‑150 (MQ=255)

  aCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCa                                                                                                                                                     >  1:3344342/1‑150 (MQ=255)

   cGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCAc                                                                                                                                                    >  1:2995466/1‑150 (MQ=255)

     cAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCa                                                                                                                                                  >  2:2779889/1‑150 (MQ=255)

      aaTCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAg                                                                                                                                                 >  2:239028/1‑150 (MQ=255)

       aTCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTATTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGa                                                                                                                                                <  1:1530532/150‑1 (MQ=255)

        tCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGaa                                                                                                                                               >  2:4071652/1‑150 (MQ=255)

         cAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAAt                                                                                                                                              >  1:1389395/1‑150 (MQ=255)

            cccTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCaa                                                                                                                                           >  1:2656246/1‑150 (MQ=255)

             ccTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCaaa                                                                                                                                          >  2:1268328/1‑150 (MQ=255)

             ccTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCaaa                                                                                                                                          <  2:1820487/150‑1 (MQ=255)

             ccTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGAACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCaaa                                                                                                                                          <  1:39922/150‑1 (MQ=255)

              cTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAAc                                                                                                                                         >  2:3855310/1‑150 (MQ=255)

                   agcaCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTa                                                                                                                                    >  1:1228325/1‑150 (MQ=255)

                    gcaCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTaa                                                                                                                                   <  2:1157790/150‑1 (MQ=255)

                      aCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAAt                                                                                                                                 >  2:3216803/1‑150 (MQ=255)

                      aCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAAt                                                                                                                                 >  2:2856828/1‑150 (MQ=255)

                        cAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGc                                                                                                                               >  1:3505949/1‑150 (MQ=255)

                             ccAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCgcgc                                                                                                                          <  1:3872616/150‑1 (MQ=255)

                                gCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCaa                                                                                                                       >  2:1853073/1‑150 (MQ=255)

                                   ggTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCg                                                                                                                    >  2:2266190/1‑150 (MQ=255)

                                   ggTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCg                                                                                                                    >  1:482354/1‑150 (MQ=255)

                                       gAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCaa                                                                                                                <  1:419685/150‑1 (MQ=255)

                                       gAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCaa                                                                                                                <  1:590972/150‑1 (MQ=255)

                                       gAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCaa                                                                                                                >  2:2695756/1‑150 (MQ=255)

                                       gAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGAATTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCaa                                                                                                                <  2:1025192/150‑1 (MQ=255)

                                        aGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAg                                                                                                               >  2:715546/1‑150 (MQ=255)

                                         gCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAgg                                                                                                              <  1:2200367/150‑1 (MQ=255)

                                            ttttCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAt                                                                                                           <  1:3802269/150‑1 (MQ=255)

                                                cAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGGTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTCAATGCCGCGCCAATCGGCAAGGGAtcggt                                                                                                       >  2:2451654/1‑150 (MQ=255)

                                                 aGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGGTGCACGATGCCCTGATTTGTACACCCGGAAAAGGTGTCCACCAGAATCAAACCGGGAAATGCCGCGCCCATCGGCAAGGGCtcgggc                                                                                                      >  1:3564456/1‑148 (MQ=255)

                                                   cACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggtcgg                                                                                                    <  1:3139806/150‑1 (MQ=255)

                                                   cACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggtcgg                                                                                                    <  1:1319050/150‑1 (MQ=255)

                                                    aCCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCggg                                                                                                   <  1:960524/150‑1 (MQ=255)

                                                    aCCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCggg                                                                                                   >  1:2307873/1‑150 (MQ=255)

                                                    aCCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCggg                                                                                                   <  2:905078/150‑1 (MQ=255)

                                                    aCCGCCGTTAACAGCGCCCCCATTTGCCCATCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCAGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCggg                                                                                                   >  1:537726/1‑150 (MQ=255)

                                                     ccgccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGc                                                                                                  >  2:478355/1‑150 (MQ=255)

                                                      cgccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTAGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGcc                                                                                                 >  1:99297/1‑150 (MQ=255)

                                                          gTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCg                                                                                             >  2:3286803/1‑150 (MQ=255)

                                                           ttAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGa                                                                                            >  2:3383756/1‑150 (MQ=255)

                                                               cAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGa                                                                                        <  2:1928969/150‑1 (MQ=255)

                                                                  cgcCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGtt                                                                                     >  1:1978763/1‑150 (MQ=255)

                                                                    cccccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTa                                                                                   >  2:3815911/1‑150 (MQ=255)

                                                                    cccccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTa                                                                                   >  1:3455330/1‑150 (MQ=255)

                                                                    cccccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTa                                                                                   >  1:414580/1‑150 (MQ=255)

                                                                      cccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAAt                                                                                 <  2:2050527/150‑1 (MQ=255)

                                                                      cccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAAt                                                                                 >  2:1722642/1‑150 (MQ=255)

                                                                       ccTATTGCCCATAAGCCAGGCTTATTGAAGCCAGCGATCGTTGCACCATCCCCTTATTTTTAAACCAGTAAAACGTGTTCACCACAATCAAACAGGAAAAAGACGCGCCAATCGGCAAGGGATAGGTCGGGCCACCCAAAGAGATTAAtt                                                                                >  2:3122248/1‑150 (MQ=255)

                                                                          tttGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTtct                                                                             <  1:1009756/150‑1 (MQ=255)

                                                                           ttGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTtctc                                                                            <  2:2615436/150‑1 (MQ=255)

                                                                               ccTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCg                                                                        >  2:2945229/1‑150 (MQ=255)

                                                                                cTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGt                                                                       >  1:2063761/1‑150 (MQ=255)

                                                                                  tCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGa                                                                     >  2:1371737/1‑150 (MQ=255)

                                                                                        aGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtg                                                               <  1:1604580/150‑1 (MQ=255)

                                                                                          gCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtggt                                                             >  2:1067215/1‑150 (MQ=255)

                                                                                          gCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtggt                                                             >  2:2567893/1‑150 (MQ=255)

                                                                                            tttttGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGa                                                           <  2:3481183/150‑1 (MQ=255)

                                                                                               ttGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGtt                                                        <  1:2695756/150‑1 (MQ=255)

                                                                                                  aCGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCt                                                     >  1:2771623/1‑150 (MQ=255)

                                                                                                      cAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAAt                                                 <  2:2132349/150‑1 (MQ=255)

                                                                                                         cgtgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                              <  2:334895/150‑1 (MQ=255)

                                                                                                         cgtgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                              <  1:3721161/150‑1 (MQ=255)

                                                                                                         cgtgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                              <  2:1570668/150‑1 (MQ=255)

                                                                                                          gtgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGccc                                             <  1:253768/150‑1 (MQ=255)

                                                                                                           tgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCt                                            <  1:3074239/150‑1 (MQ=255)

                                                                                                           tgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCt                                            <  1:329326/150‑1 (MQ=255)

                                                                                                            gcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTg                                           <  2:575767/150‑1 (MQ=255)

                                                                                                            gcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTg                                           >  2:3549572/1‑150 (MQ=255)

                                                                                                              gtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                                         <  2:3621800/150‑1 (MQ=255)

                                                                                                               ttGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCAt                                        <  1:468201/150‑1 (MQ=255)

                                                                                                                     gATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCt                                  <  1:668533/150‑1 (MQ=255)

                                                                                                                     gATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCt                                  <  1:1762693/150‑1 (MQ=255)

                                                                                                                       tGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGc                                >  2:2222235/1‑150 (MQ=255)

                                                                                                                         cccTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGccc                              <  2:1671777/150‑1 (MQ=255)

                                                                                                                          ccTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCt                             >  1:2860153/1‑150 (MQ=255)

                                                                                                                          ccTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCt                             >  1:2427726/1‑150 (MQ=255)

                                                                                                                           cTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCtt                            >  2:816868/1‑150 (MQ=255)

                                                                                                                            tGATTTGTACACCAGTAAAAGGTGTTCACCAGACTCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTg                           <  1:1677882/150‑1 (MQ=255)

                                                                                                                               tttGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGccc                        >  2:3987323/1‑150 (MQ=255)

                                                                                                                               tttGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGccc                        <  2:146769/150‑1 (MQ=255)

                                                                                                                                  gTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATc                     <  1:3358014/150‑1 (MQ=255)

                                                                                                                                   tACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATcg                    <  1:4039994/150‑1 (MQ=255)

                                                                                                                                   tACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATcg                    >  2:4154361/1‑150 (MQ=255)

                                                                                                                                    acacCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATcgc                   <  1:1067642/150‑1 (MQ=255)

                                                                                                                                      acCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGa                 >  1:1147578/1‑150 (MQ=255)

                                                                                                                                      acCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGa                 >  1:194682/1‑150 (MQ=255)

                                                                                                                                        cAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGAtt               >  1:2499368/1‑150 (MQ=255)

                                                                                                                                         aGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTa              <  1:787719/150‑1 (MQ=255)

                                                                                                                                          gTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAg             <  1:4252898/150‑1 (MQ=255)

                                                                                                                                            aaaaGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGc           >  2:3939478/1‑150 (MQ=255)

                                                                                                                                             aaaGGTGTCCACCAGAATCAAACCGGTAAATGCTGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGcc          <  1:3413532/150‑1 (MQ=255)

                                                                                                                                             aaaGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGa                 >  2:4252569/1‑143 (MQ=255)

                                                                                                                                              aaGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCa         >  2:653566/1‑150 (MQ=255)

                                                                                                                                                   gtTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATAc    <  1:3203066/150‑1 (MQ=255)

                                                                                                                                                    ttCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATAcc   >  1:1850358/1‑150 (MQ=255)

                                                                                                                                                    ttCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATAcc   >  1:3650254/1‑150 (MQ=255)

                                                                                                                                                     tCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATACCg  <  2:2332884/150‑1 (MQ=255)

AAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATACCG  >  NC_000913/3859205‑3859503

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