BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,113,475	G→A	16.2%	L61L (TTG→TTA) ‡	uspD →	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,475	0	G	A	16.2%	516.0 / 76.9	218	L61L (TTG→TTA)	uspD	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base G (93/88); new base A (20/15); total (115/103)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

Only 100 of 232 total aligned reads displayed.

TAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACA  >  NC_000913/4113379‑4113574

tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                  >  1:434476/1‑100 (MQ=255)

tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                  >  1:1125545/1‑100 (MQ=255)

tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                  >  1:3877462/1‑100 (MQ=255)

tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                  <  2:3208725/100‑1 (MQ=255)

 aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                  <  2:1351995/99‑1 (MQ=255)

   tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTgaag                                                                                               <  2:168610/100‑1 (MQ=255)

   tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTgaag                                                                                               <  2:4032792/100‑1 (MQ=255)

   tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTgaag                                                                                               <  1:733612/100‑1 (MQ=255)

   tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTgaag                                                                                               <  2:576970/100‑1 (MQ=255)

     aCGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTgaagaa                                                                                             <  1:863671/100‑1 (MQ=255)

      cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAAt                                                                                            <  2:3759529/100‑1 (MQ=255)

      cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAAt                                                                                            >  2:3015927/1‑100 (MQ=255)

          caacTGACGTTAATTCATAGTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACGTCCTTTAAACAGAAGATATGCTTCAAC‑AATTCAAGAATa                                                                                           <  1:2495134/94‑1 (MQ=255)

          cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGt                                                                                        <  2:4051715/100‑1 (MQ=255)

          cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAC‑AATTGAAGAATatct                                                                                        <  2:1415794/100‑4 (MQ=38)

           aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTc                                                                                       >  2:3589699/1‑100 (MQ=255)

            ccTGACTTTAATTCATATTGATGATGCCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAC‑AATTGAAGAATatcttc                                                                                      <  1:1884523/100‑6 (MQ=37)

            ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCg                                                                                      >  2:1918410/1‑100 (MQ=255)

            ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCg                                                                                      <  1:611450/100‑1 (MQ=255)

             cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCgg                                                                                     <  2:1017226/100‑1 (MQ=255)

              tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCgg                                                                                     <  1:147037/99‑1 (MQ=255)

              tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGa                                                                                    >  1:816380/1‑100 (MQ=255)

                 cGTTAATTCATAATGAAGATGCCTTAAGCGAGTTGTATCCGGAAATCTACTTCCCTGCAACAGAAGATATACTTCAAT‑TGTTGAAGAATAAGTCgg                                                                                     <  1:1520341/96‑1 (MQ=255)

                   ttAATTCATATTGATGATGGTTTAAGCTAGTTGTATTCTTGTATTTACTTTTCTGCAACAGAAGATATTCTTTAAC‑AATTGAAGAATatct                                                                                        <  1:3841335/91‑4 (MQ=255)

                           tatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGtat                                                                       >  1:4113267/1‑100 (MQ=255)

                            atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                      <  1:1942226/100‑1 (MQ=255)

                             ttGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGtatag                                                                     >  2:2378603/1‑99 (MQ=255)

                             ttGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATaa                                                                     >  1:1543853/1‑100 (MQ=255)

                              tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGGTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                      <  1:192963/98‑1 (MQ=255)

                              tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                    <  1:882077/100‑1 (MQ=255)

                               gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCGGTATAAAc                                                                   >  2:1869366/1‑100 (MQ=255)

                                atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                  <  2:1988761/100‑1 (MQ=255)

                                 tgatgGCTTAAGCGAGTTGTATCCGGGGATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                                 <  1:306043/100‑1 (MQ=255)

                                 tgatgGCTTAAGCAAGTTGTATCCGGGTATCGACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                                 >  2:2973602/1‑100 (MQ=255)

                                  gatgGCTTAAGCGAGTTGTACCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGa                                                                >  2:2649908/1‑100 (MQ=255)

                                      gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                            >  1:2778934/1‑100 (MQ=255)

                                       cTTGAGCGAGTTGTGTTCGAGTATCCAGGTCCCGGCTACAGAAGATATTCTGCAAT‑TGTTGATGAGGAAGTCGGATAACAAGCGGTATAAAATGACGaa                                                            >  1:2564403/1‑99 (MQ=255)

                                       cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaa                                                           <  1:474638/100‑1 (MQ=255)

                                            gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                      >  2:491208/1‑100 (MQ=255)

                                            gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                      >  2:3122749/1‑100 (MQ=255)

                                              gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGGATAATATTc                                                    <  1:3443081/100‑1 (MQ=255)

                                              gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                    <  2:2427826/100‑1 (MQ=255)

                                               aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                   <  1:2701704/100‑1 (MQ=255)

                                                gTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                                  >  2:3504410/1‑100 (MQ=255)

                                                gTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                                  >  2:1759192/1‑100 (MQ=255)

                                                  tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                                >  1:993171/1‑100 (MQ=255)

                                                    tATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGc                                              <  2:1217185/100‑1 (MQ=255)

                                                    tATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGc                                              <  1:2485373/100‑1 (MQ=255)

                                                     gtCCGGATATCTACTTCGCTGCTGCAGATCATATTTTTCGATGGGCTGAA‑CATCAGCGGGGTAGCTAGCTGCATTAACTGACGAAAAATATTCAATAGcc                                             >  1:1676334/2‑100 (MQ=255)

                                                     aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                             >  2:2703348/1‑100 (MQ=255)

                                                       gcGGGTATCTACTTCCCTGCACCCGAAGAGATTGTCCAAT‑ATTTTTTCAATCATTCGGATAACAAGCTGTAACAACTGACGAGAAGTATTCAAt                                                 >  1:3694856/2‑94 (MQ=255)

                                                        cGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGaa                                          >  2:3821581/1‑100 (MQ=255)

                                                         gggTATCTAGTTCCATGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                         <  1:353229/100‑1 (MQ=255)

                                                          ggTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                        >  2:4252555/1‑100 (MQ=255)

                                                            tATCTACTTCCCTGCAACAGAAGATATTCTGCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATAGCCGAAGACa                                      >  1:2896981/1‑100 (MQ=255)

                                                              tCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGACTAAGTCGGATAACAAGCTGTATAAACCGACGAAAAATATTCAATGGCCGAAGACaaa                                    <  2:4188116/100‑1 (MQ=255)

                                                               cTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaaa                                   <  2:707884/100‑1 (MQ=255)

                                                                tACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                  <  1:1100687/100‑1 (MQ=255)

                                                                 aCTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGACAGAAAACAAAACt                                 >  1:3814695/1‑100 (MQ=255)

                                                                   ttCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                               <  1:2122986/100‑1 (MQ=255)

                                                                    tCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                              >  1:3134699/1‑100 (MQ=255)

                                                                    tCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                              >  2:2945751/1‑100 (MQ=255)

                                                                    tCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGACGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaaaat                                 >  1:4274824/1‑95 (MQ=255)

                                                                     cccTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                             <  1:2175993/100‑1 (MQ=255)

                                                                      ccTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                            >  2:3949408/1‑100 (MQ=255)

                                                                       cTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                           >  2:476467/1‑100 (MQ=255)

                                                                       cTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                           <  2:3877664/100‑1 (MQ=255)

                                                                       cTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                           >  1:2327610/1‑100 (MQ=255)

                                                                       cTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                           >  2:2729010/1‑100 (MQ=255)

                                                                        tGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                          <  2:1639914/100‑1 (MQ=255)

                                                                         gCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTg                         >  2:2820563/1‑100 (MQ=255)

                                                                            atccgacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                      >  2:639871/9‑100 (MQ=37)

                                                                            aCAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                      <  1:168804/100‑1 (MQ=255)

                                                                              ccgacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                    >  1:2029479/7‑100 (MQ=38)

                                                                              agaagaTATTCTTGAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                    <  1:1069173/100‑1 (MQ=255)

                                                                              agaagaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                    <  1:2229394/100‑1 (MQ=255)

                                                                              agaagaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                    >  1:694021/1‑100 (MQ=255)

                                                                               gaagaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                   <  1:2584328/100‑1 (MQ=255)

                                                                                gacttatTCTTCAAC‑AATTGAAGAATCAGTCCGATACCAACCTGTATAAACTGACGAAAAATCTCCAACGCCCGAAGCCAACCCTCCGTATTGAAcgc                    >  1:2408224/5‑98 (MQ=16)

                                                                                gacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                   >  2:1095057/5‑99 (MQ=38)

                                                                                gacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                   >  1:2043207/5‑99 (MQ=38)

                                                                                gacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  >  1:1217439/5‑100 (MQ=38)

                                                                                gacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  >  1:3234989/5‑100 (MQ=38)

                                                                                 agaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 >  2:851494/1‑100 (MQ=255)

                                                                                 acttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 >  2:3311468/4‑100 (MQ=38)

                                                                                 acttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 >  1:1487762/4‑100 (MQ=38)

                                                                                  cttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag                >  1:3664195/3‑100 (MQ=255)

                                                                                    tatTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAAGCGGAGaa              >  1:4070321/1‑100 (MQ=255)

                                                                                      ttcttcAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCAGAGAAAt            <  2:2313334/100‑1 (MQ=255)

                                                                                          tcAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg        <  1:2820547/100‑1 (MQ=255)

                                                                                           cAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg       >  2:3734185/1‑100 (MQ=255)

                                                                                           cAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg       >  2:2096871/1‑100 (MQ=255)

                                                                                           cAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg       <  1:4194391/100‑1 (MQ=255)

                                                                                            aaT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa      <  1:3410365/100‑1 (MQ=255)

                                                                                            aaT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa      <  1:3754047/100‑1 (MQ=255)

                                                                                             aT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa     >  2:3418725/1‑100 (MQ=255)

                                                                                             aT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa     >  1:4244273/1‑100 (MQ=255)

                                                                                                tgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAac   <  1:2458642/100‑1 (MQ=255)

                                                                                                tgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAac   <  2:1543858/100‑1 (MQ=255)

                                                                                                 gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaaaga  >  2:1454609/1‑98 (MQ=255)

TAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACA  >  NC_000913/4113379‑4113574

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