BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,244,374	A→G	7.1%	G39G (GGT→GGC) ‡	rnt ←	ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,244,374	0	A	G	7.1%	277.8 / 7.7	102	G39G (GGT→GGC) ‡	rnt	ribonuclease T
Reads supporting (aligned to +/- strand): ref base A (49/44); new base G (2/5); total (51/51)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 108 total aligned reads displayed.

TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTT  >  NC_002947/1244228‑1244523

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

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

tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCg                                                                                                                                                    >  1:1127414/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:1438459/150‑1 (MQ=255)

                         cgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCa                                                                                                                           <  1:1163399/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)

                                                          tATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGc                                                                                          <  1:171189/150‑1 (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)

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

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

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

                                                                           aaaCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCt                                                                         <  2:1485192/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:2257569/1‑150 (MQ=255)

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

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

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

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

                                                                                            gCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtcgtc                                                        <  1:496376/150‑1 (MQ=255)

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

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

                                                                                               aaTGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                     <  2:1826755/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)

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

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

                                                                                                                tttCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCa                                    >  1:130127/1‑150 (MQ=255)

                                                                                                                  tCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCAcg                                  >  2:1352812/1‑150 (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)

                                                                                                                              ctCGGTCGCGCTGGTGAAACCTCCGGGTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTcc                       <  2:510321/149‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                  gTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCcagca                   <  1:1009408/150‑1 (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)

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

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

                                                                                                                                           ttGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgc          >  2:1032589/1‑150 (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)

                                                                                                                                              aaaCCACCAGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCa       >  1:2313880/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)

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

TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGC‑ATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTT  >  NC_002947/1244228‑1244523

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