BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,927	G→C	9.9%	L95L (CTG→CTC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,927	0	G	C	9.9%	606.5 / 45.2	192	L95L (CTG→CTC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base G (87/86); new base C (13/6); total (100/92)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 196 total aligned reads displayed.

GCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCA  >  NC_002947/1549778‑1550072

gcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTg                                                                                                                                                   >  2:1002034/1‑150 (MQ=255)

 cgcCGACCATCACCCAGGACGGCGTTTCCGTCGCCAAAGAACTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTgg                                                                                                                                                  <  2:516871/150‑1 (MQ=255)

  gcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGc                                                                                                                                                 >  1:2908653/1‑150 (MQ=255)

     gACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCa                                                                                                                                              <  1:1867578/150‑1 (MQ=255)

     acgcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCa                                                                                                                                              >  1:3685903/4‑150 (MQ=255)

      aCCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAg                                                                                                                                             <  1:5339013/150‑1 (MQ=255)

      aCCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAg                                                                                                                                             <  2:4189310/150‑1 (MQ=255)

       gcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAgg                                                                                                                                            <  1:4216420/149‑1 (MQ=255)

       gcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAgg                                                                                                                                            >  1:2093184/2‑150 (MQ=255)

       ccATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAgg                                                                                                                                            <  2:4027879/150‑1 (MQ=255)

       ccATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAgg                                                                                                                                            <  1:184726/150‑1 (MQ=255)

            aCCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATc                                                                                                                                       >  1:1434979/1‑150 (MQ=255)

            aCCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATc                                                                                                                                       <  1:3305332/150‑1 (MQ=255)

            aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                                                                                                                                       <  1:1967786/150‑1 (MQ=255)

            aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                                                                                                                                       >  2:1390641/1‑150 (MQ=255)

            aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                                                                                                                                       <  2:2783581/150‑1 (MQ=255)

                  gACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGttcag                                                                                                                                 >  1:1087694/1‑146 (MQ=255)

                     ggCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGaa                                                                                                                              >  1:1465819/1‑150 (MQ=255)

                      gCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAg                                                                                                                             <  2:2964907/150‑1 (MQ=255)

                          ttCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCt                                                                                                                         <  2:5018223/150‑1 (MQ=255)

                              gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                     <  2:3339704/150‑1 (MQ=255)

                              gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                     <  1:4020978/150‑1 (MQ=255)

                              gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                     <  1:2074508/150‑1 (MQ=255)

                                  cctaaGAATTCGAGCTGAAAGTCGCTTTCGAATACCTGGGCGCCCAGCTGGTCTAGGATGTTGCTTCCAAGGCCAACGACGGTGCCGGTGACGGCAGCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCg                                                                                                                 <  2:2002704/147‑1 (MQ=255)

                                   gcaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                                                                                                <  2:4907492/148‑1 (MQ=255)

                                   cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCGCCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                                                                                                >  2:4522324/1‑150 (MQ=255)

                                    aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                                                                                               <  1:4070381/150‑1 (MQ=255)

                                             gAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGc                                                                                                      >  2:4107260/1‑150 (MQ=255)

                                                cTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATg                                                                                                   >  2:5045400/1‑150 (MQ=255)

                                                 tGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGa                                                                                                  >  2:4323395/1‑150 (MQ=255)

                                                      gACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCg                                                                                             >  2:3768277/1‑150 (MQ=255)

                                                       acaagttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAAcccaa                                                                                            >  1:54209/6‑148 (MQ=255)

                                                       aCGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGa                                                                                            >  2:782035/1‑150 (MQ=255)

                                                        caagttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAAt                                                                                           <  2:4702801/146‑1 (MQ=255)

                                                         gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                          >  2:1853486/1‑150 (MQ=255)

                                                         gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                          <  1:782035/150‑1 (MQ=255)

                                                         gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                          <  1:55358/150‑1 (MQ=255)

                                                         gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                          <  2:2341306/150‑1 (MQ=255)

                                                         gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                          <  1:2744422/150‑1 (MQ=255)

                                                         aagttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATg                                                                                          >  2:465880/4‑150 (MQ=255)

                                                          agttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATgg                                                                                         <  1:4775466/148‑1 (MQ=255)

                                                           tttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGa                                                                                        >  2:2863731/1‑150 (MQ=255)

                                                            ttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGaa                                                                                                 >  1:1174108/1‑140 (MQ=255)

                                                            ttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATGGAc                                                                                       >  1:2847167/1‑150 (MQ=255)

                                                             tCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAcc                                                                                      >  2:2087280/1‑150 (MQ=255)

                                                               gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTg                                                                                    >  1:1246359/1‑150 (MQ=255)

                                                               gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTg                                                                                    >  1:4321437/1‑150 (MQ=255)

                                                               gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAGGGACCTg                                                                                    >  1:932950/1‑150 (MQ=255)

                                                                 aaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGaa                                                                                  >  2:4595444/1‑150 (MQ=255)

                                                                     aTGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAgcgc                                                                              >  1:3873803/1‑150 (MQ=255)

                                                                       gggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGc                                                                                                          >  1:4456802/1‑120 (MQ=255)

                                                                       gggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGc                                                                                                          <  2:4456802/120‑1 (MQ=255)

                                                                           gcCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATc                                                                        <  1:3953869/150‑1 (MQ=255)

                                                                           gcCCAGCTGGCCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATc                                                                        <  2:1678546/150‑1 (MQ=255)

                                                                            cccAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCg                                                                       >  2:4883305/1‑150 (MQ=255)

                                                                                 cTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAg                                                                  <  1:2214659/150‑1 (MQ=255)

                                                                                  tGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAgg                                                                 <  2:4078260/150‑1 (MQ=255)

                                                                                  tGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAgg                                                                 <  1:2863731/150‑1 (MQ=255)

                                                                                   ggTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                <  2:549646/150‑1 (MQ=255)

                                                                                    gTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGcc                                                               >  1:266473/1‑150 (MQ=255)

                                                                                     tCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCa                                                                                                     <  2:4493767/111‑1 (MQ=255)

                                                                                     tCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCa                                                                                                     >  1:4493755/1‑111 (MQ=255)

                                                                                      cAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCt                                                                                                                         <  2:214223/90‑1 (MQ=255)

                                                                                         ggAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCActgc                                                          >  2:3592190/1‑150 (MQ=255)

                                                                                         ggAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCActgc                                                          >  1:3286465/1‑150 (MQ=255)

                                                                                                 cTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGt                                                                                                                                     <  1:1729751/67‑1 (MQ=255)

                                                                                                 cTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGt                                                                                                                                     >  2:1729751/1‑67 (MQ=255)

                                                                                                  ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                 >  1:1693006/1‑150 (MQ=255)

                                                                                                  ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                 >  1:3213694/1‑150 (MQ=255)

                                                                                                   tCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGtt                                                >  2:1955200/1‑150 (MQ=255)

                                                                                                   tCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGtt                                                <  1:226352/150‑1 (MQ=255)

                                                                                                   tCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGtt                                                >  2:2927720/1‑150 (MQ=255)

                                                                                                    ccAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCAc                                                             >  1:4857979/1‑136 (MQ=255)

                                                                                                        ggCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGACCCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGa                                           <  1:3621561/150‑1 (MQ=255)

                                                                                                        ggCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGa                                           <  1:2199462/150‑1 (MQ=255)

                                                                                                          ccAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAAc                                         <  2:2720040/150‑1 (MQ=255)

                                                                                                          ccAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAAc                                         <  1:1232408/150‑1 (MQ=255)

                                                                                                              cgacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaa                                     <  1:5373765/150‑1 (MQ=255)

                                                                                                               gacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaag                                    <  2:1264296/150‑1 (MQ=255)

                                                                                                                acgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAATGCCACTGCTGCCGTCGTTGCCGAACTgaaga                                   >  1:2459749/1‑150 (MQ=255)

                                                                                                                 cgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaagaa                                  <  1:3443913/150‑1 (MQ=255)

                                                                                                                 cgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaagaa                                  >  1:4112534/1‑150 (MQ=255)

                                                                                                                 cgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaagaa                                  <  2:1522263/150‑1 (MQ=255)

                                                                                                                    tGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCt                               >  1:3915529/1‑150 (MQ=255)

                                                                                                                      ccGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                >  2:4971735/1‑115 (MQ=255)

                                                                                                                      ccGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                <  1:4971936/115‑1 (MQ=255)

                                                                                                                          tGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                         <  1:1492818/150‑1 (MQ=255)

                                                                                                                          tGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                         <  1:860358/150‑1 (MQ=255)

                                                                                                                           gACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAg                        <  2:3990380/150‑1 (MQ=255)

                                                                                                                              ggCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCa                     >  1:4101378/1‑150 (MQ=255)

                                                                                                                              ggCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCa                     >  1:347278/1‑150 (MQ=255)

                                                                                                                              ggCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCa                     >  1:357938/1‑150 (MQ=255)

                                                                                                                                 accaccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATgc                  >  1:3342864/1‑150 (MQ=255)

                                                                                                                                   caccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATgcgc                >  1:4567916/1‑150 (MQ=255)

                                                                                                                                    accaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGcc               <  1:2958613/150‑1 (MQ=255)

                                                                                                                                     ccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCTGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCg              >  2:1789742/1‑150 (MQ=255)

                                                                                                                                     ccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCg              >  1:1527855/1‑150 (MQ=255)

                                                                                                                                            tACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCaa       <  2:3383552/150‑1 (MQ=255)

                                                                                                                                            tACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCaa       <  2:932950/150‑1 (MQ=255)

                                                                                                                                              ccGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAgg     <  1:3614934/150‑1 (MQ=255)

GCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCA  >  NC_002947/1549778‑1550072

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