BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,953	T→C	5.1%	L104P (CTG→CCG) ‡	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,953	0	T	C	5.1%	641.2 / 13.6	197	L104P (CTG→CCG) ‡	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base T (97/90); new base C (7/3); total (104/93)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 202 total aligned reads displayed.

TTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACT  >  NC_002947/1549804‑1550102

 tCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAACACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTg                                                                                                                                                      <  2:842396/150‑1 (MQ=255)

 tCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTg                                                                                                                                                      >  1:4458275/1‑150 (MQ=255)

  ccGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                                     <  1:2738230/150‑1 (MQ=255)

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

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

       gCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACAACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGcc                                                                                                                                                >  1:2219153/1‑150 (MQ=255)

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

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

         cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                                                                                                                              >  1:4284876/1‑150 (MQ=255)

          aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                                                                                                                             >  1:1236379/1‑150 (MQ=255)

              aaaTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTg                                                                                                                                         >  2:4402455/1‑150 (MQ=255)

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

                         aaaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAAc                                                                                                                              >  2:5299394/1‑150 (MQ=255)

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

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

                             aCGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTCCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGa                                                                                                                          >  1:3587997/1‑150 (MQ=255)

                              cGCTTTCGAAACCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAt                                                                                                                         <  2:645752/150‑1 (MQ=255)

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

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

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

                               aagttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATg                                                                                                                        >  2:465880/4‑150 (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:1303961/1‑150 (MQ=255)

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

                                           aTGGGCGCCCAGCTGGTCAAGGACGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAgcgc                                                                                                            <  2:4558128/150‑1 (MQ=255)

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

                                            tGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGcgcg                                                                                                           <  2:3971086/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                             aaGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCAct                                                                                          <  2:5036373/150‑1 (MQ=255)

                                                              aGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCActg                                                                                         >  1:4341707/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                          ccAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTg                                                                             >  1:936506/1‑150 (MQ=255)

                                                                           cAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGc                                                                            >  2:4505378/1‑150 (MQ=255)

                                                                            aaGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGcc                                                                           >  1:404316/1‑150 (MQ=255)

                                                                             aGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCg                                                                          <  2:4685471/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                     gCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCAt                                                  <  1:143522/150‑1 (MQ=255)

                                                                                                       accaccaccGCTTCCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATgc                                                <  2:4018257/150‑1 (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)

                                                                                                            caccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCTACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGa                                           <  1:1707257/150‑1 (MQ=255)

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

                                                                                                                       ttCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTg                                                                             <  1:5128357/105‑1 (MQ=255)

                                                                                                                       ttCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTg                                                                             >  2:5128357/1‑105 (MQ=255)

                                                                                                                       ttCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCa                                >  1:763820/1‑150 (MQ=255)

                                                                                                                       acggctttCAGGCCTTCGTTGACGATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCa                                >  1:1878569/8‑150 (MQ=255)

                                                                                                                            tttCAGGCCTTCGTTGACGATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGcc                           >  2:3479130/3‑150 (MQ=255)

                                                                                                                             cTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGCTGGACGTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGCCCTGGAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCCTCGccc                          >  2:3409113/1‑150 (MQ=255)

                                                                                                                              tCAGGCCTTCGTTGACGATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGc                                                     >  2:1215909/1‑122 (MQ=255)

                                                                                                                              tCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCa                         <  2:4015790/150‑1 (MQ=255)

                                                                                                                               cAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCag                        >  1:4709174/1‑150 (MQ=255)

                                                                                                                                 ggCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggt                      <  1:2883467/150‑1 (MQ=255)

                                                                                                                                 ggCCATCGTCAACGAAGGCCTG‑‑AGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTg                                                                             >  2:4944048/1‑93 (MQ=255)

                                                                                                                                  gCCTTCGTTGACGATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtc                                                                                 >  2:4197762/1‑90 (MQ=255)

                                                                                                                                  gCCTTCGTTGACGATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGcc                                             <  1:1635480/126‑1 (MQ=255)

                                                                                                                                  gCCTTCGTTGACGATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGcc                                             >  2:1635480/1‑126 (MQ=255)

                                                                                                                                   ccATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggtag                    <  2:1352239/150‑1 (MQ=255)

                                                                                                                                    cATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggtagg                   >  1:894757/1‑150 (MQ=255)

                                                                                                                                      tCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggtaggtc                 >  2:362882/1‑149 (MQ=255)

                                                                                                                                         tCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCa              >  2:2649368/1‑150 (MQ=255)

                                                                                                                                                aGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGc       <  2:2170819/150‑1 (MQ=255)

                                                                                                                                                aGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGc       <  2:274733/150‑1 (MQ=255)

                                                                                                                                                   ccTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCaa    <  2:3161241/150‑1 (MQ=255)

                                                                                                                                                    cTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAAc   >  2:2355848/1‑150 (MQ=255)

TTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACT  >  NC_002947/1549804‑1550102

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