BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,244,376	C→T	5.1%	G39S (GGT→AGT) ‡	rnt ←	ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,244,376	0	C	T	5.1%	555.8 / 11.7	176	G39S (GGT→AGT) ‡	rnt	ribonuclease T
Reads supporting (aligned to +/- strand): ref base C (85/81); new base T (6/3); total (91/85)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTC  >  NC_002947/1244228‑1244525

tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCg                                                                                                                                                      >  2:290198/1‑150 (MQ=255)

tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCg                                                                                                                                                      <  2:2207002/150‑1 (MQ=255)

  aaGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGt                                                                                                                                                    >  1:3403058/1‑150 (MQ=255)

   aGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGtt                                                                                                                                                   >  1:4598210/1‑150 (MQ=255)

    gTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGttt                                                                                                                                                  >  1:224032/1‑150 (MQ=255)

      gCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCa                                                                                                                                                <  2:2487664/150‑1 (MQ=255)

         gccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCGGTGGTTTCAACa                                                                                                                                             <  2:1951587/150‑1 (MQ=255)

         gccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCGGTGGTTTCAACa                                                                                                                                             <  1:5165257/150‑1 (MQ=255)

         gccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACa                                                                                                                                             <  2:3535058/150‑1 (MQ=255)

          ccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACAt                                                                                                                                            <  2:4667712/150‑1 (MQ=255)

            ggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATcc                                                                                                                                          >  1:170432/1‑150 (MQ=255)

                cGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATccacca                                                                                                                                      >  2:3915710/1‑150 (MQ=255)

                 gATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCcaccac                                                                                                                                     <  2:5398153/150‑1 (MQ=255)

                 gATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCcaccac                                                                                                                                     <  2:2171355/150‑1 (MQ=255)

                        gcgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTGGGGGAACAGGTAACCTTTCTGTTCACTGCCAATGGTGACGGCGGCAATTTGCAGAAATGCATCGGTAGCGATGTTGAAACCACCTGTTTCAACATCCACCACGACCGGc                                                                                                                              <  1:5150782/150‑1 (MQ=255)

                         cgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCa                                                                                                                             <  2:532592/150‑1 (MQ=255)

                             cTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATa                                                                                                                         >  1:940188/1‑150 (MQ=255)

                                    ggCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAATCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCgg                                                                                                                  >  2:4255601/1‑150 (MQ=255)

                                    ggCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGGTGAAACCACCGGGTTCAACATCCACCACGACCGGCAGATCGCCGCgg                                                                                                                  >  2:2981886/1‑150 (MQ=255)

                                      cTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                                <  1:4107097/150‑1 (MQ=255)

                                      cTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                                >  1:326482/1‑150 (MQ=255)

                                      cTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                                <  1:1424174/150‑1 (MQ=255)

                                       tCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTTCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAg                                                                                                               <  2:2861657/150‑1 (MQ=255)

                                          aCTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAgcgc                                                                                                            >  2:1280210/1‑150 (MQ=255)

                                            tCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTc                                                                                                          >  1:4479954/1‑150 (MQ=255)

                                            tCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCAGGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTc                                                                                                          >  2:5211164/1‑150 (MQ=255)

                                                      tAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGACGCGCTCGGCCATCggg                                                                                                >  2:1322361/1‑150 (MQ=255)

                                                         gTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCg                                                                                                         >  1:5264418/1‑138 (MQ=255)

                                                         gTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTgg                                                                                             >  1:2624247/1‑150 (MQ=255)

                                                          tATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGc                                                                                            >  2:710763/1‑150 (MQ=255)

                                                             gTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAg                                                                                         <  2:3849822/150‑1 (MQ=255)

                                                              ttCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAgg                                                                                        >  2:880197/1‑150 (MQ=255)

                                                                   ggAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgc                                                                                   >  2:1225374/1‑150 (MQ=255)

                                                                       cAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgctgct                                                                               <  1:719046/150‑1 (MQ=255)

                                                                        aGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCGGTGGTTTCaa                                                                                                                                               >  1:3723106/1‑85 (MQ=255)

                                                                        aGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCGGTGGTTTCaa                                                                                                                                               <  2:3723106/85‑1 (MQ=255)

                                                                         ggAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCtgctgc                                                                             >  1:2593725/1‑150 (MQ=255)

                                                                          gAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCcgcg                                                                                                                   >  2:5164290/1‑111 (MQ=255)

                                                                           aaaCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCt                                                                           >  1:4063169/1‑150 (MQ=255)

                                                                             aCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTg                                                                         <  2:1629826/150‑1 (MQ=255)

                                                                             aCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTg                                                                         >  2:1377162/1‑150 (MQ=255)

                                                                              ccTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                        <  1:4366483/150‑1 (MQ=255)

                                                                              ccTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                        <  1:1621220/150‑1 (MQ=255)

                                                                              ccTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                        <  1:4276034/150‑1 (MQ=255)

                                                                                tttCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCCCCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGgt                                                                      <  2:1849364/150‑1 (MQ=255)

                                                                                tttCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGgt                                                                      <  1:4226321/150‑1 (MQ=255)

                                                                                    tcGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTcc                                                                  <  1:1225374/150‑1 (MQ=255)

                                                                                     cGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCa                                                                 <  2:4129877/150‑1 (MQ=255)

                                                                                     cGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCa                                                                 <  1:5427425/150‑1 (MQ=255)

                                                                                      gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAg                                                                >  1:4279712/1‑150 (MQ=255)

                                                                                         cATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtc                                                             >  2:589041/1‑150 (MQ=255)

                                                                                          aTGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTc                                                                                                                                                 >  1:4987987/1‑65 (MQ=255)

                                                                                            gCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtcgtc                                                          <  1:3922368/150‑1 (MQ=255)

                                                                                              cAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCtt                                                        >  1:615889/1‑150 (MQ=255)

                                                                                               aaTGGTGTCGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                       <  2:918675/150‑1 (MQ=255)

                                                                                               aaTGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                       <  2:4063169/150‑1 (MQ=255)

                                                                                               aaTGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                       <  2:1370882/150‑1 (MQ=255)

                                                                                               aaTGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                       <  2:2107651/150‑1 (MQ=255)

                                                                                                    tGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTaat                                                  >  1:617034/1‑149 (MQ=255)

                                                                                                    tGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTaaa                                                  <  2:2579083/150‑1 (MQ=255)

                                                                                                     gACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAg                                                 <  1:4694412/150‑1 (MQ=255)

                                                                                                          cggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCtt                                            <  1:2049212/150‑1 (MQ=255)

                                                                                                          cggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCtt                                            <  1:787104/150‑1 (MQ=255)

                                                                                                           ggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTc                                           >  2:4548718/1‑150 (MQ=255)

                                                                                                               aTTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTc                                       >  2:2761661/1‑150 (MQ=255)

                                                                                                                tttCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCa                                      >  2:104118/1‑150 (MQ=255)

                                                                                                                tttCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCa                                      >  1:2915042/1‑150 (MQ=255)

                                                                                                                 ttCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCAc                                     <  2:170432/150‑1 (MQ=255)

                                                                                                                     aGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGt                                 >  2:467865/1‑150 (MQ=255)

                                                                                                                     aGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGt                                 <  1:880197/150‑1 (MQ=255)

                                                                                                                      gCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCCCGCGtt                                >  1:3632326/1‑150 (MQ=255)

                                                                                                                      gCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGtt                                >  1:181167/1‑150 (MQ=255)

                                                                                                                      gCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGtt                                >  2:4981772/1‑150 (MQ=255)

                                                                                                                       cAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCTGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGttt                               <  1:4560737/150‑1 (MQ=255)

                                                                                                                            cATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTtctc                          <  1:3172997/150‑1 (MQ=255)

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

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

                                                                                                                                        tGTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCa              >  2:3175494/1‑150 (MQ=255)

                                                                                                                                        tGTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCa              >  2:3768854/1‑150 (MQ=255)

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

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

                                                                                                                                         gTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAg             >  2:5319455/1‑150 (MQ=255)

                                                                                                                                         gTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAg             <  1:1265375/150‑1 (MQ=255)

                                                                                                                                          ttGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgc            >  2:3813707/1‑150 (MQ=255)

                                                                                                                                            gAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgcgc          >  2:3859551/1‑150 (MQ=255)

                                                                                                                                             aaaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGATGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCa         >  2:3319275/1‑150 (MQ=255)

                                                                                                                                             aaaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCa         <  1:2162806/150‑1 (MQ=255)

                                                                                                                                             aaaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCa         <  1:1572836/150‑1 (MQ=255)

                                                                                                                                              aaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAg        >  2:4422138/1‑150 (MQ=255)

                                                                                                                                              aaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAg        >  1:3830306/1‑150 (MQ=255)

                                                                                                                                               accaccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGt       >  1:3238569/1‑150 (MQ=255)

                                                                                                                                               accaccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGt       >  1:1783708/1‑150 (MQ=255)

                                                                                                                                               accaccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGt       >  1:4795603/1‑150 (MQ=255)

                                                                                                                                                ccaccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGtt      <  2:3948567/150‑1 (MQ=255)

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

                                                                                                                                                  accGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGtttt    <  2:3403058/150‑1 (MQ=255)

                                                                                                                                                   ccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCTAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGttttt   <  1:1117131/150‑1 (MQ=255)

                                                                                                                                                    cGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTc  <  1:326941/150‑1 (MQ=255)

                                                                                                                                                    cGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTc  <  1:3505726/150‑1 (MQ=255)

                                                                                                                                                    cGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTc  <  2:3030664/150‑1 (MQ=255)

TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTC  >  NC_002947/1244228‑1244525

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