BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,642,279	A→G	V30A (GTC→GCC)	b1560 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,642,279	0	A	G	99.6%	719.2 / NA	224	V30A (GTC→GCC)	b1560	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (131/92); minor base T (0/1); total (131/93)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

GATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACG  >  NC_000913/1642189‑1642368

 aTGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGgcca                                                                                        <  2:3318153/93‑1 (MQ=255)

  tGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGgccac                                                                                       <  2:4087517/93‑1 (MQ=255)

  tGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGgccac                                                                                       <  1:6264945/93‑1 (MQ=255)

   ggCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGTATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGgccacc                                                                                      >  1:2268793/1‑93 (MQ=255)

    gCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGccacca                                                                                     >  1:1557331/1‑93 (MQ=255)

    gCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGccacca                                                                                     <  2:1859179/93‑1 (MQ=255)

    gCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGccacca                                                                                     >  1:3075581/1‑93 (MQ=255)

      aGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGc                                                                                   >  2:3688535/1‑93 (MQ=255)

       gATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGcc                                                                                  <  1:5178668/93‑1 (MQ=255)

       gATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGcc                                                                                  <  1:5515373/93‑1 (MQ=255)

        aTTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGccc                                                                                 <  1:4039600/93‑1 (MQ=255)

        aTTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGccc                                                                                 >  1:3753833/1‑93 (MQ=255)

         ttACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCa                                                                                >  1:3269695/1‑93 (MQ=255)

            cGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGc                                                                             <  1:6271628/93‑1 (MQ=255)

             gCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCt                                                                            >  1:4380797/1‑93 (MQ=255)

              cATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTc                                                                           >  1:5435763/1‑93 (MQ=255)

                     tttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCAcc                                                                    <  2:2951947/93‑1 (MQ=255)

                      ttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCa                                                                   <  1:264820/93‑1 (MQ=255)

                      ttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCa                                                                   >  1:539612/1‑93 (MQ=255)

                      ttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCa                                                                   <  2:5703893/93‑1 (MQ=255)

                         gttccggttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGaa                                                                >  2:3627364/1‑93 (MQ=255)

                            ccggttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAg                                                             <  2:5528211/93‑1 (MQ=255)

                              ggttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGt                                                           >  2:6178592/1‑93 (MQ=255)

                               gttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTc                                                          <  1:646711/93‑1 (MQ=255)

                               gttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTc                                                          >  1:3826621/1‑93 (MQ=255)

                                ttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCt                                                         <  2:3062480/93‑1 (MQ=255)

                                ttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCt                                                         <  1:1419093/93‑1 (MQ=255)

                                  ccaccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTaa                                                       <  2:1440302/93‑1 (MQ=255)

                                  ccaccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTaa                                                       >  1:1238049/1‑93 (MQ=255)

                                   caccaGCATTCGAGGATTTTGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAg                                                      <  2:3299187/93‑1 (MQ=255)

                                   caccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAg                                                      >  2:2144117/1‑93 (MQ=255)

                                   caccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAg                                                      <  2:472947/93‑1 (MQ=255)

                                    accaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGt                                                     >  2:2360356/1‑93 (MQ=255)

                                    accaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGt                                                     >  2:2849196/1‑93 (MQ=255)

                                    accaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGt                                                     >  1:4920467/1‑93 (MQ=255)

                                     ccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTa                                                    >  2:6020890/1‑93 (MQ=255)

                                     ccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTa                                                    >  1:3644547/1‑93 (MQ=255)

                                      caGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTaa                                                   <  1:3105954/93‑1 (MQ=255)

                                      caGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTaa                                                   >  2:5608412/1‑93 (MQ=255)

                                      caGCACTCGAGGATTATGAAATATCTGTATGGATGGACGGCCCGGCTTAAGGTCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTaa                                                   <  1:1352494/93‑1 (MQ=255)

                                        gCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATa                                                 >  2:3140830/1‑93 (MQ=255)

                                        gCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATa                                                 >  1:4315326/1‑93 (MQ=255)

                                         cACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATAc                                                >  1:5656474/1‑93 (MQ=255)

                                          aCTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATAcc                                               <  1:5151303/93‑1 (MQ=255)

                                          aCTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATAcc                                               <  1:4722467/93‑1 (MQ=255)

                                           cTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATAccc                                              <  1:405913/93‑1 (MQ=255)

                                           cTCGAGGATGATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATAccc                                              <  1:3930439/93‑1 (MQ=255)

                                            tCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCg                                             >  2:5222577/1‑93 (MQ=255)

                                             cGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGc                                            <  1:5993350/93‑1 (MQ=255)

                                              gAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCa                                           >  2:444272/1‑93 (MQ=255)

                                              gAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCa                                           >  2:2992898/1‑93 (MQ=255)

                                                ggATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACg                                         >  1:4140224/1‑93 (MQ=255)

                                                ggATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACg                                         >  2:2493084/1‑93 (MQ=255)

                                                   ttATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTa                                      <  2:3871644/93‑1 (MQ=255)

                                                    tATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTAc                                     >  1:19283/1‑93 (MQ=255)

                                                    tATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTAc                                     >  2:5458425/1‑93 (MQ=255)

                                                    tATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTAc                                     >  2:4811836/1‑93 (MQ=255)

                                                      tGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCt                                   >  1:2961995/1‑93 (MQ=255)

                                                      tGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCt                                   >  1:28447/1‑93 (MQ=255)

                                                       gAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTc                                  <  1:1636447/93‑1 (MQ=255)

                                                         aaTATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCa                                >  2:778181/1‑93 (MQ=255)

                                                          atatCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAg                               >  1:982777/1‑93 (MQ=255)

                                                           tatCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGa                              <  2:1049541/93‑1 (MQ=255)

                                                            atCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGAt                             <  1:3130223/93‑1 (MQ=255)

                                                            atCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGAt                             >  1:5036856/1‑93 (MQ=255)

                                                             tCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCTGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATc                            <  1:4552762/93‑1 (MQ=255)

                                                               tGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCg                          >  1:2116856/1‑93 (MQ=255)

                                                                   tGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTCCCTCCAGATCCGTTGc                      <  2:2130890/93‑1 (MQ=255)

                                                                   tGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGc                      >  1:2701230/1‑93 (MQ=255)

                                                                    ggATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCt                     >  1:588032/1‑93 (MQ=255)

                                                                       ttCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGTCACCAGGACAGGTCTAAGCAATACCCGCTCGTTATCTCCAGATCCGTTGCTGGa                  >  2:2513143/1‑93 (MQ=255)

                                                                       ttCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTgag                  >  1:6051964/1‑91 (MQ=255)

                                                                       ttCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGa                  >  1:912340/1‑93 (MQ=255)

                                                                       ttCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGa                  >  2:3592305/1‑93 (MQ=255)

                                                                        tCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGaa                 >  1:4224749/1‑93 (MQ=255)

                                                                        tCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGaa                 >  1:353565/1‑93 (MQ=255)

                                                                         cACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGaaa                <  1:1829811/93‑1 (MQ=255)

                                                                         cACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGaaa                <  2:1310540/93‑1 (MQ=255)

                                                                          aCGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAg               >  2:627036/1‑93 (MQ=255)

                                                                          aCGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAg               >  2:5697566/1‑93 (MQ=255)

                                                                           cGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGt              <  2:2096396/93‑1 (MQ=255)

                                                                           cGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGt              <  1:464198/93‑1 (MQ=255)

                                                                            ggCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTg             >  2:4005289/1‑93 (MQ=255)

                                                                            ggCCCGGCTGAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTg             <  2:1754897/93‑1 (MQ=255)

                                                                             gCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGc            <  2:5072488/93‑1 (MQ=255)

                                                                             gCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGc            >  2:4396076/1‑93 (MQ=255)

                                                                              cccGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCg           <  1:5344982/93‑1 (MQ=255)

                                                                              cccGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCg           <  1:1079874/93‑1 (MQ=255)

                                                                               ccGGCTTAAGGGCCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCgg          <  2:4722945/93‑1 (MQ=255)

                                                                                cGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCggg         >  1:3729432/1‑93 (MQ=255)

                                                                                 ggCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGa        <  1:527981/93‑1 (MQ=255)

                                                                                 ggCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGa        <  1:50305/93‑1 (MQ=255)

                                                                                 ggCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGa        <  2:575304/93‑1 (MQ=255)

                                                                                  gCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGAc       >  1:3797493/1‑93 (MQ=255)

                                                                                   cTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACg      <  1:4237622/93‑1 (MQ=255)

                                                                                    ttAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGc     >  2:2660364/1‑93 (MQ=255)

                                                                                     tAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCa    >  2:5156771/1‑93 (MQ=255)

                                                                                     tAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCa    <  2:1180500/93‑1 (MQ=255)

                                                                                     tAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCa    >  1:3461042/1‑93 (MQ=255)

                                                                                     tAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAAGGCGGGACGCa    >  2:3980960/1‑93 (MQ=255)

GATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACG  >  NC_000913/1642189‑1642368

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