BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,550,196	A→C	5.1%	N185T (AAC→ACC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,550,196	0	A	C	5.1%	535.0 / 12.8	157	N185T (AAC→ACC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base A (70/79); new base C (5/3); total (75/82)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 172 total aligned reads displayed.

CAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGC  >  NC_002947/1550047‑1550345

cAAGCCATGCGCCGCCCCCACCGCCCTCGCCCACGTAGGTACCCTCTCTGCCAACTCCGACCACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGaaaa                                                                                                                                                       <  1:3409113/150‑1 (MQ=255)

cAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGaaaa                                                                                                                                                       >  2:3730789/1‑150 (MQ=255)

cAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGATCGGGCCTGGaaaa                                                                                                                                                       >  2:3730760/1‑150 (MQ=255)

 aaGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAAc                                                                                                                                                      >  1:2061347/1‑150 (MQ=255)

   gCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACgt                                                                                                                                                    <  1:877112/150‑2 (MQ=255)

    ccATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGaa                                                                                                                                                   <  2:3286465/150‑1 (MQ=255)

    ccATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGaa                                                                                                                                                   <  2:356346/150‑1 (MQ=255)

     cATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAAc                                                                                                                                                  <  1:1919548/150‑1 (MQ=255)

      aTGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACt                                                                                                                                                 >  1:3217828/1‑150 (MQ=255)

       tGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTg                                                                                                                                                <  2:763820/150‑1 (MQ=255)

          gcCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGGCTGGAAAACGAACTGTcg                                                                                                                                             >  2:3525968/1‑149 (MQ=255)

          gcCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTcc                                                                                                                                             >  2:2010047/1‑150 (MQ=255)

                tCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTa                                                                                                                                       >  2:5299811/1‑150 (MQ=255)

                tCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTa                                                                                                                                       <  2:3217828/150‑1 (MQ=255)

                  cAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGa                                                                                                                                     >  2:3978714/1‑150 (MQ=255)

                   aaGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGTTGAAATCATCGCCGAAGCCATCTAAAAAGTCGGTAAAGAAGGCGTGATCACAGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCcgtcgaatac                                                                                                                                    >  2:524131/1‑145 (MQ=255)

                    aGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAg                                                                                                                                   <  1:5299720/150‑1 (MQ=255)

                    aGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAg                                                                                                                                   <  1:5299550/150‑1 (MQ=255)

                       ccATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCa                                                                                                                                <  2:4112534/150‑1 (MQ=255)

                        cATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCAt                                                                                                                               <  2:4775834/150‑1 (MQ=255)

                         aTCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATg                                                                                                                              <  2:3342615/150‑1 (MQ=255)

                         aTCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATg                                                                                                                              <  2:3897575/150‑1 (MQ=255)

                         aTCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATg                                                                                                                              <  2:3342864/150‑1 (MQ=255)

                             cccAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGt                                                                                                                          >  1:3025222/1‑150 (MQ=255)

                               cAGGTAGGTACCCTCTCTGCCACCTCCGACCACTCCATCGGTGACAGCCTCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGACCTGTCCGTCGTAGAAGGCATGCAGTTc                                                                                                                        <  2:1551251/150‑1 (MQ=255)

                                     ggtaCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGaa                                                                                                                                    >  2:5001543/1‑132 (MQ=255)

                                     ggtaCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGaa                                                                                                                                    <  2:3829909/132‑1 (MQ=255)

                                     ggtaCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGaa                                                                                                                                    >  1:3829909/1‑132 (MQ=255)

                                     ggtaCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGaa                                                                                                                                    <  1:5001540/132‑1 (MQ=255)

                                     ggtaCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACcgc                                                                                                                  <  1:2355848/150‑1 (MQ=255)

                                     gggaCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACcgc                                                                                                                  <  1:3568002/147‑1 (MQ=255)

                                       taCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCgg                                                                                                                <  2:347278/150‑1 (MQ=255)

                                        aCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGc                                                                                                               >  2:2914554/1‑150 (MQ=255)

                                        aCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGc                                                                                                               >  2:4638977/1‑150 (MQ=255)

                                        aCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGc                                                                                                               >  2:4638657/1‑150 (MQ=255)

                                          cATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTa                                                                                                             >  1:1483914/1‑150 (MQ=255)

                                           aTCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGTAGTTCGACCGCGGCTAc                                                                                                            >  1:4194460/1‑150 (MQ=255)

                                            tctcTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAAc                                                                                                                                                      <  1:5304162/107‑1 (MQ=255)

                                            tctcTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTAcc                                                                                                           >  1:1315697/1‑150 (MQ=255)

                                                tGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACcgcg                                                                                                                 >  2:4100186/1‑140 (MQ=255)

                                                tGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTc                                                                                                       >  2:2454467/1‑150 (MQ=255)

                                                 gCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCg                                                                                                      <  2:4143668/150‑1 (MQ=255)

                                                   cAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGcc                                                                                                    <  2:352090/150‑1 (MQ=255)

                                                   cAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGcc                                                                                                    >  1:2512778/1‑150 (MQ=255)

                                                   cAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGcc                                                                                                    <  2:4186418/150‑1 (MQ=255)

                                                    aaCTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCg                                                                                                   <  1:362882/150‑1 (MQ=255)

                                                     aCTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCg                                                                                                   <  2:4490689/149‑1 (MQ=255)

                                                       tCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTAc                                                                                                <  1:4196342/150‑1 (MQ=255)

                                                          gACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTc                                                                                             <  1:2914554/150‑1 (MQ=255)

                                                          gACAACTCCAGCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTc                                                                                             <  2:226155/150‑1 (MQ=255)

                                                                tCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaac                                                                                       <  1:4610996/150‑1 (MQ=255)

                                                                tCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaac                                                                                       <  1:20025/150‑1 (MQ=255)

                                                                 ccATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaaca                                                                                      <  1:4681393/150‑1 (MQ=255)

                                                                  cATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGaaaa                                                                                                                                                       >  2:665542/1‑84 (MQ=255)

                                                                  cATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGaaaa                                                                                                                                                       <  1:665542/84‑1 (MQ=255)

                                                                   aTCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAg                                                                                    >  1:2662799/1‑150 (MQ=255)

                                                                        tGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGa                                                                               <  2:1483914/150‑1 (MQ=255)

                                                                            atcatcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACAcc                                                                           <  1:4638657/150‑1 (MQ=255)

                                                                            atcatcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACAcc                                                                           <  1:4638977/150‑1 (MQ=255)

                                                                               atcGCTGAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGacaagatg                                                                        >  2:3704644/1‑145 (MQ=255)

                                                                               atcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATg                                                                        <  1:4957426/150‑1 (MQ=255)

                                                                                tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                       >  1:2598555/1‑150 (MQ=255)

                                                                                tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                       >  2:4373908/1‑150 (MQ=255)

                                                                                tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                       >  1:1405263/1‑150 (MQ=255)

                                                                                tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                       >  1:4838727/1‑150 (MQ=255)

                                                                                  gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                     >  2:1640099/4‑150 (MQ=255)

                                                                                  gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                     >  1:375046/4‑150 (MQ=255)

                                                                                  gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                     <  1:708264/147‑1 (MQ=255)

                                                                                  gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                     >  2:1640083/4‑150 (MQ=255)

                                                                                   ccGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTg                                                                    >  2:4870210/1‑150 (MQ=255)

                                                                                    tgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGTTGc                                                                   <  2:1206031/149‑1 (MQ=255)

                                                                                    tgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGTTGc                                                                   >  2:15247/2‑150 (MQ=255)

                                                                                     gAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGcc                                                                  >  1:2190167/1‑150 (MQ=255)

                                                                                      aaGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCg                                                                 >  1:1799413/1‑150 (MQ=255)

                                                                                         ccATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGc                                                              >  1:4543638/1‑150 (MQ=255)

                                                                                         ccATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGc                                                              >  1:4543610/1‑150 (MQ=255)

                                                                                         ccATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGc                                                              <  2:2207475/150‑1 (MQ=255)

                                                                                                 aaaGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTGCCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGc                                                      >  1:3266636/1‑150 (MQ=255)

                                                                                                   aGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGccc                                                    <  2:2434435/150‑1 (MQ=255)

                                                                                                     tCGGTAAAGAAGGCGTGATTACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgc                                                  >  1:5029348/1‑150 (MQ=255)

                                                                                                     tCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgc                                                  <  2:2719932/150‑1 (MQ=255)

                                                                                                      cGGTTAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgct                                                 <  1:4461255/150‑1 (MQ=255)

                                                                                                        gTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgctgc                                               >  2:5070327/1‑150 (MQ=255)

                                                                                                           aagaagGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACt                                                                                                                                                 >  2:547825/1‑49 (MQ=255)

                                                                                                           aagaagGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgctgctgc                                            <  1:1042920/150‑1 (MQ=255)

                                                                                                                ggCGAGATCCCCGTTGAAGAAGGCTCGGGCCGGGAAAACGAACTGTCCGACGTAGAAGGCAAGCAGTTCGACCGCGGCAACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGCTTGCCGAGCTGGAAGGCCCGctgctgctgctg                                       >  2:412156/1‑150 (MQ=255)

                                                                                                                         aCCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAg                              >  1:2734522/1‑150 (MQ=255)

                                                                                                                          ccGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGa                             >  1:1020613/1‑150 (MQ=255)

                                                                                                                            gTTGAAGTAGGCTAGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGTCATGCAGTTCGACCGCGGCTACCTGTCGACGTACTTCGTCATCCAGCCGGACTCCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCAGCTGGTCGACAAGaaa                           <  2:841994/150‑1 (MQ=255)

                                                                                                                            gTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGaaa                           >  1:509875/1‑150 (MQ=255)

                                                                                                                                aagaagGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACACGAAAAtct                       >  2:1212196/1‑150 (MQ=255)

                                                                                                                                      gCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACa                 >  2:1517862/1‑150 (MQ=255)

                                                                                                                                       cTCGGGTCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACAt                <  2:670751/150‑1 (MQ=255)

                                                                                                                                        tCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATc               <  2:2169622/150‑1 (MQ=255)

                                                                                                                                             ccTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGa          <  2:3435918/150‑1 (MQ=255)

                                                                                                                                               tGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAgc        <  2:3025222/150‑1 (MQ=255)

                                                                                                                                                   aaaCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAgctgct    <  1:837723/150‑1 (MQ=255)

                                                                                                                                                     aCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAg                                                               <  2:1376258/89‑1 (MQ=255)

                                                                                                                                                     aCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCtgctgc  <  2:1700604/150‑1 (MQ=255)

                                                                                                                                                     aCGAACTGGCAGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGGTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCACCATCCGTGAGCtgctgc  <  1:3730760/150‑1 (MQ=255)

CAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGC  >  NC_002947/1550047‑1550345

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