BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,550,205	C→A	5.1%	S188Y (TCC→TAC) ‡	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,550,205	0	C	A	5.1%	523.4 / 12.7	156	S188Y (TCC→TAC) ‡	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base C (69/79); new base A (5/3); total (74/82)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 159 total aligned reads displayed.

CGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCT  >  NC_002947/1550056‑1550352

cgcCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTc                                                                                                                                                     <  2:1693006/150‑1 (MQ=255)

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

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

  ccGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCcg                                                                                                                                                   >  2:4217718/1‑150 (MQ=255)

       tCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGTACTGTCCGTCGTa                                                                                                                                              >  2:4569508/1‑150 (MQ=255)

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

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

       tCCAAGGCCATCGCCCAGGGAGGTACCATCTCTGCCAACTCCGACACCTCCATCGGTGAAATCATCGCCGCAGCCATGGACAAAGGCGGTCAAGAAGGCGGGATCACCGTGGACGAAGGCTCGGGCCTGGCACACGAACTGGCcgtcggc                                                                                                                                              >  2:3787619/1‑148 (MQ=255)

         ggtagCCCTCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTCGa                                                                                                                                            >  2:2573778/5‑150 (MQ=255)

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

           aGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAg                                                                                                                                          <  1:5299550/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)

                    cccAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAg                                                                                                                                          <  1:1908126/141‑1 (MQ=255)

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

                     ccAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGtt                                                                                                                                <  2:1527855/150‑1 (MQ=255)

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

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

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

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

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

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

                             gtaCCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACcgcg                                                                                                                        >  1:5091704/1‑150 (MQ=255)

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

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

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

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

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

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

                                  aTCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTAc                                                                                                                   <  1:3479130/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                             cTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTa                                                                                                        <  1:3360372/150‑1 (MQ=255)

                                              gccGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTAc                                                                                                       <  1:4733265/149‑1 (MQ=255)

                                                 gACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTc                                                                                                    >  1:839711/1‑150 (MQ=255)

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

                                                      cTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaa                                                                                               <  2:747816/150‑1 (MQ=255)

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

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

                                                              gTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCgg                                                                                       >  1:463132/1‑150 (MQ=255)

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

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

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

                                                                   atcatcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACAcc                                                                                  >  1:4257527/1‑150 (MQ=255)

                                                                    tgatcGCTGAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGacaaga                                                                                 <  2:2804242/148‑4 (MQ=255)

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

                                                                      atcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATg                                                                               >  2:3248376/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                                            gAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGcc                                                                         <  2:357938/150‑1 (MQ=255)

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

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

                                                                                 cATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCt                                                                    >  1:2072622/1‑150 (MQ=255)

                                                                                     gAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGaa                                                                <  2:3365230/150‑1 (MQ=255)

                                                                                     gAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGaa                                                                <  1:1549250/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                             aTCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGGACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGAc                                        <  2:411710/150‑1 (MQ=255)

                                                                                                               cACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACaa                                      >  1:4585825/1‑150 (MQ=255)

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

                                                                                                                 ccGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGa                                    <  2:2061347/150‑1 (MQ=255)

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

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

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

                                                                                                                      gaagaaGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAAtc                               <  1:877663/150‑1 (MQ=255)

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

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

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

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

                                                                                                                               tCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTATCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATc                      >  1:2278669/1‑150 (MQ=255)

                                                                                                                                cGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATcc                     <  2:1514774/150‑1 (MQ=255)

                                                                                                                                cGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATcc                     <  2:2244103/150‑1 (MQ=255)

                                                                                                                                  ggCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGt                   <  1:2593070/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                                aCTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGctgctgctgctg                                              <  1:4974126/109‑1 (MQ=255)

                                                                                                                                                  tGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTc   <  1:4569508/150‑1 (MQ=255)

                                                                                                                                                  tGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTc   >  1:817423/1‑150 (MQ=255)

                                                                                                                                                   gTCCGTCGTAGAAGGCATGAAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCt  <  1:524131/150‑1 (MQ=255)

CGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCT  >  NC_002947/1550056‑1550352

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