BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,030,523	A→G	R402R (CGA→CGG)	ygfU →	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,030,523	0	A	G	82.3%	475.9 / 99.9	203	R402R (CGA→CGG)	ygfU	putative permease
Reads supporting (aligned to +/- strand): ref base A (14/22); new base G (79/88); total (93/110)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

Only 100 of 209 total aligned reads displayed.

AAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGACGCTCTCTC  >  NC_000913/3030432‑3030614

aaaTGGCGGTGCTGGTAGGCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGaa                                                                                            >  1:6413424/1‑93 (MQ=255)

 aaTGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGAt                                                                                           <  2:3186271/93‑1 (MQ=255)

 aaTGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGAt                                                                                           >  2:4406721/1‑93 (MQ=255)

  aTGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTTGCGACAGGGATTCGGAtt                                                                                          >  2:3787739/1‑93 (MQ=255)

  aTGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGAtt                                                                                          <  2:3896465/93‑1 (MQ=255)

   tGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGggattcggattc                                                                                         <  2:2950170/93‑1 (MQ=255)

        gTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTcg                                                                                    >  1:936611/1‑93 (MQ=255)

        gTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTcg                                                                                    <  1:2936400/93‑1 (MQ=255)

             ggTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTg                                                                               <  1:5693656/93‑1 (MQ=255)

             ggTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTg                                                                               <  1:4950496/93‑1 (MQ=255)

             ggTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTg                                                                               <  1:5209243/93‑1 (MQ=255)

             ggTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGAGACAGGGATTCGGATTCTGTCGCGCTg                                                                               >  2:6635742/1‑93 (MQ=255)

               tAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTa                                                                             <  1:6249185/93‑1 (MQ=255)

               tAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTa                                                                             <  2:4533907/93‑1 (MQ=255)

                aGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGAGACAGGGATTCGGATTCTGTCGCGCTGTaa                                                                            >  2:5346263/1‑93 (MQ=255)

                 gCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAAc                                                                           <  1:12076/93‑1 (MQ=255)

                    tCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTac                                                                        >  1:2026928/1‑93 (MQ=255)

                    tCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTac                                                                        <  2:1625515/93‑1 (MQ=255)

                      cATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTacac                                                                      >  2:4773168/1‑93 (MQ=255)

                      cATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTacac                                                                      <  1:5007364/93‑1 (MQ=255)

                       aTTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTAcacc                                                                     <  2:2416528/93‑1 (MQ=255)

                       aTTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTAcacc                                                                     <  2:333299/93‑1 (MQ=255)

                        ttCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTAcacca                                                                    >  2:2358862/1‑93 (MQ=255)

                            gCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCaa                                                                >  1:1119088/1‑93 (MQ=255)

                              aaTTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAAcc                                                              <  2:5635312/93‑1 (MQ=255)

                              aaTTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAAcc                                                              >  1:1013308/1‑93 (MQ=255)

                               aTTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCg                                                             >  2:6605802/1‑93 (MQ=255)

                                  tgtgCTGGGCGGCGCTGGTCTAGTGGTGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGtca                                                          <  2:6421227/93‑3 (MQ=255)

                                  tgtgCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTa                                                          <  1:812789/93‑1 (MQ=255)

                                    tgCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACa                                                        <  2:2239935/93‑1 (MQ=255)

                                     gCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACaa                                                       <  2:6080301/93‑1 (MQ=255)

                                      cTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAAc                                                      <  2:3988075/93‑1 (MQ=255)

                                        gggCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACct                                                    >  2:1197928/1‑93 (MQ=255)

                                         ggcggcGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACctc                                                   <  2:5863827/93‑1 (MQ=255)

                                         ggcggcGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACctc                                                   >  2:3580219/1‑93 (MQ=255)

                                          gcggcgCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACctct                                                  <  2:2362627/93‑1 (MQ=255)

                                           cggcgCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCta                                                 >  2:6301525/1‑93 (MQ=255)

                                           cggcgCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCta                                                 <  1:377784/93‑1 (MQ=255)

                                             gcgcTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCtata                                               <  1:5105377/93‑1 (MQ=255)

                                             gcgcTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCtata                                               <  2:1776410/93‑1 (MQ=255)

                                             gcgcTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCtata                                               >  1:1303012/1‑93 (MQ=255)

                                               gcTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATAtt                                             >  2:5108172/1‑93 (MQ=255)

                                               gcTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATAtt                                             <  1:3123290/93‑1 (MQ=255)

                                                cTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATtg                                            <  2:5282888/93‑1 (MQ=255)

                                                  ggTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATtgtg                                          >  2:3693767/1‑93 (MQ=255)

                                                   gTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTgg                                         <  2:885632/93‑1 (MQ=255)

                                                    tCTAGTGATGTTCGGCATGGTACTGGAGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGc                                        >  1:3979163/1‑93 (MQ=255)

                                                     cTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCg                                       <  1:1461859/93‑1 (MQ=255)

                                                     cTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCg                                       <  1:2868032/93‑1 (MQ=255)

                                                      tAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGa                                      <  1:3922045/93‑1 (MQ=255)

                                                      tAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGa                                      >  2:4063736/1‑93 (MQ=255)

                                                       aGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGAt                                     >  1:6086856/1‑93 (MQ=255)

                                                       aGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGAt                                     <  2:5011510/93‑1 (MQ=255)

                                                        gTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATc                                    <  2:6203892/93‑1 (MQ=255)

                                                        gTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATc                                    <  2:5155374/93‑1 (MQ=255)

                                                        gTGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATc                                    <  2:2743407/93‑1 (MQ=255)

                                                         tGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCa                                   <  2:1446454/93‑1 (MQ=255)

                                                         tGATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCa                                   <  1:3449576/93‑1 (MQ=255)

                                                         tGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCa                                   <  2:3088308/93‑1 (MQ=255)

                                                          gATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAg                                  <  1:3231528/93‑1 (MQ=255)

                                                          gATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAg                                  <  1:5447090/93‑1 (MQ=255)

                                                          gATGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAg                                  >  1:5759183/1‑93 (MQ=255)

                                                          gATGTTCGGCATGGTACTGGAGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAg                                  >  1:1344049/1‑93 (MQ=255)

                                                            tggtCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGtc                                <  1:2786166/90‑1 (MQ=255)

                                                            tGTTCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGtc                                <  2:2338305/93‑1 (MQ=255)

                                                              ttCGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGtctc                              <  2:3803644/93‑1 (MQ=255)

                                                                cGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCgg                            <  1:1087639/93‑1 (MQ=255)

                                                                cGGCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCgg                            <  2:1829615/93‑1 (MQ=255)

                                                                 ggCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGc                           >  2:2902600/1‑93 (MQ=255)

                                                                 ggCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGc                           <  2:4499909/93‑1 (MQ=255)

                                                                  gCATGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCg                          <  1:4241337/93‑1 (MQ=255)

                                                                    aTGGTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGtt                        >  1:1654494/1‑93 (MQ=255)

                                                                     tGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTg                       >  2:1865557/1‑93 (MQ=255)

                                                                       gTACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGc                     <  2:784142/93‑1 (MQ=255)

                                                                        tACTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCa                    <  2:1949499/93‑1 (MQ=255)

                                                                          cTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATg                  >  1:5265861/1‑93 (MQ=255)

                                                                          cTGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATg                  <  2:2103452/93‑1 (MQ=255)

                                                                           tGGCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGa                 <  2:5773554/93‑1 (MQ=255)

                                                                            ggCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGAc                <  1:4546549/93‑1 (MQ=255)

                                                                            ggCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGAc                >  1:1998355/1‑93 (MQ=255)

                                                                             gCGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACt               >  2:5940577/1‑93 (MQ=255)

                                                                              cGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTc              <  2:1255620/93‑1 (MQ=255)

                                                                              cGACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTc              <  2:1847182/93‑1 (MQ=255)

                                                                              cGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTc              <  1:5601922/93‑1 (MQ=255)

                                                                              cGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTc              >  1:1718331/1‑93 (MQ=255)

                                                                               gACAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTcc             >  1:2449753/1‑93 (MQ=255)

                                                                               gACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTcc             <  2:4910008/93‑1 (MQ=255)

                                                                                aCAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCg            <  2:3394985/93‑1 (MQ=255)

                                                                                 cAGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGa           >  1:4758249/1‑93 (MQ=255)

                                                                                  aGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGAc          <  1:4673993/93‑1 (MQ=255)

                                                                                  aGGGATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGAc          <  2:1302211/93‑1 (MQ=255)

                                                                                  aGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGAc          >  1:2708884/1‑93 (MQ=255)

                                                                                  aGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGAc          <  1:6605802/93‑1 (MQ=255)

                                                                                     gATTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGACGct       <  1:3669636/93‑1 (MQ=255)

                                                                                      aTTCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGACGctc      >  2:4676976/1‑93 (MQ=255)

                                                                                      aTTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGACGctc      >  2:4251322/1‑93 (MQ=255)

                                                                                        tCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGACGctctc    <  2:1827443/93‑1 (MQ=255)

                                                                                        tCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGACGctctc    >  1:5494238/1‑93 (MQ=255)

                                                                                        tCGGATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGACGctctc    <  1:2165959/93‑1 (MQ=255)

                                                                                          gAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGACGCtctctc  >  1:1817994/1‑93 (MQ=255)

AAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATATTGTGGCGATCAGTCTCGGCGTTGGCATGACTCCGACGCTCTCTC  >  NC_000913/3030432‑3030614

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