BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,244,375	C→G	6.9%	G39A (GGT→GCT) ‡	rnt ←	ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,244,375	0	C	G	6.9%	316.6 / 12.2	102	G39A (GGT→GCT) ‡	rnt	ribonuclease T
Reads supporting (aligned to +/- strand): ref base C (48/46); new base G (2/5); total (51/51)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 109 total aligned reads displayed.

TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTT  >  NC_002947/1244228‑1244524

tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCg                                                                                                                                                     >  1:1127414/1‑150 (MQ=255)

tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCg                                                                                                                                                     >  1:2036523/1‑150 (MQ=255)

tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCg                                                                                                                                                     >  1:2033532/1‑150 (MQ=255)

tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCg                                                                                                                                                     <  1:1208606/150‑1 (MQ=255)

   aGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGGTGAAACCACCGGtt                                                                                                                                                  <  1:1662159/150‑1 (MQ=255)

     tGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACAACCGGTTTc                                                                                                                                                <  2:1355189/150‑1 (MQ=255)

         gccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGTGAACAGGAAACCTTTCTCGTCCATGCCAATCGTGACGGCGGAAACTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCGGTGGTTTTAACa                                                                                                                                            <  1:51949/150‑1 (MQ=255)

         gccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCGGTGGTTTCAACa                                                                                                                                            <  1:2375741/150‑1 (MQ=255)

         gccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCCGTTGAAACCACCGGTTTCAACa                                                                                                                                            <  2:1040739/150‑1 (MQ=255)

         gccggcTCGATGTTGGCGCCCCCCCACGGCTCCACCCGGGTGAAGTATGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCGGTGGTTTCAACa                                                                                                                                            <  1:2375748/150‑1 (MQ=255)

          ccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACAt                                                                                                                                           <  1:2051643/150‑1 (MQ=255)

           cgggTCGATGTTGGCGCCCTTGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCGGTGGTTTCAACagc                                                                                                                                          <  2:792021/150‑3 (MQ=255)

              cTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATccac                                                                                                                                       >  2:2350962/1‑150 (MQ=255)

                 gATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCcaccac                                                                                                                                    >  2:545679/1‑150 (MQ=255)

                        gcgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGc                                                                                                                             <  2:194612/150‑1 (MQ=255)

                         cgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCa                                                                                                                            <  2:424693/150‑1 (MQ=255)

                         cgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCa                                                                                                                            <  1:1163399/150‑1 (MQ=255)

                         cgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCa                                                                                                                            <  2:1438459/150‑1 (MQ=255)

                         cgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCa                                                                                                                            <  1:1180154/150‑1 (MQ=255)

                           cccTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGa                                                                                                                          >  2:1429384/1‑150 (MQ=255)

                             cTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCCCCACGACAGGCAAATa                                                                                                                        >  2:151/1‑150 (MQ=255)

                             cTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATa                                                                                                                        >  2:248960/1‑150 (MQ=255)

                              tCGAACGGCTCCACTCGGTGGAAGTAGGTATGCTCGGGGAACAGGAAACCTTTCTAGTAGCTCCCAATGGTGAGAATAACCGCCTGGTG‑AATGCTCGCCATAACGCTGTTGAAACCACCGGTTTCAACATCCACAACGACCGGCAGATAg                                                                                                                       >  1:544325/1‑150 (MQ=255)

                                 aaCGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGa                                                                                                                          >  1:1607817/1‑144 (MQ=255)

                                 aaCGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGa                                                                                                                          <  2:1607817/144‑1 (MQ=255)

                                  aCGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGACAATGTTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCcgc                                                                                                                   <  2:72826/150‑1 (MQ=255)

                                      cTCCACTCGGTGGAAGTATGTATGTTCGCTGAACAGGATTACTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                               <  1:379/150‑1 (MQ=255)

                                      cTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                               >  1:1166744/1‑150 (MQ=255)

                                      cTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                               <  2:540768/150‑1 (MQ=255)

                                             cGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGTAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCg                                                                                                        <  2:670814/150‑1 (MQ=255)

                                             cGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCg                                                                                                        <  2:2036523/150‑1 (MQ=255)

                                              ggTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCgg                                                                                                       >  1:510321/1‑150 (MQ=255)

                                                    aGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATccacc                                                                                                                                      <  2:1312621/113‑1 (MQ=255)

                                                    aGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATccacc                                                                                                                                      >  1:1312621/1‑113 (MQ=255)

                                                          tATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGc                                                                                           <  1:171189/150‑1 (MQ=255)

                                                          gaTGGTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCCTCGGGTGGc                                                                                           >  2:1543503/2‑150 (MQ=255)

                                                                    gAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCg                                                                                                        <  2:1214617/127‑1 (MQ=255)

                                                                    gAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCg                                                                                                        >  1:1214617/1‑127 (MQ=255)

                                                                     aaCAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTACTGAAACCACCGGTTTCAACATCCTCCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgctg                                                                                <  1:1334894/150‑1 (MQ=255)

                                                                        aGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCGGTGGTTTCaa                                                                                                                                              <  1:693290/85‑1 (MQ=255)

                                                                        aGGAAACCTTTCTCGTACATGCCAATGGTGACGGCGGCAATTTCCAGGAATGC‑ATCGGTAGCGCTGTTGCAACCACCGGTTTCAACATCCACCACGACCGACAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACAGctgctg                                                                             >  2:1227131/1‑150 (MQ=255)

                                                                          ggaaCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGcc                                                                           >  1:1901750/3‑150 (MQ=255)

                                                                          gAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAg                                                                                                                       <  2:319302/106‑1 (MQ=255)

                                                                          gAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAg                                                                                                                       >  1:319302/1‑106 (MQ=255)

                                                                           aaaCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCt                                                                          <  2:1485192/150‑1 (MQ=255)

                                                                            aaCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCtt                                                                         <  1:545679/150‑1 (MQ=255)

                                                                             aCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTg                                                                        >  2:1993906/1‑150 (MQ=255)

                                                                              ccTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCAGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                       >  2:755336/1‑150 (MQ=255)

                                                                               cTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTggg                                                                      <  1:1765755/150‑1 (MQ=255)

                                                                                      gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTTAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCAACCTTCAGGGGGCGAGGACCGCTGCTGCCTTGGGTGTCCAg                                                               >  2:1379908/1‑150 (MQ=255)

                                                                                      gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAg                                                               >  1:1140280/1‑150 (MQ=255)

                                                                                      gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAg                                                               >  2:859754/1‑150 (MQ=255)

                                                                                      gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAg                                                               <  2:1127414/150‑1 (MQ=255)

                                                                                      gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAg                                                               >  2:2257569/1‑150 (MQ=255)

                                                                                       tCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTc                                                                                                         <  2:1412446/107‑1 (MQ=255)

                                                                                       tCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTc                                                                                                         >  1:1412446/1‑107 (MQ=255)

                                                                                            gCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGCCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtcgtc                                                         <  1:1379908/150‑1 (MQ=255)

                                                                                               aaTGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCCGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                      <  1:1993906/150‑1 (MQ=255)

                                                                                                   gTGACAGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTaa                                                  >  1:343072/1‑150 (MQ=255)

                                                                                                     gACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAg                                                <  1:2350962/150‑1 (MQ=255)

                                                                                                     cgcgtcggtAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGTCCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAg                                                <  2:1306572/145‑1 (MQ=255)

                                                                                                         gcggcAATTTCCAGCAATGC‑ATCGGTAGCGCTGATGATACCACCGGTTTCAACATCCACCACGACCGGCAGATAGACGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCt                                            >  1:789120/1‑150 (MQ=255)

                                                                                                         gcggcAATTTCCAGCAATGC‑ATCGGTAGCGCTGATGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGGGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCt                                            >  1:789191/1‑150 (MQ=255)

                                                                                                          cggcAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCtt                                           <  1:2257569/150‑1 (MQ=255)

                                                                                                           ggcAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTc                                          >  1:1120957/1‑150 (MQ=255)

                                                                                                              aaTTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCt                                       >  1:1925263/1‑150 (MQ=255)

                                                                                                               aTTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTgg                                                                                            >  1:801351/1‑96 (MQ=255)

                                                                                                               aTTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTgg                                                                                            <  2:801351/96‑1 (MQ=255)

                                                                                                                   ccAGCAATGC‑ATCGGTAGTGATGCTGAAACCTCCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCAcgc                                  <  2:1099354/150‑1 (MQ=255)

                                                                                                                    cAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCtt                                                                         >  2:1166473/1‑110 (MQ=255)

                                                                                                                    cAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCtt                                                                         <  1:1166473/110‑1 (MQ=255)

                                                                                                                      gCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGtt                               <  1:79878/150‑1 (MQ=255)

                                                                                                                              aTCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTcc                        >  2:76336/1‑150 (MQ=255)

                                                                                                                              aTCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTcc                        >  1:1692159/1‑150 (MQ=255)

                                                                                                                               tCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCca                       <  1:2445314/150‑1 (MQ=255)

                                                                                                                               ccGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCca                       <  2:544325/149‑1 (MQ=255)

                                                                                                                                  gTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCcagca                    <  1:1009408/150‑1 (MQ=255)

                                                                                                                                  gTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCcagca                    >  2:2402399/1‑150 (MQ=255)

                                                                                                                                     gcgcTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGgc                 <  1:859754/150‑1 (MQ=255)

                                                                                                                                        cTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGTTCTTCGCTCACGCGTTTTTCTCCAGCAGgcgac              >  1:401645/1‑148 (MQ=255)

                                                                                                                                        cTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGcc              >  2:1625590/1‑150 (MQ=255)

                                                                                                                                         tGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCa             >  1:530723/1‑150 (MQ=255)

                                                                                                                                          gTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAg            >  1:2418188/1‑150 (MQ=255)

                                                                                                                                          gTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtc                                                            >  2:1384167/1‑102 (MQ=255)

                                                                                                                                          gTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtc                                                            <  1:1384167/102‑1 (MQ=255)

                                                                                                                                           ttGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTg                                                                        <  1:1107854/89‑1 (MQ=255)

                                                                                                                                           ttGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTg                                                                        >  2:1107854/1‑89 (MQ=255)

                                                                                                                                           ttGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgc           >  1:2087510/1‑150 (MQ=255)

                                                                                                                                           ttGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgc           >  1:1650056/1‑150 (MQ=255)

                                                                                                                                              aaaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCa        >  2:1955497/1‑150 (MQ=255)

                                                                                                                                              aaaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCa        >  1:287028/1‑150 (MQ=255)

                                                                                                                                               aaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                      >  2:2276734/1‑103 (MQ=255)

                                                                                                                                               aaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                      <  1:2276734/103‑1 (MQ=255)

                                                                                                                                               aaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAg       <  1:2174742/150‑1 (MQ=255)

                                                                                                                                                accaccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCCTCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                       >  2:1285038/1‑85 (MQ=255)

                                                                                                                                                accaccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                       <  1:1285038/85‑1 (MQ=255)

                                                                                                                                                accaccGGTTTCAACATCCACCACGACCGGCAGATAGCAGCTGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGt      <  1:21895/150‑1 (MQ=255)

                                                                                                                                                  caccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGttt    >  2:953328/1‑150 (MQ=255)

                                                                                                                                                   accGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGtttt   >  1:486527/1‑150 (MQ=255)

                                                                                                                                                    ccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGttttt  >  2:1283213/1‑150 (MQ=255)

TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTT  >  NC_002947/1244228‑1244524

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