BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,550,203	G→C	6.2%	L187L (CTG→CTC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,550,203	0	G	C	6.2%	517.7 / 15.0	157	L187L (CTG→CTC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base G (71/76); new base C (5/5); total (76/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

TGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCT  >  NC_002947/1550054‑1550352

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

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

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

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

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

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

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

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

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

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

                  aTCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCagg                                                                                                                                     >  2:3529528/1‑148 (MQ=255)

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

                   tCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGc                                                                                                                                    >  2:472586/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                              aCTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGt                                                                                                         <  2:2919858/150‑1 (MQ=255)

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

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

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

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

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

                                                      aaCTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTc                                                                                                 <  1:4924865/150‑1 (MQ=255)

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

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

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

                                                          ccATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaaca                                                                                             >  2:586285/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                             cGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGc                                                                          >  1:754393/1‑150 (MQ=255)

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

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

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

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

                                                                                       gAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGaa                                                                >  2:372604/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                        acgAAAACCAACTCTCCGTCCTAGAACGCAACCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCCGACACCATGCTAGCCGAGCTGGAAGGCCCGCTGCTGCTGCACGTCGACAACAAAATCTCCAACATCCGTGAgc               <  2:3490645/148‑1 (MQ=255)

                                                                                                                                            aaaCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAgctgct           >  2:346314/1‑150 (MQ=255)

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

                                                                                                                                             aacgaacgGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGctgctg          <  1:3787619/150‑1 (MQ=255)

                                                                                                                                             aaCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGctgctg          >  1:5181232/1‑150 (MQ=255)

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

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

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

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

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

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

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

TGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCTGCCAACTCCGACAACTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCT  >  NC_002947/1550054‑1550352

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