BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,882,422	T→G	5.1%	A355A (GCT→GCG)	PP_4289 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,882,422	0	T	G	5.1%	540.3 / 12.0	158	A355A (GCT→GCG)	PP_4289	membrane protein
Reads supporting (aligned to +/- strand): ref base T (70/80); new base G (6/2); total (76/82)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 203 total aligned reads displayed.

CGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATCCAGGCCCA  >  NC_002947/4882275‑4882571

 gCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCcaggccgctgcccaggccgctgcc                                                                                                                                                    <  2:3628293/150‑1 (MQ=255)

  cAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCaggccgctgcccaggccgctgccc                                                                                                                                                   >  2:2779641/1‑150 (MQ=255)

   aGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAggccgctgcccaggccgctgccca                                                                                                                                                  <  2:1795160/150‑1 (MQ=255)

      gAAGATCGAGTACCAGCCGCTGCCGGCACCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGACGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCcgctgcccaggccgctgcccaggc                                                                                                                                               >  2:5440537/1‑150 (MQ=255)

      gAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCcgctgcccaggccgctgcccaggc                                                                                                                                               <  1:570933/150‑1 (MQ=255)

      gAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCcgctgcccaggccgctgcccaggc                                                                                                                                               >  1:3708160/1‑150 (MQ=255)

              aGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTg                                                                                                                                       <  1:1517711/150‑1 (MQ=255)

              aGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTg                                                                                                                                       >  1:3218709/1‑150 (MQ=255)

                  ccAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGa                                                                                                                                   >  1:1310700/1‑150 (MQ=255)

                   cAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGaa                                                                                                                                  >  1:1265085/1‑150 (MQ=255)

                    aGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAAt                                                                                                                                 <  2:239408/150‑1 (MQ=255)

                     gCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGACtt                                                                                                                                >  1:4969079/1‑150 (MQ=255)

                     gCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAAtt                                                                                                                                <  2:1310472/150‑1 (MQ=255)

                      ccGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTc                                                                                                                               >  2:2967005/1‑150 (MQ=255)

                      ccGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTc                                                                                                                               <  2:3737378/150‑1 (MQ=255)

                        gCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGa                                                                                                                             <  1:2967005/150‑1 (MQ=255)

                          tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                           <  1:2157997/150‑1 (MQ=255)

                          tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                           <  2:5262056/150‑1 (MQ=255)

                          tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                           <  1:232953/150‑1 (MQ=255)

                          tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGACGAATTCGACa                                                                                                                           <  2:2647509/150‑1 (MQ=255)

                           gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                          <  1:263027/150‑1 (MQ=255)

                           gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                          <  1:263252/150‑1 (MQ=255)

                           gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                          <  2:2318507/150‑1 (MQ=255)

                              ggCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGAt                                                                                                                       >  1:3828215/1‑150 (MQ=255)

                                 aaCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCca                                                                                                                    <  2:1421109/150‑1 (MQ=255)

                                 aaCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCca                                                                                                                    <  2:4110233/150‑1 (MQ=255)

                                  aCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCcac                                                                                                                   <  2:298389/150‑1 (MQ=255)

                                    tGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACaa                                                                                                                 >  1:17079/1‑150 (MQ=255)

                                    tGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACaa                                                                                                                 >  1:79246/1‑150 (MQ=255)

                                    tGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACaa                                                                                                                 >  1:79156/1‑150 (MQ=255)

                                         tAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCa                                                                                                            >  1:189377/1‑150 (MQ=255)

                                            gcggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATcc                                                                                                         >  1:4634319/1‑150 (MQ=255)

                                              ggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAg                                                                                                       >  2:2110102/1‑150 (MQ=255)

                                                   gACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACg                                                                                                  >  1:1480570/1‑150 (MQ=255)

                                                     cGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCt                                                                                                <  2:879663/150‑1 (MQ=255)

                                                     cGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCt                                                                                                <  1:965769/150‑1 (MQ=255)

                                                             gAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtg                                                                                        >  2:3312784/1‑150 (MQ=255)

                                                              agcagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtgt                                                                                       <  2:1123522/150‑1 (MQ=255)

                                                               gcagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtgtg                                                                                      <  1:4316512/150‑1 (MQ=255)

                                                                    gTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACt                                                                                 <  2:821028/150‑1 (MQ=255)

                                                                     tGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTc                                                                                <  2:3708070/150‑1 (MQ=255)

                                                                       ccGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGt                                                                              <  1:3034914/150‑1 (MQ=255)

                                                                        cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCTAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                             <  1:835546/150‑1 (MQ=255)

                                                                        cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                             <  2:3897017/150‑1 (MQ=255)

                                                                        cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                             <  2:4972793/150‑1 (MQ=255)

                                                                        cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGGGTGCCACTCGTc                                                                             >  2:367596/1‑150 (MQ=255)

                                                                         gAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTcc                                                                            >  2:4366945/1‑150 (MQ=255)

                                                                          aGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCa                                                                           >  1:5211106/1‑150 (MQ=255)

                                                                               gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                      >  2:4242487/1‑150 (MQ=255)

                                                                               gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                      >  2:2617199/1‑150 (MQ=255)

                                                                               gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                      >  2:5150091/1‑150 (MQ=255)

                                                                               gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                      >  2:1117188/1‑150 (MQ=255)

                                                                                ccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGa                                                                     <  1:5296674/150‑1 (MQ=255)

                                                                                  gccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCATTGTGTGCCACTCGTCCAAGCCGAcc                                                                   >  1:5058930/1‑150 (MQ=255)

                                                                                        gCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCa                                                             >  1:4558929/1‑150 (MQ=255)

                                                                                            ccGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCTGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGACGACCAGCCCACTGt                                                         >  2:1182337/1‑150 (MQ=255)

                                                                                            ccGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGACTTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGt                                                         >  2:5072244/1‑150 (MQ=255)

                                                                                             cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgt                                                                                         <  1:1537917/117‑1 (MQ=255)

                                                                                             cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgt                                                                                         >  2:1537917/1‑117 (MQ=255)

                                                                                             cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                                        <  1:4866514/150‑1 (MQ=255)

                                                                                             cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                                        <  1:2511482/150‑1 (MQ=255)

                                                                                             cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCAAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                                        >  1:1663886/1‑150 (MQ=255)

                                                                                              gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                                       <  2:4378418/150‑1 (MQ=255)

                                                                                              gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                                       <  2:258439/150‑1 (MQ=255)

                                                                                               aaGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCa                                                      <  2:1671629/150‑1 (MQ=255)

                                                                                               aaGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCa                                                      <  1:2617199/150‑1 (MQ=255)

                                                                                                  cACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCa                                                   >  1:3132807/1‑150 (MQ=255)

                                                                                                     ccgccgCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCg                                                >  1:4233710/1‑150 (MQ=255)

                                                                                                        ccgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCa                                                                                                        <  2:2601678/91‑1 (MQ=255)

                                                                                                         cgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGcccc                                            >  1:2803021/1‑150 (MQ=255)

                                                                                                         cgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGcccc                                            <  2:1864530/150‑1 (MQ=255)

                                                                                                          gcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCt                                           >  2:1329979/1‑150 (MQ=255)

                                                                                                               tGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCgg                                      >  2:3909985/1‑150 (MQ=255)

                                                                                                                  gccacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                                   >  2:2206550/8‑150 (MQ=25)

                                                                                                                  gcAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                                   <  1:4396325/150‑1 (MQ=255)

                                                                                                                    cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                                 >  2:3458151/6‑150 (MQ=25)

                                                                                                                    aaGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                                 <  2:4634319/150‑1 (MQ=255)

                                                                                                                    aaGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                                 >  2:4788489/1‑150 (MQ=255)

                                                                                                                       cACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGt                              <  1:5090722/150‑1 (MQ=255)

                                                                                                                        aCCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGtt                             >  1:2811484/1‑150 (MQ=255)

                                                                                                                        aCCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGtt                             >  1:1673641/1‑150 (MQ=255)

                                                                                                                         ccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                      <  1:2562005/108‑1 (MQ=21)

                                                                                                                         ccGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTc                            <  1:2336661/150‑1 (MQ=255)

                                                                                                                         ccGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTc                            <  1:1309089/150‑1 (MQ=255)

                                                                                                                           ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                          <  1:4831957/150‑1 (MQ=255)

                                                                                                                           ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                          <  2:3809560/150‑1 (MQ=255)

                                                                                                                           ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                          <  2:1265133/150‑1 (MQ=255)

                                                                                                                                      ctgcccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCAcagc               <  1:3948304/150‑1 (MQ=255)

                                                                                                                                       tgcccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCAcagca              <  2:2811484/150‑1 (MQ=255)

                                                                                                                                         cccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGc                                               <  2:663732/115‑1 (MQ=255)

                                                                                                                                          ccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGAt           >  1:3373920/1‑150 (MQ=255)

                                                                                                                                          ccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGAt           >  1:665034/1‑150 (MQ=255)

                                                                                                                                           caggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAg                                                                 <  1:592391/95‑1 (MQ=255)

                                                                                                                                            aggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                          <  2:4748909/133‑1 (MQ=255)

                                                                                                                                            aggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                          >  1:4748909/1‑133 (MQ=255)

                                                                                                                                             ggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATcca        <  2:4558929/150‑1 (MQ=255)

                                                                                                                                                   tgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGc                                                                <  1:2435034/88‑1 (MQ=255)

                                                                                                                                                   tgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGc                                                                >  2:2435034/1‑88 (MQ=255)

                                                                                                                                                   tgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATccaggccca  >  1:1595877/1‑150 (MQ=255)

                                                                                                                                                   tgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATccaggccca  >  1:1036508/1‑150 (MQ=255)

CGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATCCAGGCCCA  >  NC_002947/4882275‑4882571

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