BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	723,954	G→A	G83G (GGC→GGT)	kdpC ←	KdpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	723,954	0	G	A	98.2%	378.8 / NA	113	G83G (GGC→GGT)	kdpC	KdpC
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (53/58); total (54/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

TGGCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCGCCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTAA  >  NC_000913/723866‑724043

gggCATCCGGGTAAGCGGCGCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAg                                                                                       <  2:3302710/92‑1 (MQ=255)

  gCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGaa                                                                                     <  1:2309355/93‑1 (MQ=255)

  gCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGaa                                                                                     <  1:966338/93‑1 (MQ=255)

  gCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGaa                                                                                     <  2:3957105/93‑1 (MQ=255)

    aTCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGc                                                                                   <  2:2811953/93‑1 (MQ=255)

    aTCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGc                                                                                   >  1:1832045/1‑93 (MQ=255)

      ccGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCt                                                                                 <  1:2421689/93‑1 (MQ=255)

      ccGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCt                                                                                 <  2:1092283/93‑1 (MQ=255)

       cGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtg                                                                                <  2:3358600/93‑1 (MQ=255)

        gggTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtgt                                                                               >  1:320700/1‑93 (MQ=255)

        gggTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtgt                                                                               <  2:4214843/93‑1 (MQ=255)

         ggTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtgtg                                                                              <  1:1803829/93‑1 (MQ=255)

         ggTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtgtg                                                                              <  2:2928451/93‑1 (MQ=255)

           ttAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGa                                                                            <  2:3698865/93‑1 (MQ=255)

            tAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGAt                                                                           <  1:924899/93‑1 (MQ=255)

              gCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATta                                                                         >  1:2549596/1‑93 (MQ=255)

                ggCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATtata                                                                       <  2:4449009/93‑1 (MQ=255)

                 gCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAg                                                                      <  1:3852592/93‑1 (MQ=255)

                 gCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAg                                                                      <  1:818837/93‑1 (MQ=255)

                  cccGTAATGCAGCAACGCGTGCAGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAgg                                                                     <  2:2660355/93‑1 (MQ=255)

                   ccGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAggg                                                                    >  2:397054/1‑93 (MQ=255)

                   ccGTAATGCAGCAACGCGTGCGGCTATTAGCTTATCCAGCTCAGGGGTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAgcg                                                                    >  2:3318584/1‑91 (MQ=255)

                      tAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAt                                                                 >  1:1010169/1‑93 (MQ=255)

                        aTGCAGCATGGCGTACGCCTATTAGTTTATCCAGCTCAGGATTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAttt                                                               <  2:3116498/93‑1 (MQ=255)

                        aTGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAttt                                                               >  2:615343/1‑93 (MQ=255)

                        aTGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAttt                                                               >  2:572090/1‑93 (MQ=255)

                         tGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTc                                                              <  1:2877765/93‑1 (MQ=255)

                          gcagcaACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCt                                                             <  2:258781/93‑1 (MQ=255)

                            agcaACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCGCCAGAAGCCTGTGGATTATAGGGCATTTCTGc                                                           >  1:2290399/1‑93 (MQ=255)

                            agcaACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGc                                                           <  2:2451319/93‑1 (MQ=255)

                             acaACGCATGCGGATATTAGTTTATCCAGCGCAGGGTTACTGACCGCCAGATTGCTCCAACCAGAAGCCTGTGGATTATAGGGCATTTCTGcc                                                          <  2:1626685/92‑1 (MQ=255)

                              caACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCg                                                         <  2:3913818/93‑1 (MQ=255)

                               aaCGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGt                                                        >  1:2751580/1‑93 (MQ=255)

                                aCGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGtt                                                       <  2:1063569/93‑1 (MQ=255)

                                aCGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGtt                                                       >  2:3063597/1‑93 (MQ=255)

                                 cgcgTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTg                                                      <  1:22904/93‑1 (MQ=255)

                                  gcgTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGc                                                     <  2:1662446/93‑1 (MQ=255)

                                   cgTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGcc                                                    <  2:255283/93‑1 (MQ=255)

                                    gTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCg                                                   >  1:1491509/1‑93 (MQ=255)

                                     tGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGa                                                  <  1:1389280/93‑1 (MQ=255)

                                        ggCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACgg                                               <  1:1602256/93‑1 (MQ=255)

                                         gCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGGTTATAGGGCATTTCTGCCGTTGCCGACggg                                              >  1:3280308/1‑93 (MQ=255)

                                         gCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACggg                                              >  2:1749993/1‑93 (MQ=255)

                                            tttAGTTTATCCAGCGCAGGGTTACTCACCGGCAGATTGCTCCCACCAGAAGCCTGTGGATTATTGGGCATTTCTGCCGTTGCCGACGGGCGa                                           <  2:464804/92‑1 (MQ=255)

                                            aTTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGa                                           >  2:2429943/1‑93 (MQ=255)

                                            aTTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGa                                           >  1:2994157/1‑93 (MQ=255)

                                             ttAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGAc                                          <  2:1601479/93‑1 (MQ=255)

                                                gTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCAt                                       <  2:1571877/93‑1 (MQ=255)

                                                gTTGATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCAt                                       <  1:2223618/93‑1 (MQ=255)

                                                 tttATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATg                                      >  1:3342389/1‑93 (MQ=255)

                                                  ttATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGa                                     <  1:931783/93‑1 (MQ=255)

                                                  ttATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGa                                     >  1:1304640/1‑93 (MQ=255)

                                                      ccAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAAt                                 <  2:719069/93‑1 (MQ=255)

                                                       cAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATa                                <  2:3640575/93‑1 (MQ=255)

                                                       cAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATa                                >  1:1819465/1‑93 (MQ=255)

                                                       cAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGAGAACATGAAAATa                                >  1:4438267/1‑93 (MQ=255)

                                                       cAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTAGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATa                                <  1:2173994/93‑1 (MQ=255)

                                                         gCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGc                              <  1:1270713/93‑1 (MQ=255)

                                                         gCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGc                              >  1:3441510/1‑93 (MQ=255)

                                                         gCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGc                              <  2:3079553/93‑1 (MQ=255)

                                                         gCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGc                              <  1:1785342/93‑1 (MQ=255)

                                                          cTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGcc                             >  1:400881/1‑93 (MQ=255)

                                                          cTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGcc                             >  2:3066848/1‑93 (MQ=255)

                                                           tCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCg                            >  2:106709/1‑93 (MQ=255)

                                                           tCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCg                            >  2:2503845/1‑93 (MQ=255)

                                                            cANGGTTACTGACCGCCAGATTNCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGt                           <  1:3202918/93‑1 (MQ=255)

                                                            cAGGGTTACTGACCGCCAGATTGCTCCCGCCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGt                           <  2:1460025/93‑1 (MQ=255)

                                                            cAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGt                           <  2:2707151/93‑1 (MQ=255)

                                                            cAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGt                           >  2:3449535/1‑93 (MQ=255)

                                                             aGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGtt                          <  2:3413597/93‑1 (MQ=255)

                                                             aGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGtt                          >  2:4485459/1‑93 (MQ=255)

                                                              gggTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTg                         >  1:3269551/1‑93 (MQ=255)

                                                               ggTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGACGAAGGGCGACCATGAAAATAGCCGTTGc                        >  2:3320057/1‑93 (MQ=255)

                                                                gTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGcc                       <  1:4293907/93‑1 (MQ=255)

                                                                gTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGcc                       >  1:303415/1‑93 (MQ=255)

                                                                      gACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTaaa                 >  2:3455814/1‑93 (MQ=255)

                                                                       gccGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTaaaa                <  2:2137924/92‑1 (MQ=255)

                                                                       aCCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTaaaa                <  1:4056041/93‑1 (MQ=255)

                                                                       aCCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTaaaa                <  1:4042165/93‑1 (MQ=255)

                                                                        ccGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCAGGTAAAAt               <  2:54745/93‑1 (MQ=255)

                                                                         cGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAAtt              >  1:4126373/1‑93 (MQ=255)

                                                                          gCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTc             >  1:3394805/1‑93 (MQ=255)

                                                                          gCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTc             >  2:4144085/1‑93 (MQ=255)

                                                                          gCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTc             <  1:3139677/93‑1 (MQ=255)

                                                                          gCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTc             >  1:1923372/1‑93 (MQ=255)

                                                                           ccAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCCGTAAAATTCt            <  2:2867535/93‑1 (MQ=255)

                                                                              gATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGcc         >  2:566487/1‑93 (MQ=255)

                                                                              gATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGcc         >  1:4218290/1‑93 (MQ=255)

                                                                              gATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGcc         >  1:4158693/1‑93 (MQ=255)

                                                                              gATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGcc         >  2:4539024/1‑93 (MQ=255)

                                                                                 tGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGa      <  2:815666/93‑1 (MQ=255)

                                                                                   cTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGAtt    >  2:4239002/1‑93 (MQ=255)

                                                                                   cTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGAtt    >  1:1344997/1‑93 (MQ=255)

                                                                                   cTCCCACCAGAAGCCTGTGGAGTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGAtt    <  1:734842/93‑1 (MQ=255)

                                                                                    tCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTa   >  2:1258700/1‑93 (MQ=255)

                                                                                    tCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTa   >  2:2089224/1‑93 (MQ=255)

                                                                                    tCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTa   >  1:1141458/1‑93 (MQ=255)

                                                                                     cccACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTaa  >  2:2670035/1‑93 (MQ=255)

                                                                                     cccACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTaa  <  2:2214529/93‑1 (MQ=255)

                                                                                     cccACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTaa  <  1:3946175/93‑1 (MQ=255)

TGGCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCGCCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTAA  >  NC_000913/723866‑724043

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