BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	144	A→G	intergenic (–/‑46)	– / → thrL	–/thr operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	144	0	A	G	98.6%	1188.2 / 2.6	349	intergenic (–/‑46)	–/thrL	–/thr operon leader peptide
Reads supporting (aligned to +/- strand): ref base A (4/1); new base G (201/142); total (205/144)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 363 total aligned reads displayed.

AGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTAC  >  NC_000913/53‑233

 gagTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATggg                                                                                         >  2:8018305/1‑93 (MQ=255)

   gtgtCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAGATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCa                                                                                       <  2:133364/93‑1 (MQ=255)

   gtgtCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCa                                                                                       >  2:5726302/1‑93 (MQ=255)

     gtCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATa                                                                                     >  2:4964574/1‑93 (MQ=255)

      tCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAg                                                                                    <  1:3286689/93‑1 (MQ=255)

        tGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAgcg                                                                                  <  2:349476/93‑1 (MQ=255)

        tGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAgcg                                                                                  <  2:375089/93‑1 (MQ=255)

         gATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAgcgc                                                                                 >  1:1059176/1‑93 (MQ=255)

         gATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAgcgc                                                                                 <  1:1318148/93‑1 (MQ=255)

          aTAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGca                                                                                >  1:2917312/1‑93 (MQ=255)

          aTAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGca                                                                                >  1:2451232/1‑93 (MQ=255)

           tAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGcac                                                                               >  2:7205690/1‑93 (MQ=255)

            agcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGcaca                                                                              <  2:7739997/93‑1 (MQ=255)

            agcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGcaca                                                                              >  1:7709453/1‑93 (MQ=255)

            agcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGcaca                                                                              <  1:6667324/93‑1 (MQ=255)

             gcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacag                                                                             >  1:3171223/1‑93 (MQ=255)

             gcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacag                                                                             >  1:5317803/1‑93 (MQ=255)

              cagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacaga                                                                            <  2:8471526/93‑1 (MQ=255)

               agcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagac                                                                           >  1:1713167/1‑93 (MQ=255)

               agcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagac                                                                           >  1:7425986/1‑93 (MQ=255)

               agcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagac                                                                           >  1:5431270/1‑93 (MQ=255)

                gcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagaca                                                                          >  1:4338591/1‑93 (MQ=255)

                gcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagaca                                                                          >  2:3009512/1‑93 (MQ=255)

                  ttCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCAcagacaga                                                                        <  2:2693627/93‑1 (MQ=255)

                  ttCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCAcagacaga                                                                        >  2:6733200/1‑93 (MQ=255)

                  ttCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCAcagacaga                                                                        <  1:5734538/93‑1 (MQ=255)

                   tCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGAt                                                                       >  1:5891849/1‑93 (MQ=255)

                   tCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGAt                                                                       >  1:2744728/1‑93 (MQ=255)

                     tGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATaa                                                                     <  1:8697374/93‑1 (MQ=255)

                     tGAACTGGTTACCTGCCGCGGGTAAAATAAAATTTTATTGACTTAGGGCACTTAATACTTTAACCAATATAGGCATACCGCAATGACAGATaa                                                                     >  1:6746447/1‑93 (MQ=255)

                      gAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATaaa                                                                    <  1:4925032/93‑1 (MQ=255)

                       aaCTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATaaaa                                                                   <  1:899372/93‑1 (MQ=255)

                            gTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTAc                                                              <  1:7157829/93‑1 (MQ=255)

                            gTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTAc                                                              >  2:3943796/1‑93 (MQ=255)

                             ttACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACa                                                             <  2:3085279/93‑1 (MQ=255)

                               aCCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACaga                                                           <  1:6023098/93‑1 (MQ=255)

                                ccTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACagag                                                          <  1:6255893/93‑1 (MQ=255)

                                ccTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACagag                                                          <  2:5413477/93‑1 (MQ=255)

                                ccTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACagag                                                          >  2:7550088/1‑93 (MQ=255)

                                 cTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGt                                                         <  1:1889441/93‑1 (MQ=255)

                                  tGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGAGCACAGACAGATAAAAATTACAGAGTa                                                        >  2:6047334/1‑93 (MQ=255)

                                    ccGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTaca                                                      <  1:8753089/93‑1 (MQ=255)

                                     cGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTacac                                                     <  1:703397/93‑1 (MQ=255)

                                       tgggTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaa                                                   >  2:6004303/4‑93 (MQ=255)

                                        gAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaac                                                  >  1:4680749/1‑93 (MQ=255)

                                        gAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaac                                                  >  2:8735852/1‑93 (MQ=255)

                                         aGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaaca                                                 >  1:3616933/1‑93 (MQ=255)

                                         aGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaaca                                                 <  1:3082713/93‑1 (MQ=255)

                                         aGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaaca                                                 <  2:1408024/93‑1 (MQ=255)

                                          gTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACAt                                                <  1:3846910/93‑1 (MQ=255)

                                          gTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACAt                                                >  2:1142181/1‑93 (MQ=255)

                                          gTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACAt                                                >  2:3151106/1‑93 (MQ=255)

                                           tAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATc                                               <  1:6611842/93‑1 (MQ=255)

                                           tAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATc                                               >  2:6200843/1‑93 (MQ=255)

                                           tAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATc                                               >  2:4177751/1‑93 (MQ=255)

                                            aaaTTAAAATTTTATTGACTTAGGTCACTAAATCCTTTAACCAATATGGGCATATCGCACAGACAGATAAAAATTACAGAGTACACAACATcc                                              <  1:5722781/93‑1 (MQ=255)

                                            aaaTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATcc                                              <  2:5272773/93‑1 (MQ=255)

                                            aaaTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATcc                                              >  2:5308078/1‑93 (MQ=255)

                                             aaTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCa                                             >  1:760690/1‑93 (MQ=255)

                                               ttAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATg                                           <  2:8315504/93‑1 (MQ=255)

                                                 aaaaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGaa                                         >  2:2424284/1‑93 (MQ=255)

                                                   aaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAAc                                       >  2:216984/1‑93 (MQ=255)

                                                   aaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAAc                                       >  2:7042692/1‑93 (MQ=255)

                                                     ttttATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGc                                     <  2:1994295/93‑1 (MQ=255)

                                                     ttttATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGc                                     >  2:4498791/1‑93 (MQ=255)

                                                     ttttATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGc                                     >  1:8510723/1‑93 (MQ=255)

                                                       ttATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCAt                                   <  1:5926635/93‑1 (MQ=255)

                                                       ttATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCAt                                   <  1:230116/93‑1 (MQ=255)

                                                        tATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCAtt                                  <  1:5107401/93‑1 (MQ=255)

                                                         aTTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTa                                 >  1:6476252/1‑93 (MQ=255)

                                                            gACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGca                              <  2:1450749/93‑1 (MQ=255)

                                                            gACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGca                              >  1:4205000/1‑93 (MQ=255)

                                                            gACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGca                              >  1:2647775/1‑93 (MQ=255)

                                                            gACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGca                              >  2:5162526/1‑93 (MQ=255)

                                                             aCTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTAAAGAGTACACAACATACATGAAACGCATTAGcac                             >  2:6561108/1‑93 (MQ=255)

                                                              cTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGcacc                            >  2:1319131/1‑93 (MQ=255)

                                                              cTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGaacc                            >  2:2915287/1‑93 (MQ=255)

                                                                 aGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCaccacc                         <  1:2797886/93‑1 (MQ=255)

                                                                    tCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCAtt                      <  2:4717247/93‑1 (MQ=255)

                                                                       cTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTacc                   <  1:1111054/93‑1 (MQ=255)

                                                                       cTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTacc                   >  2:1891105/1‑93 (MQ=255)

                                                                       cTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTacc                   >  1:826332/1‑93 (MQ=255)

                                                                        tAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTacca                  >  1:5588969/1‑93 (MQ=255)

                                                                         aaaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccac                 <  2:976886/93‑1 (MQ=255)

                                                                         aaaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccac                 >  2:4345481/1‑93 (MQ=255)

                                                                         aaaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccac                 >  2:6320380/1‑93 (MQ=255)

                                                                         aaaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccac                 >  2:1584435/1‑93 (MQ=255)

                                                                         aaaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccac                 >  2:6882090/1‑93 (MQ=255)

                                                                         aaaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccac                 <  1:8616144/93‑1 (MQ=255)

                                                                          aaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccacc                >  1:8437359/1‑93 (MQ=255)

                                                                          aaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccacc                >  2:876739/1‑93 (MQ=255)

                                                                             aCTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccaccacc             <  1:8048625/93‑1 (MQ=255)

                                                                               tttAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACcaccat           <  2:5075315/93‑1 (MQ=255)

                                                                                   aCCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACcaccatcacc       >  2:4028096/1‑93 (MQ=255)

                                                                                    ccAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACcaccatcacca      <  1:1388560/93‑1 (MQ=255)

                                                                                    ccAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACcaccatcacca      <  2:4550395/93‑1 (MQ=255)

                                                                                      aaTATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCAtt    <  2:804190/93‑1 (MQ=255)

                                                                                       atatGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTa   <  1:2047053/93‑1 (MQ=255)

                                                                                        tatgGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTAc  >  1:344875/1‑93 (MQ=255)

                                                                                        tatgGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTAc  >  2:4852916/1‑93 (MQ=255)

AGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTAC  >  NC_000913/53‑233

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