BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,449	C→T	6.5%	intergenic (+17/+97)	mdfA → / ← ampG	multidrug MFS transporter/muropeptide permease AmpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,449	0	C	T	6.5%	620.4 / 22.3	203	intergenic (+17/+97)	mdfA/ampG	multidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand): ref base C (32/156); new base T (5/8); total (37/166)
Fisher's exact test for biased strand distribution p-value = 6.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 216 total aligned reads displayed.

CTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGC  >  NC_002947/1549301‑1549585

gtgGCGATGATTGGCGCGGGTGCCAGCTTGCTGCTTTTTGCCCGGGGGGGGGGGGTGGGGGGGGCGTTGGCGCTGTGGCCGGTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGggcca                                                                                                                                             <  2:5271273/149‑2 (MQ=255)

cgggCGATGTTGGGCGCCGGTGCCAGCTTGCTGCATTTGGCCGGGGGGGGGGGGGGGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGggccg                                                                                                                                             <  2:1994252/148‑1 (MQ=255)

  ggCGTTGTTGGGCGCCGGTGCCCGCTTGCTGATTTTGGCCCTGGCGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGggccggc                                                                                                                                           <  1:6171823/150‑1 (MQ=255)

  ggCGATTATGGGCGCCGGTGCCAGCTTGCTGCTTTTTGCCCGGGCGGGGGGGGTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGggccggc                                                                                                                                           <  2:2075465/150‑1 (MQ=255)

     gttgttgggCGCCGGTGCCAGCTTGCTGCTTTTTGCCCTGGGGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCctc                                                                                                                                        <  1:6442549/143‑1 (MQ=255)

     gatgttGGGCGCGGGTGCCTGGTTGCTGCATTTTGCCCGGGGGGGGGGGGGGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCctc                                                                                                                                        <  1:2960527/150‑1 (MQ=255)

      ttgttgggCGCCGGTGCCAGGTTGCTGCATTTGGCCCTGGGGGGGGGGGGGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCctct                                                                                                                                       <  2:4616905/144‑1 (MQ=255)

      ttgatGGCCGCCGGTGCCAGCTGGCGCCATTTGGCCCGGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCGGTGGTTGGGGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGccccggcaata                                                                                                                                       <  1:3087268/149‑5 (MQ=255)

      ctgatTGGCGCGGGGCCCAGTTTGCTGCATTTGGCCCGGGGGGGGGGGGGGGCGGGGGCGTGGGCGGTGGGGCCGCTGTGGTTTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGccccggcaata                                                                                                                                       <  1:7084030/149‑5 (MQ=255)

        gtttGGCGCCGGTGCCCGCTTGCGGCATTTTGCCCGGGGGGGGGGGGTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTc                                                                                                                                     <  2:7363713/148‑1 (MQ=255)

           tGGCGCCGGTGCCTGCTTGCTGCTTTTTGCCGGGGCGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTcg‑c‑‑g                                                                                                                               <  2:4718607/150‑1 (MQ=255)

            ggggcCGGTGCCAGCTTGCTGCATTTGGCCGTGGCGGGGGGGGGGGGGGGGGCGTGGGCGCGGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑gg                                                                                                                              <  1:4052258/147‑1 (MQ=255)

            ggCGCCGGGGCCCGCTGGCTGCATTGGGCCCTGGCGGGGGGGGGGGGGGGGGCGTTGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑gg                                                                                                                              <  2:558833/150‑1 (MQ=255)

             gcgcCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑ggg                                                                                                                             >  2:3055806/1‑150 (MQ=255)

             gcgcCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑ggg                                                                                                                             <  2:1857891/150‑1 (MQ=255)

                 cGGTGCCAGCTTGCTGAATTTTGCCCTGGGGGGGGGGGTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTaaa                                                                                                                         <  1:116719/150‑1 (MQ=255)

                  ggTGCCAGCTTGCTGAATTTTGCCGTGGGGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAAc                                                                                                                        <  1:6029032/150‑1 (MQ=255)

                   gcgcacgCTTGGGGCATTTGGCCGGGGCGGGGGGGGGGGCGGGGGCGTTGGGGCTGTGGCCGCGGGGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAAcc                                                                                                                       <  1:368754/144‑1 (MQ=255)

                    tGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAAacc                                                                                                                      >  2:5712472/1‑150 (MQ=255)

                     gCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCg                                                                                                                     <  2:2435966/150‑1 (MQ=255)

                        aGCTTGCTGAATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTc                                                                                                                  <  1:5270391/150‑1 (MQ=255)

                        aGCTTGCTGAATTTGGCCTTGGGGGGGGGTGTGGCGGGGGCGTTGGCGCTTTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTc                                                                                                                  <  1:6164688/150‑1 (MQ=255)

                        aGCTGGCTGCATTTTGCCGGGGCGGGGGGGGGGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTc                                                                                                                  <  2:1054093/150‑1 (MQ=255)

                         ggttGCTGCATTTGGCCGGGGCGGGGGGGGGGGCGGGGGCGTGGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTcc                                                                                                                 <  1:6311331/148‑1 (MQ=255)

                           ttGGTGCTTTTTGCCCTGGGGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTa                                                                                                               <  2:3944859/150‑1 (MQ=255)

                           ttGCTGCTTTTGGCGCGGGCGGGGGGGGGGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTa                                                                                                               <  2:3253378/150‑1 (MQ=255)

                           ttGCTGCATTTTGCCCTGGGGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTa                                                                                                               <  2:1457347/150‑1 (MQ=255)

                           ttGCTGCATTTGGCCCTGGGGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTa                                                                                                               <  2:495701/150‑1 (MQ=255)

                           ttGCTGATTTTTGCCCTGGGGGGGGGGGGGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTa                                                                                                               <  2:3499582/150‑1 (MQ=255)

                             gctgcATTTTGCCCTGGCGGCGGGGGTGGCGGGGCCGTTGGCGCTGTGCCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTaca                                                                                                             <  2:2241465/150‑1 (MQ=255)

                              ctgcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGGCTCTTCG‑C‑‑GGGGAAACACGGTccaagac                                                                                                            >  2:4979929/1‑145 (MQ=255)

                              ctgcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGGCTCTACG‑CCGGGGAAAACCCggccaaaa                                                                                                              >  1:1624428/1‑143 (MQ=255)

                                gcATTTTGCCCTGGCGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAg                                                                                                          <  1:4966076/150‑1 (MQ=255)

                                gcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCACTTCG‑C‑‑GGGTAAAACCGCTCCAACACAg                                                                                                          >  2:4084898/1‑150 (MQ=255)

                                gcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCACTTAG‑C‑‑CGGTAAAACCGCACCAACAAAg                                                                                                          >  2:2136346/1‑150 (MQ=255)

                                       gCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGgccgcg                                                                                                   <  2:6069790/150‑1 (MQ=255)

                                        cccGGGGGGCGGGGGGGGGGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGgccgcgc                                                                                                  <  1:965146/150‑1 (MQ=39)

                                            ggcgggGGGGGTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCt                                                                                              <  1:5824764/150‑1 (MQ=255)

                                            ggcgggGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCggcgcggccct                                                                                              <  2:4477053/150‑6 (MQ=38)

                                             gggggggggggggggggggggTGGGCGGTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCtt                                                                                             <  1:6916272/145‑1 (MQ=39)

                                              cgggggGGGGGGCGGGGGCGTTGGCGCTGTGGCCGCTGGGGTGTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCttt                                                                                            <  2:7171520/150‑1 (MQ=255)

                                                ggggggggggCGGGGGGGTGGGCGCTGTGGCCGCTGTGTTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTgtgt                                                                                          <  1:12532/148‑1 (MQ=25)

                                                gcggGGGGGGCGGGGGCGTGGGCGCTGTGGCCGCGGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                          <  2:86621/150‑1 (MQ=255)

                                                 gggggggggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                         <  1:4288322/149‑1 (MQ=37)

                                                 gggggggggggggggCGTTGGCGCTGTGGCCGCTGGGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                         <  2:5136686/149‑1 (MQ=255)

                                                 gggggggggCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                         <  2:1723310/149‑1 (MQ=37)

                                                 gggggggggCGGGGGCGTGGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                         <  1:239423/149‑1 (MQ=255)

                                                 ggggggTGGGGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                         <  1:2204246/149‑1 (MQ=37)

                                                 ggggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                         <  2:2490185/149‑1 (MQ=255)

                                                  ggggggggggggggCGTGGGCGCGGTGGCCGCTGGGGTTTGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAg                                                                                        <  1:487069/150‑1 (MQ=25)

                                                  gggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                        <  2:3517699/150‑1 (MQ=255)

                                                  gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCTGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAg                                                                                        >  1:6966330/1‑150 (MQ=37)

                                                  gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCTGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAg                                                                                        >  1:3028999/1‑150 (MQ=37)

                                                  gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGGGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAg                                                                                        <  2:1324463/150‑1 (MQ=37)

                                                  gggggTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTATGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                        <  2:5308026/150‑1 (MQ=255)

                                                   ggggggggggggGCGTTGGCGCTGGGGCCGCGGTGGTTTGGGGTGGGTGGGCGGCCTGGGGACGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAgg                                                                                       <  1:3408600/150‑1 (MQ=25)

                                                   ggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAgg                                                                                       <  1:2594221/150‑1 (MQ=37)

                                                    gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                      <  2:2944042/150‑1 (MQ=255)

                                                    gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                      <  2:1238779/150‑1 (MQ=255)

                                                     ggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAg                                                                                     <  2:1572703/150‑1 (MQ=255)

                                                      ggggggggggCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                    <  1:3173454/145‑1 (MQ=255)

                                                       gggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGC‑g                                                                                  <  1:1344144/149‑1 (MQ=37)

                                                        gggggggggGTGGGGGCTGGGGCCGCTGTGGTTTGGGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑gg                                                                                 <  1:1056391/147‑1 (MQ=38)

                                                        ggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGC‑gg                                                                                 <  2:774774/147‑1 (MQ=25)

                                                        ggggggggCGTGGGCGCGGTGGCCGCTGTGGTTTGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGC‑gg                                                                                 <  1:530795/147‑1 (MQ=25)

                                                        ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑gg                                                                                 <  1:2008613/150‑1 (MQ=255)

                                                        ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑gg                                                                                 <  1:3463414/150‑1 (MQ=255)

                                                        ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGC‑gg                                                                                 <  2:6339096/150‑1 (MQ=37)

                                                        ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAAGCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑gg                                                                                 <  2:2063659/150‑1 (MQ=255)

                                                        ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGGAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑gg                                                                                 <  2:1682543/150‑1 (MQ=255)

                                                        ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTGGTGGTGGTTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGC‑gg                                                                                 <  1:7039977/150‑1 (MQ=25)

                                                        ggCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGC‑gg                                                                                 <  2:3922236/150‑1 (MQ=25)

                                                         gggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGC‑ggg                                                                                <  2:4804691/148‑1 (MQ=25)

                                                         gggggggCGTGGGCGCTGTGGCCGCTGGGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGC‑ggg                                                                                <  1:1052293/148‑1 (MQ=25)

                                                         gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑ggg                                                                                <  2:1788323/150‑1 (MQ=255)

                                                           gggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑GGGtt                                                                              <  1:464347/150‑1 (MQ=255)

                                                           gggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCggcgcggccct                                                                                              <  2:1540707/135‑6 (MQ=37)

                                                           gggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCggcgcggccct                                                                                              >  1:1540707/1‑130 (MQ=37)

                                                            ggggCGTTGGCGCTGTGGCCGCTGTGTTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTAGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGC‑GGGttt                                                                             <  2:2023328/150‑1 (MQ=25)

                                                            ggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTTAAACCCCTCCGACACCGGGCCGCCCCGGGGGGTAGGGGC‑GGGttt                                                                             >  1:5509198/1‑150 (MQ=255)

                                                                 gTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAAACCGCTCCTACACAGGGCCGCGCCGCGTTGTAGGAGC‑GGGGTTAcccc                                                                        >  2:6097600/1‑149 (MQ=255)

                                                                  ttGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTAG‑C‑‑GGGTAAAACCCCTCCGACACACGGCCCCCCCGCGTTGTTTGGGC‑GGGGTTAacccc                                                                       >  2:6447710/1‑148 (MQ=255)

                                                                  ttGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCGCATAGCC‑‑GGGTAAAACC‑CCCCCACACACCGCCCCCCCCCCGTGTGGGGGCGGGGTTTAcccc                                                                        >  1:2920362/1‑149 (MQ=255)

                                                                     gcgcTGGGGCCGCTGTGGTTTGGGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATACTCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑GGGTTTACCCGCGaa                                                                    <  2:6740500/150‑1 (MQ=255)

                                                                         tgtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                         >  2:4112916/1‑126 (MQ=25)

                                                                             gCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                         <  1:5752474/122‑1 (MQ=25)

                                                                                      gTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACACGCaa                                                   >  2:3621667/1‑150 (MQ=255)

                                                                                             gtgggtggGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGc                                                                                    <  1:5493064/111‑1 (MQ=14)

                                                                                             gtgggtggGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGc                                                                                    >  2:5493503/1‑111 (MQ=14)

                                                                                             gggggggggCGGCCGGGGGACGGGGCTGTCGCTCTATATTGGGTATGGCGGGGGCCGGCCTCTTCG‑C‑‑GGTTAACCCCGCTCCTACACTGGGCCGCGCCGCTTTGTTGGTGG‑GGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAg                                            <  2:3736025/148‑1 (MQ=25)

                                                                                                 gggggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACCGGGCCGCGCCGCTTTGTCGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTc                                        <  2:602624/148‑1 (MQ=255)

                                                                                                                    ggCTGTCGCTCGATAATCGGTATGGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGa                     <  1:5577637/150‑1 (MQ=255)

                                                                                                                    ggCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGATcc                                                                                                                 <  2:541397/59‑1 (MQ=255)

                                                                                                                    ggCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGATcc                                                                                                                 >  1:541397/1‑59 (MQ=255)

                                                                                                                            cTCGATAATCGGTATTGCCGGGGCCGGCCTCTTAG‑C‑‑GGGTAAAACCGCTACAAAACAGGGCCGCGCCGCGTTGTAGGAGC‑GGGTTTAACCGCGAAGAGGCCAGCACAGGCAATACAAAGCGTCAGGGCCGCTCAGCCCCCGAGTGgtcgt             >  2:4320521/1‑150 (MQ=38)

                                                                                                                             tCGATAATCGGTATTGCCGGGGCCGGCCACTTCG‑C‑‑GGGTAAAACCGCCCCGACACAGGGCCGCGCCGCTGTGTAGGCGC‑GGGTTTACCCGCGCAGCGGCCAGCACAGGCCATACAAAGCGTCAGGGCCGCGCAGCCACCGAGTGggcgtc            >  2:4272065/1‑150 (MQ=255)

                                                                                                                                    tttgTATTGCCTGTGCTGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAgg                                                                                       >  1:4338869/4‑69 (MQ=14)

                                                                                                                                    acGGTATGGCCGGGGCCGGCCCCTCCG‑C‑‑GGCAAACCCCGCACCTACACCGGGCCGCGCCGCTCTGTCGCCGC‑GGGTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCAcc     <  2:715402/149‑1 (MQ=255)

                                                                                                                                        tATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑gg                                                                                 >  2:5820673/1‑70 (MQ=255)

                                                                                                                                              cGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                         <  1:7322142/57‑1 (MQ=255)

CTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCG‑C‑‑GGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGC  >  NC_002947/1549301‑1549585

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