BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,272,817	G→A	5.8%	V23I (GTC→ATC)	PP_3743 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,272,817	0	G	A	5.8%	359.8 / 8.6	121	V23I (GTC→ATC)	PP_3743	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (59/54); new base A (2/5); total (61/60)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 121 total aligned reads displayed.

AAGGGAACTACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTATGTG  >  NC_002947/4272669‑4272960

aaGGGAACTACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGt                                                                                                                                                <  1:688694/150‑1 (MQ=255)

 aGGGAACTACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTCTTTTACTATGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTc                                                                                                                                               >  1:1240622/1‑150 (MQ=255)

       cTACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCg                                                                                                                                         >  2:518954/1‑150 (MQ=255)

        tACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTc                                                                                                                                               <  1:34655/143‑1 (MQ=255)

        tACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTc                                                                                                                                               >  2:34655/1‑143 (MQ=255)

          ccccAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGc                                                                                                                                      >  1:1384219/1‑150 (MQ=255)

           cccAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCg                                                                                                                                     >  1:898309/1‑150 (MQ=255)

           cccAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCg                                                                                                                                     >  1:1102391/1‑150 (MQ=255)

             cAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCATCGACGCGGGccag                                                                                                                                   <  1:498908/150‑4 (MQ=255)

                 gCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTCGCCCGCATCGACGCGGGccagggca                                                                                                                               <  2:1950916/150‑8 (MQ=255)

                 gCCCCCGCCGTCTTAGATACAGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTCGCCCGCATCGACGCGGGccagggca                                                                                                                               <  2:771675/150‑8 (MQ=255)

                   ccccGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCACGTCGATGCGGGCGACGTCGcc                                                                                                                             >  2:1199252/1‑150 (MQ=255)

                             ttAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCg                                                                                                                   >  2:357699/1‑150 (MQ=255)

                               aGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGAt                                                                                                                                            <  2:1260287/123‑1 (MQ=255)

                                   aCTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGt                                                                                                                                 >  1:276813/1‑130 (MQ=255)

                                   aCTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGt                                                                                                                                 <  2:276813/130‑1 (MQ=255)

                                       aGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGc                                                                                                         >  1:1008688/1‑150 (MQ=255)

                                       aGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGc                                                                                                         <  2:1708316/150‑1 (MQ=255)

                                         tGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGcc                                                                                                                    <  1:1338538/137‑1 (MQ=255)

                                           gTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCg                                                                                                     >  2:1057344/1‑150 (MQ=255)

                                           gTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCCGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCg                                                                                                     >  2:771423/1‑150 (MQ=255)

                                                tttCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTa                                                                                                <  1:554928/150‑1 (MQ=255)

                                                      tgtgCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                          <  1:2057068/150‑1 (MQ=255)

                                                      tgtgCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGAAAATTCTGTCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                          <  2:1797473/150‑1 (MQ=255)

                                                      tgtgCTGCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                          >  1:772956/1‑150 (MQ=255)

                                                           ttCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCgg                                                                                                                                        >  2:41368/1‑99 (MQ=255)

                                                           ttCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCgg                                                                                                                                        <  1:41368/99‑1 (MQ=255)

                                                           ttCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCa                                                                                     >  1:1472568/1‑150 (MQ=255)

                                                                gCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgac                                                                                >  1:1762235/1‑150 (MQ=255)

                                                                gCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCCGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgac                                                                                <  1:771423/150‑1 (MQ=255)

                                                                 cTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgacg                                                                               >  1:1509772/1‑150 (MQ=255)

                                                                   aCTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgacgac                                                                             <  1:1837662/150‑1 (MQ=255)

                                                                     tGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGc                                                                                                                                          >  1:2292379/1‑87 (MQ=255)

                                                                     tGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGc                                                                                                                                          <  2:2292379/87‑1 (MQ=255)

                                                                     tGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACaa                                                                           >  2:105859/1‑150 (MQ=255)

                                                                      ggAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGt                                                                                                                                                >  1:955312/1‑80 (MQ=255)

                                                                       gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGTTGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAAc                                                                         <  1:743365/150‑1 (MQ=255)

                                                                       gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTc                                                                                                      <  1:1950352/121‑1 (MQ=255)

                                                                       gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTc                                                                                                      >  2:1950352/1‑121 (MQ=255)

                                                                       gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                          >  1:1387568/1‑133 (MQ=255)

                                                                       gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                          <  2:1387568/133‑1 (MQ=255)

                                                                       gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAAc                                                                         <  2:1762235/150‑1 (MQ=255)

                                                                           acacGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGat                                                                     <  2:1310780/150‑1 (MQ=255)

                                                                           acacGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTACTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGat                                                                     <  2:1596799/150‑1 (MQ=255)

                                                                            cacGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGatc                                                                    >  1:513651/1‑150 (MQ=255)

                                                                            cacGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGatc                                                                    <  1:583117/150‑1 (MQ=255)

                                                                            cacGACATGCGCCACCTACTGCTCGCACCTGCGCAACTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGatc                                                                    >  1:789021/1‑150 (MQ=255)

                                                                             acGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGatca                                                                   >  2:2373984/1‑150 (MQ=255)

                                                                               gACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGatcatc                                                                 <  1:2373984/150‑1 (MQ=255)

                                                                               gACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTCGCCCGCATCGACGCGGGccag                                                                                                                                   >  2:1556131/1‑81 (MQ=38)

                                                                               gACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTCGCCCGCATCGACGCGGGccag                                                                                                                                   <  1:1556131/84‑4 (MQ=38)

                                                                                 cATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCaa                                                                                    <  1:1247339/129‑1 (MQ=255)

                                                                                 cATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATACGACCTTCaa                                                                                    >  2:1247339/1‑129 (MQ=255)

                                                                                      gcCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGCTGATCCCTGccc                                                          >  1:1566710/1‑150 (MQ=255)

                                                                                       ccACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCg                                                         >  2:1290536/1‑150 (MQ=255)

                                                                                        cACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGt                                                        >  2:1437881/1‑150 (MQ=255)

                                                                                          ccTACTGCTCGCCCCTTCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGACACCACGGCAGCCATTCCCCCTTCGGTGTTTTCCGTCGTCAAGTACGACAAACGGATCATCCCTGCCCGTGa                                                      <  1:1749215/150‑1 (MQ=255)

                                                                                          ccTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgac                                                                                <  2:2132374/124‑1 (MQ=255)

                                                                                          ccTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGa                                                      <  2:1472568/150‑1 (MQ=255)

                                                                                           cTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGAc                                                     <  1:1057344/150‑1 (MQ=255)

                                                                                              cTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                          <  2:1340987/110‑1 (MQ=255)

                                                                                              cTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                          >  1:1340987/1‑110 (MQ=255)

                                                                                                gCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCTTCGATGCCGGCGACGTCGCCACCTCGACCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAACCGGATCATCCCTGCCCGTGACGAACt                                                <  2:898309/150‑1 (MQ=255)

                                                                                                 cTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGa                                                                      >  1:817472/1‑127 (MQ=255)

                                                                                                 cTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGa                                                                      >  2:283441/1‑127 (MQ=255)

                                                                                                 cTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGa                                                                      <  2:817472/127‑1 (MQ=255)

                                                                                                  tCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACttt                                              >  2:2005562/1‑150 (MQ=255)

                                                                                                    gCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCa                                            <  2:1384219/150‑1 (MQ=255)

                                                                                                      cccTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGga                                                                                 <  2:448123/111‑1 (MQ=255)

                                                                                                      cccTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGc                                          <  2:1743562/150‑1 (MQ=255)

                                                                                                          gcgcTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTtc                                      <  2:405006/150‑1 (MQ=255)

                                                                                                             ctgctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAg                                                                                   >  2:2030963/1‑102 (MQ=255)

                                                                                                             ctgctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAg                                                                                   <  1:2030963/102‑1 (MQ=255)

                                                                                                             ctgctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTtcgtc                                   <  1:1544953/150‑1 (MQ=255)

                                                                                                               gctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGc                                 <  1:518954/150‑1 (MQ=255)

                                                                                                               gctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGc                                 <  2:897336/150‑1 (MQ=255)

                                                                                                                 tgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCa                               >  1:801144/1‑150 (MQ=255)

                                                                                                                  gctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAg                              >  2:2088464/1‑150 (MQ=255)

                                                                                                                   ctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAgc                             >  2:1772469/1‑150 (MQ=255)

                                                                                                                      cTGATGCCTGCCGGTGCTGCCCTGGCCCGAGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAgcggc                          >  2:646927/1‑150 (MQ=255)

                                                                                                                       tGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACCAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAgcggcg                         <  2:336063/150‑1 (MQ=255)

                                                                                                                       tGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAgcggcg                         >  1:424777/1‑150 (MQ=255)

                                                                                                                       tGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAgcggcg                         >  1:1722543/1‑150 (MQ=255)

                                                                                                                        gATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAg                                           <  1:785421/131‑1 (MQ=255)

                                                                                                                        gATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAg                                           >  2:785421/1‑131 (MQ=255)

                                                                                                                         aTGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGt                       >  1:882289/1‑150 (MQ=255)

                                                                                                                           gcctgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGc                                                                                                    >  1:2376838/1‑71 (MQ=255)

                                                                                                                           gcctgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGc                                                                                                    <  2:2376838/71‑1 (MQ=255)

                                                                                                                           acctgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTtc                                      <  1:1609914/132‑1 (MQ=255)

                                                                                                                             ctgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCa                   >  2:851096/1‑150 (MQ=255)

                                                                                                                              tgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTANCAt                  <  2:1509772/150‑1 (MQ=255)

                                                                                                                               gccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATc                 <  2:424777/150‑1 (MQ=255)

                                                                                                                               gccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATc                 <  1:357699/150‑1 (MQ=255)

                                                                                                                               gccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATc                 <  1:358021/150‑1 (MQ=255)

                                                                                                                                         cccTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATcc                                                                >  2:775062/1‑93 (MQ=255)

                                                                                                                                         cccTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCt       >  1:2069515/1‑150 (MQ=255)

                                                                                                                                            tcgCCCGCATCGACGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTAtg    >  2:154809/3‑150 (MQ=255)

                                                                                                                                             ggCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTAtgt   <  1:2094968/150‑1 (MQ=255)

                                                                                                                                              gCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTAtgtg  >  1:2202312/1‑150 (MQ=255)

                                                                                                                                              gCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTAtgtg  <  2:1008688/150‑1 (MQ=255)

AAGGGAACTACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTATGTG  >  NC_002947/4272669‑4272960

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