BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,293,270	C→A	6.5%	A289A (GCC→GCA)	PP_1128 →	OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,293,270	0	C	A	6.5%	327.9 / 7.8	105	A289A (GCC→GCA)	PP_1128	OmpA family protein
Reads supporting (aligned to +/- strand): ref base C (47/51); new base A (3/4); total (50/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 105 total aligned reads displayed.

CCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGCCAAC  >  NC_002947/1293123‑1293417

ccGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGCTGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGc                                                                                                                                                   >  1:519191/1‑150 (MQ=255)

  gcgcCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCtg                                                                                                                                                 >  1:978653/1‑150 (MQ=255)

   cgcCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCtgg                                                                                                                                                <  1:1162389/150‑1 (MQ=255)

   cgcCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCtgg                                                                                                                                                <  1:1455536/150‑1 (MQ=255)

    gcCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCGGCCTTGGCAGCtggt                                                                                                                                               <  1:2257427/150‑1 (MQ=255)

     cccTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCtggtg                                                                                                                                              <  1:1780345/150‑1 (MQ=255)

      ccTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCtggtgg                                                                                                                                             <  2:2114677/150‑1 (MQ=255)

       cTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCGGCCTTGGCAGCtggtgcc                                                                                                                                            <  1:819027/150‑3 (MQ=255)

        tGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCa                                                                                                                                           >  1:1733869/1‑150 (MQ=255)

         gAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCaa                                                                                                                                          >  2:1913804/1‑150 (MQ=255)

            ccgccgCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAAcc                                                                                                                                       >  2:1459528/1‑150 (MQ=255)

                cgcgTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGcc                                                                                                                                   <  2:2215198/150‑1 (MQ=255)

                     ggAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCaa                                                                                                                              <  1:1123104/150‑1 (MQ=255)

                     ggAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCaa                                                                                                                              <  2:978653/150‑1 (MQ=255)

                     ggAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCaa                                                                                                                              <  2:1024633/150‑1 (MQ=255)

                      gAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAg                                                                                                                             >  1:1035389/1‑150 (MQ=255)

                      gAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAg                                                                                                                             <  2:970710/150‑1 (MQ=255)

                      gAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAg                                                                                                                             <  2:503849/150‑1 (MQ=255)

                       aaaTGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCGGCCTTGGCAGCtggtgccttg                                                                                                                                         >  1:1766683/1‑132 (MQ=255)

                       aaaTGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCGGCCTTGGCAGCtggtgccttg                                                                                                                                         <  2:1766683/137‑6 (MQ=255)

                          tGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGACCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGc                                                                                                                           <  2:1523750/148‑1 (MQ=255)

                          tGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGACCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGc                                                                                                                           >  1:1523750/1‑148 (MQ=255)

                           gTTGCTGACCCCGCAAACCAACATGGCCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTgc                                                                                                                        <  2:2360803/150‑1 (MQ=255)

                           gTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTgc                                                                                                                        >  2:1855447/1‑150 (MQ=255)

                            ttGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCCAGGCCGCTGGTGGCAAACCTCACGTCAAGGCTgcg                                                                                                                       >  1:496241/1‑150 (MQ=255)

                             tGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTgcgg                                                                                                                      <  2:2374432/150‑1 (MQ=255)

                              gCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTgcggc                                                                                                                     >  1:1350239/1‑150 (MQ=255)

                                tGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCga                                                                                                                   >  1:1523778/1‑150 (MQ=255)

                                 gACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATCCCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCgaa                                                                                                                  <  1:632779/150‑1 (MQ=255)

                                 gACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTTAAGGCTGCGGCgaa                                                                                                                  <  2:1639991/150‑1 (MQ=255)

                                    cccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCtggt                                                                                                                                               >  2:1823551/1‑118 (MQ=255)

                                    cccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCtggt                                                                                                                                               <  1:1823551/118‑1 (MQ=255)

                                     ccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCAGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGaaa                                                                                                              >  2:1790464/1‑150 (MQ=255)

                                      cGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTCCCGTCAGGCCTGAGGCGAAGAAAc                                                                                                             >  2:1942474/1‑150 (MQ=255)

                                         aaaaCACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTg                                                                                                          <  2:1035389/150‑1 (MQ=255)

                                           aaCACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTa                                                                                                        <  1:753474/150‑1 (MQ=255)

                                            acacCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCCGTGCCGGATG‑CAACGCACCAAGCTGCCAACGCCGGTGGTGGTAAACCTGGGGTCAAGTCTGCGGCGAAGAAACCTGTAg                                                                                                       <  2:636121/150‑1 (MQ=255)

                                              acCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGcc                                                                                                     >  2:1127763/1‑150 (MQ=255)

                                              acCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGcc                                                                                                     >  1:1585069/1‑150 (MQ=255)

                                                 aTGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAg                                                                                                  <  1:387805/150‑1 (MQ=255)

                                                  tGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAgg                                                                                                 >  2:1094925/1‑150 (MQ=255)

                                                             ctgGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCt                                                                                      >  2:476717/1‑150 (MQ=255)

                                                               ggCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGc                                                                                    >  1:2002139/1‑150 (MQ=255)

                                                                 ccAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGACGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGACCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGa                                                                                  >  2:755750/1‑150 (MQ=255)

                                                                 acAAATACCAGCAGGCCGCACCGGCACCCGCCCCGGCCTCGATAGCGGCGGTGGAGGATGCGAAGGCACC‑AGCTGCCAGTGCCGCTGGTGGCACACCTGCCTTCACGGCTCCGACGAAGAACCCTGTAGCCACGGCCAAAACGCCTGCGa                                                                                  >  2:464042/2‑150 (MQ=255)

                                                                  cAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGTGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGaa                                                                                 >  2:604428/1‑150 (MQ=255)

                                                                  cAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAAc                                                                                                                                        <  1:1280501/95‑1 (MQ=255)

                                                                  cAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAAc                                                                                                                                        >  2:1280501/1‑95 (MQ=255)

                                                                  cAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAGGAAACCTGTAGCCAAGGCCAAAACGCCTGCGaa                                                                                 >  1:725890/1‑150 (MQ=255)

                                                                   aaaTACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAAc                                                                                                             >  1:2119242/1‑121 (MQ=255)

                                                                   aaaTACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAAc                                                                                                             <  2:2119242/121‑1 (MQ=255)

                                                                    aaTACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGa                                                                               <  1:1127763/150‑1 (MQ=255)

                                                                       aCCAGCAGGCCGCCCCGGCACCCGCCACGGCCGCGAGTGTGGCGGTGCAGGATGCCAAGGCACA‑AGATGGCGTGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTg                                                                            <  1:1852759/150‑1 (MQ=255)

                                                                         cagcagGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCt                                                                                                           >  2:1877872/1‑117 (MQ=255)

                                                                         cagcagGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCt                                                                                                           <  1:1877872/117‑1 (MQ=255)

                                                                              ggCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCaa                                                                     <  1:1754846/150‑1 (MQ=255)

                                                                              ggCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCaa                                                                     <  2:29999/150‑1 (MQ=255)

                                                                              ggCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCaa                                                                     <  2:836726/150‑1 (MQ=255)

                                                                               gcctccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAg                                                                    >  2:837390/1‑150 (MQ=255)

                                                                               gccgccCCGGCACGCGCCCCGGCCTCGATGGGGGCGGTTCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGAAAACATGCCGACAAGGCTGTGGCGGAGAAAACTGTATCCAACGCCAAAACGCCTGCGAAGACTGCAGCCAAg                                                                    >  2:2330118/1‑150 (MQ=255)

                                                                               gccgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCt                                                                                      >  2:2398680/1‑132 (MQ=255)

                                                                               gccgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCt                                                                                      <  1:2398680/132‑1 (MQ=255)

                                                                               gccgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAg                                                                    >  2:1660396/1‑150 (MQ=255)

                                                                                   ccccGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGCCGAAGAAACCTGTCGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGaaaa                                                                <  1:1459528/150‑1 (MQ=255)

                                                                                    cccGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGaaaaa                                                               >  1:1181302/1‑150 (MQ=255)

                                                                                     ccGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGaaaaaa                                                              <  2:1526687/150‑1 (MQ=255)

                                                                                        gCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCt                                                           >  2:1508438/1‑150 (MQ=255)

                                                                                          aCCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCAAACGTCGCTGGTTGCAAACCTTCCGTCCAGGCTGCGGCGAAGAAACCTGTATCCAAGTCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTgc                                                         <  1:1855447/150‑1 (MQ=255)

                                                                                                cccGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCAGCGAAGAAACCTGTAGCCAAGGACAAAACGCCTGCGAAGACTGCAGCCACGAAAAAAGCTGTGCCggg                                                   >  2:782789/1‑149 (MQ=255)

                                                                                                    gCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAg                                               >  1:1520493/1‑150 (MQ=255)

                                                                                                      cTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGcc                                             >  1:2104656/1‑150 (MQ=255)

                                                                                                          aTGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCAGAGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTgcc                                         >  1:900314/1‑150 (MQ=255)

                                                                                                                    tGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACGTGTAGCCAAGGCCAAAACGCCTGCGAAAAGGGCGGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCGGCCAAGAAAg                               >  1:263477/1‑150 (MQ=255)

                                                                                                                     gCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGc                              >  1:660657/1‑150 (MQ=255)

                                                                                                                       aGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTg                            <  2:1733869/150‑1 (MQ=255)

                                                                                                                       aGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTg                            >  1:297244/1‑150 (MQ=255)

                                                                                                                             ccAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGGACTTCTTCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTg                      >  1:1171991/1‑150 (MQ=255)

                                                                                                                              tttgccACC‑AGCGGCCTTGGCAGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGtt                       >  1:439993/6‑148 (MQ=255)

                                                                                                                              tttgccACC‑AGCGGCCTTGGCAGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGtt                       <  2:439993/143‑1 (MQ=255)

                                                                                                                              cacggCACC‑AGCTCCCAAGGCCGCTGGTGGCAAACCCGCCGTCACGGCCGAGCTCCGAAAATATAGGGTTGAGGCCAAAACGACTGCCATTACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGa                     <  1:782789/147‑1 (MQ=255)

                                                                                                                              cAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTa                                                                                                        >  1:308096/1‑67 (MQ=255)

                                                                                                                              cAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTa                                                                                                        <  2:308096/67‑1 (MQ=255)

                                                                                                                              cAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGc                                              >  2:301647/1‑125 (MQ=255)

                                                                                                                              cAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGc                                              <  1:301647/125‑1 (MQ=255)

                                                                                                                              cAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGa                     <  1:604428/150‑1 (MQ=255)

                                                                                                                               aaGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGAc                    <  2:1585069/150‑1 (MQ=255)

                                                                                                                               aaGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGAc                    >  1:884547/1‑150 (MQ=255)

                                                                                                                                aGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACa                   <  1:2240728/150‑1 (MQ=255)

                                                                                                                                    aCC‑AGCGGCCTTGGCAGCGGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGAc                                                                              >  2:1296810/1‑87 (MQ=255)

                                                                                                                                     cc‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGaaa              <  1:1921331/150‑1 (MQ=255)

                                                                                                                                      c‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAg             >  1:1675246/1‑150 (MQ=255)

                                                                                                                                      c‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAg             <  2:1181302/150‑1 (MQ=255)

                                                                                                                                         gCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTa           <  1:912413/150‑1 (MQ=255)

                                                                                                                                            gCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCt        <  1:1790464/150‑1 (MQ=255)

                                                                                                                                            gCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCt        <  2:1523778/150‑1 (MQ=255)

                                                                                                                                            gCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTCCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCt        <  2:1350239/150‑1 (MQ=255)

                                                                                                                                             ccAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGaaa                                >  2:320969/1‑125 (MQ=255)

                                                                                                                                             ccAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGaaa                                <  1:320969/125‑1 (MQ=255)

                                                                                                                                              cAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGc      <  1:2330118/150‑1 (MQ=255)

                                                                                                                                                  gCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGCCAAc  >  1:1286648/1‑150 (MQ=255)

CCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACC‑AGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGCCAAC  >  NC_002947/1293123‑1293417

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`