BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,357,039	T→C	5.8%	Q46Q (CAA→CAG) ‡	PP_2950 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,357,039	0	T	C	5.8%	978.6 / 30.3	309	Q46Q (CAA→CAG) ‡	PP_2950	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (159/130); new base C (15/3); total (176/133)
Fisher's exact test for biased strand distribution p-value = 2.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

Only 100 of 317 total aligned reads displayed.

CACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTC  >  NC_002947/3356890‑3357187

 aCTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCAtt                                                                                                                                                     >  1:9079699/1‑150 (MQ=255)

 aCTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCAtt                                                                                                                                                     >  1:1906464/1‑150 (MQ=255)

 aCTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCAtt                                                                                                                                                     >  2:8040426/1‑150 (MQ=255)

  cTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTg                                                                                                                                                    >  1:3071629/1‑150 (MQ=255)

   tCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTgg                                                                                                                                                   <  1:8963658/150‑1 (MQ=255)

     gTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGgt                                                                                                                                                 <  1:475754/150‑1 (MQ=255)

     gTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGgt                                                                                                                                                 >  1:4593412/1‑150 (MQ=255)

            gaaCGGACCGCGAGGCGGTGGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTGCTTCCACTTGCGAGTTAGTAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCAcc                                                                                                                                          <  2:6912048/149‑1 (MQ=255)

              aCGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGt                                                                                                                                        >  2:7816409/1‑150 (MQ=255)

               cGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGta                                                                                                                                       <  2:1155338/150‑1 (MQ=255)

               cGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTATGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGta                                                                                                                                       >  2:1565074/1‑150 (MQ=255)

                ggACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTAGCGAGTTAGGAAGTACCGAGCTTCTTCAAGGTTGTTTGCAAACTCGGCACTGTTACTCCATTGGGTGTCCACCGtat                                                                                                                                      >  2:3539441/1‑150 (MQ=255)

                    cgcgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTgg                                                                                                                                                   <  2:7263176/133‑1 (MQ=255)

                    cgcgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGc                                                                                                                                  >  1:2821487/1‑150 (MQ=255)

                      cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                                                                                <  1:6555529/150‑1 (MQ=255)

                      cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                                                                                <  2:4435372/150‑1 (MQ=255)

                      cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                                                                                <  2:2429247/150‑1 (MQ=255)

                        aTGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAgg                                                                                                                              <  2:4380418/150‑1 (MQ=255)

                                  cGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCtttt                                                                                                                    <  2:2835405/150‑1 (MQ=255)

                                    gTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCg                                                                                                                  >  2:3747073/1‑150 (MQ=255)

                                    gTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGAATAGGCCAGGGCAATCTTTTCg                                                                                                                  >  1:6403435/1‑150 (MQ=255)

                                     tGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCgg                                                                                                                 <  2:5523737/150‑1 (MQ=255)

                                      gAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCggg                                                                                                                >  1:1812276/1‑150 (MQ=255)

                                            ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                          >  1:7934174/1‑150 (MQ=255)

                                            ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                          <  2:237990/150‑1 (MQ=255)

                                            ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                          >  1:101337/1‑150 (MQ=255)

                                                 gCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                                                     >  1:7413657/1‑150 (MQ=255)

                                                     tttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCa                                                                                                 >  1:8149253/1‑150 (MQ=255)

                                                      ttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAg                                                                                                >  1:7247732/1‑150 (MQ=255)

                                                       tGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGc                                                                                               >  2:2270646/1‑150 (MQ=255)

                                                          tCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCAt                                                                                            >  1:5255898/1‑150 (MQ=255)

                                                          tCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCAt                                                                                            >  2:6299684/1‑150 (MQ=255)

                                                            aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                                                                          >  1:2875830/1‑150 (MQ=255)

                                                            aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                                                                          >  2:8329264/1‑150 (MQ=255)

                                                             aCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCt                                                                                         >  2:1725348/1‑150 (MQ=255)

                                                             aCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCt                                                                                         >  2:9142436/1‑150 (MQ=255)

                                                                gATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAg                                                                                      >  2:8904025/1‑150 (MQ=255)

                                                                    gTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAg                                                                                  >  1:8265565/1‑150 (MQ=255)

                                                                       cAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCg                                                                               >  1:1934175/1‑150 (MQ=255)

                                                                       cAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCg                                                                               >  1:884670/1‑150 (MQ=255)

                                                                          ttctttctTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGAc                                                                            >  1:4389425/1‑150 (MQ=255)

                                                                             tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                                                                         >  2:3821212/1‑150 (MQ=255)

                                                                             cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                                                                         >  2:1425777/5‑150 (MQ=255)

                                                                              ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGTCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                                        <  1:7583285/150‑1 (MQ=255)

                                                                              ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                                        <  2:3440458/150‑1 (MQ=255)

                                                                              ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                                        <  1:2664431/150‑1 (MQ=255)

                                                                              ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                                        <  1:4866419/150‑1 (MQ=255)

                                                                              ttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCttt                                                                        >  2:3987359/4‑150 (MQ=255)

                                                                               tcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtttt                                                                       <  2:7090428/150‑1 (MQ=255)

                                                                                 ttcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCTTTTcg                                                                     >  2:845976/1‑149 (MQ=255)

                                                                                 ttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCa                                                                     >  2:3215354/1‑150 (MQ=255)

                                                                                   ccACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAAt                                                                   >  1:3839940/1‑150 (MQ=255)

                                                                                    cACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATc                                                                  <  2:139986/150‑1 (MQ=255)

                                                                                    cACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATc                                                                  >  1:5771161/1‑150 (MQ=255)

                                                                                     aCTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCg                                                                 >  2:4880557/1‑150 (MQ=255)

                                                                                       ttGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGcc                                                               >  1:4100912/1‑150 (MQ=255)

                                                                                          cGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGag                                                            >  1:2437606/1‑149 (MQ=255)

                                                                                              ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                                        >  2:5006492/1‑150 (MQ=255)

                                                                                              ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                                        <  2:6832411/150‑1 (MQ=255)

                                                                                                aGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGa                                                      <  1:6214588/150‑1 (MQ=255)

                                                                                                 ggAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCATTTCAATCGCCGATTCGCGAg                                                     >  2:7755081/1‑150 (MQ=255)

                                                                                                  gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtttt                                                                       >  1:6636007/1‑131 (MQ=255)

                                                                                                  gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                                    >  2:6693899/1‑150 (MQ=255)

                                                                                                  gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                                    <  1:5171014/150‑1 (MQ=255)

                                                                                                  gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                                    >  2:6524965/1‑150 (MQ=255)

                                                                                                  gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACATTTTCAATCGCCGATTCGCGAGt                                                    >  1:4596436/1‑150 (MQ=255)

                                                                                                      tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                <  1:1712646/150‑1 (MQ=255)

                                                                                                      tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                <  2:2607102/150‑1 (MQ=255)

                                                                                                      tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                >  2:2047827/1‑150 (MQ=255)

                                                                                                      tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                >  2:5455230/1‑150 (MQ=255)

                                                                                                        ccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAggg                                                                                                                             >  2:4198699/1‑71 (MQ=255)

                                                                                                        ccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                                                     >  2:182504/1‑95 (MQ=255)

                                                                                                         cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                                             <  2:1270700/150‑1 (MQ=255)

                                                                                                         cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                                             <  1:8494567/150‑1 (MQ=255)

                                                                                                         cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                                             >  2:4207656/1‑150 (MQ=255)

                                                                                                         cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGGGAAagga                                             >  1:5366386/1‑150 (MQ=255)

                                                                                                          gAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAaggag                                            <  2:877221/150‑1 (MQ=255)

                                                                                                          gAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGGGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAaggag                                            >  2:4685849/1‑150 (MQ=255)

                                                                                                                 ttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTgg                                     >  1:719614/1‑150 (MQ=255)

                                                                                                                 ttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTgg                                     >  1:7530820/1‑150 (MQ=255)

                                                                                                                  tcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc                                    <  2:2821490/150‑1 (MQ=255)

                                                                                                                  tcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc                                    >  2:7572396/1‑150 (MQ=255)

                                                                                                                  tcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc                                    >  2:3310167/1‑150 (MQ=255)

                                                                                                                    aCGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCgg                                  >  2:8671766/1‑150 (MQ=255)

                                                                                                                          gttTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgt                            <  1:1234662/150‑1 (MQ=255)

                                                                                                                          gttTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgt                            <  1:5006492/150‑1 (MQ=255)

                                                                                                                            ttGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgtgt                          <  1:8764398/150‑1 (MQ=255)

                                                                                                                            ttGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgtgt                          <  2:9079699/150‑1 (MQ=255)

                                                                                                                              gCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgtgtgt                        <  2:7316176/150‑1 (MQ=255)

                                                                                                                                aaaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                      >  1:2926160/1‑150 (MQ=255)

                                                                                                                                aaaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                      >  1:8997268/1‑150 (MQ=255)

                                                                                                                                aaaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                      <  1:8012039/150‑1 (MQ=255)

                                                                                                                                 aaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGAACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCg                     >  2:5861182/1‑150 (MQ=255)

                                                                                                                                    tCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTg                  >  1:7679383/1‑150 (MQ=255)

                                                                                                                                      ggCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGa                <  1:5495868/150‑1 (MQ=255)

                                                                                                                                        cACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACa              <  2:4570344/150‑1 (MQ=255)

                                                                                                                                               aCGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATggtggt       <  2:2750473/150‑1 (MQ=255)

                                                                                                                                                   cATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTGCGTTGGACATGGTGGTGCCt   >  2:6121216/1‑150 (MQ=255)

                                                                                                                                                    aTTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTc  >  1:3397264/1‑150 (MQ=255)

                                                                                                                                                    aTTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTc  >  2:2898825/1‑150 (MQ=255)

CACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTC  >  NC_002947/3356890‑3357187

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