BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,312,886	T→G	5.4%	R5R (AGG→CGG)	PP_2034 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,312,886	0	T	G	5.4%	316.2 / 5.2	109	R5R (AGG→CGG)	PP_2034	membrane protein
Reads supporting (aligned to +/- strand): ref base T (53/50); new base G (4/2); total (57/52)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 109 total aligned reads displayed.

CTACGCCCAGCAGGCCGGGCAGCACGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACAGGCTGATGTTGGCCGCC  >  NC_002947/2312739‑2313033

cTACGCCCAGCAGGCCGGGCAGCACGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAAGAGCGTGCGCCCCTgg                                                                                                                                                   >  2:455686/1‑150 (MQ=255)

    gCCCAGCAGGCCGGGCAGCACGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGAACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCa                                                                                                                                               >  2:2327822/1‑150 (MQ=255)

     cccAGCAGGCCGGGCAGCACGGGCATTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCAc                                                                                                                                              >  1:2287616/1‑150 (MQ=255)

      ccAGCAGGCCGGGCAGCACGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACg                                                                                                                                             >  2:1460259/1‑150 (MQ=255)

           agGCCGGGCAGCACGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAg                                                                                                                                        <  1:1331202/150‑1 (MQ=255)

            ggCCGGGCAGCACGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAgg                                                                                                                                       >  1:949729/1‑150 (MQ=255)

             gCCGGGCAGCACGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAggt                                                                                                                                      <  1:2099022/150‑1 (MQ=255)

               cGGGCAGCACGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAggtgg                                                                                                                                    >  2:425418/1‑150 (MQ=255)

                     gcaTGGGCAGTAGCGCCGGGATCAGCATGATGTGCATTGACCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGaaa                                                                                                                              <  1:713236/150‑1 (MQ=255)

                     gcaCGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCAGGGGCGCACgctg                                                                                                                                          <  2:93902/138‑4 (MQ=255)

                     gcaCGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCAGGGGCGCACgctg                                                                                                                                          >  1:93902/1‑135 (MQ=255)

                      caCGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAg                                                                                                                             >  2:2059180/1‑150 (MQ=255)

                           gCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGc                                                                                                                        >  2:2235679/1‑150 (MQ=255)

                           gCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGc                                                                                                                        <  2:340038/150‑1 (MQ=255)

                            cAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCt                                                                                                                       <  1:1254890/150‑1 (MQ=255)

                                  cgcCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCgg                                                                                                                 >  2:1523388/1‑150 (MQ=255)

                                   gcCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCACGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTGGCTTAGAggg                                                                                                                >  2:1056862/1‑150 (MQ=255)

                                         aTCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGgcgc                                                                                                          >  2:291882/1‑150 (MQ=255)

                                         aTCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGgcgc                                                                                                          >  2:2001419/1‑150 (MQ=255)

                                          tCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGcgcc                                                                                                         >  1:2120428/1‑150 (MQ=255)

                                            aggatgatgTACACATGGCTGACGAAGTGGACGGCAGCAATTTGGGTCAGCAGCGTGGGGCCAGCGTGGAGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGcgccgc                                                                                                       <  1:1566525/147‑1 (MQ=255)

                                              cATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACg                                                                                                                                             <  1:1941344/110‑1 (MQ=255)

                                              cATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCgcgc                                                                                                     >  2:1845855/1‑150 (MQ=255)

                                              cATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCAATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCgcgc                                                                                                     >  1:1746049/1‑150 (MQ=255)

                                                    gtgCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCAAGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCg                                                                                               <  1:852915/150‑1 (MQ=255)

                                                      gCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGTGGCATGTCGATACCTGTTGCACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCggg                                                                                             <  1:1915274/150‑1 (MQ=255)

                                                       cacaTGGCTGACGAAGTGCACGGCAGCAATTTGGTTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTc                                                                                                                          <  2:661679/120‑1 (MQ=255)

                                                       cacaTGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTc                                                                                                                          >  1:661679/1‑120 (MQ=255)

                                                        acaTGGCTGACGAAGTGTACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGACACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAATGTGGTTGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCa                                                                                           <  1:2279299/150‑1 (MQ=255)

                                                        acaTGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCa                                                                                           >  2:283348/1‑150 (MQ=255)

                                                         caTGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGgc                                                                                                            >  2:66346/1‑132 (MQ=255)

                                                         caTGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGAACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGgc                                                                                                            <  1:66346/132‑1 (MQ=255)

                                                          aTGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAgg                                                                                         >  2:465456/1‑150 (MQ=255)

                                                           tGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGc                                                                                        <  1:843394/150‑1 (MQ=255)

                                                            ggCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGATGAACAGCGTGCGCCCCTGGCGCACGGAAACGCGGTGAAAAGATCTCTTATCGGGAGTTGCCGCGCAAACCGGGCAGGCg                                                                                       >  2:843394/1‑150 (MQ=255)

                                                             gCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGa                                                                                      >  2:505139/1‑150 (MQ=255)

                                                               tGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGAcc                                                                                    >  2:1226192/1‑150 (MQ=255)

                                                               tGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCAGGGGCGCACgc                                                                                                                                            >  2:727507/1‑93 (MQ=255)

                                                                  cGAAGTGCACGGCAGGAATTTGGGGCAACACGTTGGGGCCAGCGTGGGGCATGTCGATAGCGGTTGAACAGATTGCGACCATTGGGCGCGGGACGGAGCTGGAAAGCTCGCTTAACGGGAGGCACAGCGCAAACCGGGCAGGCGACCGtt                                                                                 >  2:432773/1‑150 (MQ=255)

                                                                   gAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTc                                                                                <  2:544949/150‑1 (MQ=255)

                                                                   gAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGAACATCACGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGCGGCCGCGCAAACCTGGCAGGCGACCGTTc                                                                                >  1:1677147/1‑150 (MQ=255)

                                                                    aaGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTTTGAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCa                                                                               <  1:425418/150‑1 (MQ=255)

                                                                           cGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAAc                                                                        <  2:1440859/150‑1 (MQ=255)

                                                                             gcagcaATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGaaa                                                                                                                              <  1:904970/94‑1 (MQ=255)

                                                                             gcagcaATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGaaa                                                                                                                              >  2:904970/1‑94 (MQ=255)

                                                                              cagcaATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTg                                                                     <  1:505139/150‑1 (MQ=255)

                                                                              cagcaATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTg                                                                     <  1:2327822/150‑1 (MQ=255)

                                                                              cagcaATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTg                                                                     <  2:646078/150‑1 (MQ=255)

                                                                                  aaTTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCTATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCaaa                                                                                                  <  1:2262402/117‑1 (MQ=255)

                                                                                  aaTTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCaaa                                                                                                  >  2:2262402/1‑117 (MQ=255)

                                                                                  aaTTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGcc                                                                 >  1:1264901/1‑150 (MQ=255)

                                                                                   aTTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCca                                                                >  1:1294741/1‑150 (MQ=255)

                                                                                   aTTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCca                                                                <  1:1226192/150‑1 (MQ=255)

                                                                                     ttGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCcagc                                                              <  1:2235679/150‑1 (MQ=255)

                                                                                      tGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGTTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAg                                                                              <  1:280030/133‑1 (MQ=255)

                                                                                         gTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGAc                                                          >  2:160495/1‑150 (MQ=255)

                                                                                          tCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGAcc                                                         >  1:860363/1‑150 (MQ=255)

                                                                                             gcagcTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAAt                                                      <  2:949729/150‑1 (MQ=255)

                                                                                             gcagcTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAAt                                                      <  1:1845855/150‑1 (MQ=255)

                                                                                             gcagcTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGAGCCGCGCAAACCGGGCAGGCGACCGTTCATGACAACGTGGGCCAGCAGACCAAt                                                      <  1:1460259/150‑1 (MQ=255)

                                                                                                 cTTGGGGCCAGCGTGGGGCATGTCGATACCGGTTGAACAGAGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGCCGCCGCGCAAACCGGGCAAGCCCCCGTTCAGGACTACTTGGGCCCGCAGACCCATTAcc                                                  <  2:142085/150‑1 (MQ=255)

                                                                                                  ttGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGgtggtg                                                                                                                                  >  2:2110171/1‑69 (MQ=255)

                                                                                                  ttGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGgtggtg                                                                                                                                  <  1:2110171/69‑1 (MQ=255)

                                                                                                  ttGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCa                                                 <  1:2074165/150‑1 (MQ=255)

                                                                                                    ggagCCAGCGGGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAgg                                               >  2:1625531/4‑150 (MQ=255)

                                                                                                     gggCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGc                                              <  2:2287616/150‑1 (MQ=255)

                                                                                                      ggCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGcc                                             >  1:16275/1‑150 (MQ=255)

                                                                                                       gCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAg                                                                                                                                        <  1:45607/58‑1 (MQ=255)

                                                                                                       gCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAg                                                                                                                                        >  2:45607/1‑58 (MQ=255)

                                                                                                        ccAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGcccc                                           <  1:328293/150‑1 (MQ=255)

                                                                                                          aGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAg                                         >  1:1276189/1‑150 (MQ=255)

                                                                                                            cGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGaa                                       <  2:2120428/150‑1 (MQ=255)

                                                                                                             gTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGaaa                                      <  1:291882/150‑1 (MQ=255)

                                                                                                             gTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGATCGGTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGaaa                                      <  1:146278/150‑1 (MQ=255)

                                                                                                                 ggCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCcagca                                                             <  2:213801/123‑1 (MQ=255)

                                                                                                                 ggCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCcagca                                                             >  1:213801/1‑123 (MQ=255)

                                                                                                                   cATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCg                                <  2:2155689/150‑1 (MQ=255)

                                                                                                                    aTGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGc                               <  1:2059180/150‑1 (MQ=255)

                                                                                                                     tGTCAATACCTGTTGAACAGCGTGCGCCCCGGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAAGGACCAGGCCCCAGAAAGCcccccc                              <  2:1746049/150‑1 (MQ=255)

                                                                                                                      gTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGGTGATTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCa                             >  1:1194499/1‑150 (MQ=255)

                                                                                                                      gTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCa                             >  1:801993/1‑150 (MQ=255)

                                                                                                                       tCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCaa                            >  1:1954081/1‑150 (MQ=255)

                                                                                                                           tACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCCACGCGCGACTCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATAcc                        <  1:283348/150‑1 (MQ=255)

                                                                                                                            aCCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATAccc                       <  2:1581849/150‑1 (MQ=255)

                                                                                                                              cTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTgg                                                                    <  2:572007/103‑1 (MQ=255)

                                                                                                                              cTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCaa                     >  1:1456683/1‑150 (MQ=255)

                                                                                                                                gtggAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCGCGCTTAGCGGGAGGCGCAGCGCAAACCGGGTAGGGGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACa                   <  1:2001419/147‑1 (MQ=255)

                                                                                                                                 ttGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCa                                                                               <  1:524925/89‑1 (MQ=255)

                                                                                                                                    aaCAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAg                                                <  2:243527/117‑1 (MQ=255)

                                                                                                                                    aaCAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACAGGCt               >  1:1665305/1‑150 (MQ=255)

                                                                                                                                      cAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGaaa                                      <  2:421477/125‑1 (MQ=255)

                                                                                                                                      cAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGaaa                                      >  1:421477/1‑125 (MQ=255)

                                                                                                                                       gccGTGCGCCAGGGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCAGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACAGGCTGAt            >  1:361277/3‑150 (MQ=255)

                                                                                                                                         cGTGCGCCAGGGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGa                                                           >  1:1648286/1‑101 (MQ=255)

                                                                                                                                           tGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACAGGCTGATGTTg        <  1:1120430/150‑1 (MQ=255)

                                                                                                                                              gcCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACAGGCTGATGTTGgcc     >  1:1320417/1‑150 (MQ=255)

                                                                                                                                                 ccTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCACTGACCAGGCCCCAGAAAGCCCCGCCAATACCCATCAGGCTGATGTTGgccgcc  <  2:320685/150‑1 (MQ=255)

                                                                                                                                                 ccTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACAGGCTGATGTTGgcc     <  2:565184/147‑1 (MQ=255)

                                                                                                                                                 ccTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACAGGCTGATGTTGgcc     >  1:565184/1‑147 (MQ=255)

                                                                                                                                                 ccTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGAGCAGGCGACCGTTCAGCACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACAGGCTGATGTTgg       <  2:1518098/145‑1 (MQ=255)

CTACGCCCAGCAGGCCGGGCAGCACGGGCAGTAGCGCCGGGATCAGCATGATGTGCACATGGCTGACGAAGTGCACGGCAGCAATTTGGGTCAGCAGCTTGGGGCCAGCGTGGGGCATGTCGATACCTGTTGAACAGCGTGCGCCCCTGGCGCACGGCAAGGTGGTGGAAAGCTCGCTTAGCGGGAGGCGCCGCGCAAACCGGGCAGGCGACCGTTCAGGACAACGTGGGCCAGCAGACCAATGACCAGGCCCCAGAAAGCCCCGCCAATACCCAACAGGCTGATGTTGGCCGCC  >  NC_002947/2312739‑2313033

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