BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,284,252	C→T	8.0%	E27K (GAG→AAG)	PP_2883 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,284,252	0	C	T	8.0%	480.6 / 22.9	163	E27K (GAG→AAG)	PP_2883	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (71/79); new base T (8/5); total (79/84)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 165 total aligned reads displayed.

GCTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTG  >  NC_002947/3284105‑3284399

 cTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                                   <  2:1206297/150‑1 (MQ=255)

 cTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                                   >  1:868007/1‑150 (MQ=255)

  tGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAgg                                                                                                                                                  <  1:985700/150‑1 (MQ=255)

  tGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAgg                                                                                                                                                  <  1:985693/150‑1 (MQ=255)

    ccAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtg                                                                                                                                                >  2:3630385/1‑150 (MQ=255)

     cAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg                                                                                                                                               >  2:676262/1‑150 (MQ=255)

     cAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg                                                                                                                                               <  2:5101344/150‑1 (MQ=255)

      aGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg‑t                                                                                                                                             <  2:530764/150‑1 (MQ=255)

          ttCGAGGCCGGCCTGTTCTCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCATGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCAc                                                                                                                                         <  1:4434049/150‑1 (MQ=255)

          ttCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCAc                                                                                                                                         >  2:1221944/1‑150 (MQ=255)

          ttCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCA‑CCACCTTGAGGTGGCTGCAc                                                                                                                                         <  1:4193401/150‑1 (MQ=255)

          ttCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCA‑CCACCTTGAGGTGGCTGCAc                                                                                                                                         <  1:713427/150‑1 (MQ=255)

             gAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTc                                                                                                                                                      <  2:3574959/135‑1 (MQ=255)

                gccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAg                                                                                                                                   >  2:5169356/1‑150 (MQ=255)

                 ccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAgg                                                                                                                                  >  2:1096694/1‑150 (MQ=255)

                  cggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAggg                                                                                                                                 >  2:786502/1‑150 (MQ=255)

                  cggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAggg                                                                                                                                 <  2:292753/150‑1 (MQ=255)

                    gccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGcc                                                                                                                               >  2:4022476/1‑150 (MQ=255)

                         ttCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACAtt                                                                                                                          >  1:4842712/1‑150 (MQ=255)

                           cGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTgc                                                                                                                        >  2:2967397/1‑150 (MQ=255)

                           cGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTgc                                                                                                                        >  1:2011867/1‑150 (MQ=255)

                           cGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTgc                                                                                                                        >  2:2967633/1‑150 (MQ=255)

                            gCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTgcg                                                                                                                       <  2:4842712/150‑1 (MQ=255)

                            gCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTgcg                                                                                                                       >  2:5080330/1‑150 (MQ=255)

                              cGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCa                                                                                                                     <  1:4293111/150‑1 (MQ=255)

                               gTCGGGGGCTTGGCACCGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAg                                                                                                                    >  2:2122269/1‑150 (MQ=255)

                                tCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAgg                                                                                                                   <  1:5080330/150‑1 (MQ=255)

                                  gggggCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCa                                                                                                                 >  1:3902154/1‑150 (MQ=255)

                                       cTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGt                                                                                                            <  1:5153680/150‑1 (MQ=255)

                                        ttGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTa                                                                                                           <  1:1300861/150‑1 (MQ=255)

                                          ggCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTgtggcc                                                                                                         <  1:3110699/150‑6 (MQ=255)

                                           gCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGAc                                                                                                        >  1:4646646/1‑150 (MQ=255)

                                            cAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACa                                                                                                       >  1:338428/1‑150 (MQ=255)

                                             aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACAt                                                                                                      >  2:4463370/1‑150 (MQ=255)

                                              acgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGTGGCGTGCA‑CCACCTTGAGGTGGCTGCAc                                                                                                                                         >  1:5013703/1‑114 (MQ=38)

                                              acgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGTGGCGTGCA‑CCACCTTGAGGTGGCTGCAc                                                                                                                                         <  2:5013703/114‑1 (MQ=38)

                                               cgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATcc                                                                                                    >  2:1614434/1‑150 (MQ=255)

                                                gaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATccc                                                                                                   >  2:138535/1‑150 (MQ=255)

                                                gaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCGGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATccc                                                                                                   >  2:1536324/1‑150 (MQ=255)

                                                  acgCAAGGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCgg                                                                                                 <  2:4641838/150‑1 (MQ=255)

                                                    gCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGc                                                                                               <  2:3932988/150‑1 (MQ=255)

                                                     cAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGGAGACATCCCGGGCg                                                                                              >  2:5233980/1‑150 (MQ=255)

                                                      aaTGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGt                                                                                             >  2:4634923/1‑150 (MQ=255)

                                                      aaTGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGt                                                                                             >  2:730956/1‑150 (MQ=255)

                                                          cgcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCgg                                                                                         <  1:1536324/150‑1 (MQ=255)

                                                           gcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGt                                                                                        <  2:707965/150‑1 (MQ=255)

                                                              aGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTgg                                                                                                             <  1:5090071/126‑1 (MQ=255)

                                                              aGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTgg                                                                                                             >  2:5090134/1‑126 (MQ=255)

                                                               ggCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAAc                                                                                    <  2:902979/150‑1 (MQ=255)

                                                                   ttCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCt                                                                                <  2:2846275/150‑1 (MQ=255)

                                                                    tCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtg                                                                               <  1:1096694/150‑1 (MQ=255)

                                                                    tCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtg                                                                               <  1:1096895/150‑1 (MQ=255)

                                                                       gTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTgg                                                                            <  2:3490197/150‑1 (MQ=255)

                                                                       gTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTgg                                                                            >  1:5107216/1‑150 (MQ=255)

                                                                        tATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAACGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCGGTggg                                                                           >  1:4445070/1‑150 (MQ=255)

                                                                          tCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCTAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTc                                                                         <  2:1735594/150‑1 (MQ=255)

                                                                            cAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCaa                                                                       <  1:449206/150‑1 (MQ=255)

                                                                            cAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCA‑CCACCTTGAggtgg                                                                                                                                               <  2:5055885/78‑1 (MQ=38)

                                                                            cAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCA‑CCACCTTGAggtgg                                                                                                                                               >  1:5055885/1‑78 (MQ=38)

                                                                                 ccAACGTCAGTTGCCATCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCGAGATCa                                                                  <  2:175625/150‑1 (MQ=255)

                                                                                 ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCGAGATCa                                                                  <  2:175863/150‑1 (MQ=255)

                                                                                 ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCa                                                                  <  1:786502/150‑1 (MQ=255)

                                                                                   aaCGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTa                                                                                                           <  2:3850030/107‑1 (MQ=255)

                                                                                   aaCGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTa                                                                                                           >  1:3850030/1‑107 (MQ=255)

                                                                                   aaCGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGGGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGGTCGGTGCACGTCTGGGGGGCAAGATCAgc                                                                >  1:275114/1‑149 (MQ=255)

                                                                                      gTCAGTTGCCAGCCCTGGATGCTGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAAt                                                             >  1:1737174/1‑150 (MQ=255)

                                                                                         aGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCg                                                          >  2:32129/1‑150 (MQ=255)

                                                                                             gCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCgg                                                      <  1:58228/150‑1 (MQ=255)

                                                                                              ccAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCACGATCAGAAATTCGCCGGc                                                     <  2:867834/150‑1 (MQ=255)

                                                                                              ccAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGc                                                     <  2:868007/150‑1 (MQ=255)

                                                                                               cAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCg                                                    >  2:4668614/1‑150 (MQ=255)

                                                                                                aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgt                                                   >  2:4941575/1‑150 (MQ=255)

                                                                                                  cccTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgtgt                                                 >  1:4137063/1‑150 (MQ=255)

                                                                                                     tGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACAATGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTgg                                              >  1:677037/1‑150 (MQ=255)

                                                                                                            cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCt                                          >  1:3613310/1‑147 (MQ=255)

                                                                                                             gACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTgg                                                                            <  1:4499372/112‑1 (MQ=255)

                                                                                                             gACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTgg                                                                            >  2:4499372/1‑112 (MQ=255)

                                                                                                             gACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGCTCAGAAATTCGCCGGCGTGTTGGCGCtgatg                                      <  1:1614434/150‑1 (MQ=255)

                                                                                                                    aTCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCgg                               >  2:3069565/1‑150 (MQ=255)

                                                                                                                      cAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGcc                             <  1:81776/150‑1 (MQ=255)

                                                                                                                       aGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCt                            <  1:3630385/150‑1 (MQ=255)

                                                                                                                          ccgccATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGGTGGCGCTGATGATCTCGGCCTCAt                         >  1:3337540/1‑150 (MQ=255)

                                                                                                                           cgccATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATc                        <  1:2915185/150‑1 (MQ=255)

                                                                                                                              cATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCgg                     >  2:2420859/1‑150 (MQ=255)

                                                                                                                              cATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCgg                     >  1:1156109/1‑150 (MQ=255)

                                                                                                                               aTCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCggg                    <  2:4137063/150‑1 (MQ=255)

                                                                                                                                  ttggcgTGCA‑CCACCTTGAGGTGGCTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCg                 >  1:3435674/3‑150 (MQ=255)

                                                                                                                                  gcggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCg                 >  1:2970775/1‑150 (MQ=255)

                                                                                                                                   cggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGa                <  2:3894115/150‑1 (MQ=255)

                                                                                                                                    ggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt               >  1:3683559/1‑150 (MQ=255)

                                                                                                                                      cgTGCA‑CCACCTTGAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGtt            >  1:4030569/1‑150 (MQ=255)

                                                                                                                                       gTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGtt            <  1:5169356/150‑1 (MQ=255)

                                                                                                                                        gtgcaccACCTTGAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTc           >  2:3864666/6‑150 (MQ=255)

                                                                                                                                          cAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCtt         <  2:2189142/150‑1 (MQ=255)

                                                                                                                                          cAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCtt         <  1:1952870/150‑1 (MQ=255)

                                                                                                                                            gCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGa       >  2:2972545/1‑150 (MQ=255)

                                                                                                                                             ccACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCTGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGaa      <  1:1838431/150‑1 (MQ=255)

                                                                                                                                              cACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGaaa     <  2:3426485/150‑1 (MQ=255)

                                                                                                                                              cACCTCAAGGTGG‑TGCACGCCAAGGGCCACAGTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGaaa     <  2:4874142/150‑1 (MQ=255)

                                                                                                                                               aCCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAAc    >  2:2174376/1‑150 (MQ=255)

GCTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTG  >  NC_002947/3284105‑3284399

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