BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,369,637	T→C	R345R (CGT→CGC)	mopA →	GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,369,637	0	T	C	98.9%	578.3 / ‑4.0	183	R345R (CGT→CGC)	mopA	GroEL
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (91/90); total (93/90)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 193 total aligned reads displayed.

GAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGG  >  NC_000913/4369545‑4369726

  aGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCgcgt                                                                                         >  1:2279376/1‑93 (MQ=255)

  aGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCgcgt                                                                                         <  2:698496/93‑1 (MQ=255)

  aGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCgcgt                                                                                         <  2:3558617/93‑1 (MQ=255)

   gACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGtt                                                                                        >  1:1059706/1‑93 (MQ=255)

    aCCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTg                                                                                       >  2:278863/1‑93 (MQ=255)

     ccTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGc                                                                                      >  1:4409780/1‑93 (MQ=255)

       tGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTc                                                                                    >  2:1952834/1‑93 (MQ=255)

        gggACAGGCTAAACGTGGTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCa                                                                                   >  2:4224760/1‑93 (MQ=255)

            cAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATc                                                                               <  1:4290503/93‑1 (MQ=255)

            cAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATc                                                                               <  1:1114012/93‑1 (MQ=255)

             aGGCTAAACGTGTTGTGATCAACACAGACACCACCACTATCATCTATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGCTGCTCAGATcc                                                                              >  1:550631/1‑93 (MQ=255)

             aGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATcc                                                                              >  1:1197504/1‑93 (MQ=255)

              ggCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAGGCTGCAATCCAGGGCCGCGTTGCTCAGATGCg                                                                             >  2:1328498/1‑93 (MQ=255)

                 tAAACGTGTTGTTATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTca                                                                          <  2:3147264/93‑1 (MQ=255)

                   aaCGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTcagc                                                                        >  1:3140891/1‑93 (MQ=255)

                    aCGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTcagca                                                                       <  1:2545921/93‑1 (MQ=255)

                    aCGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTcagca                                                                       <  1:3226730/93‑1 (MQ=255)

                      gtgtTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGa                                                                     <  2:1633277/93‑1 (MQ=255)

                        gttgtGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGAtt                                                                   >  1:414366/1‑93 (MQ=255)

                        gttgtGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGAtt                                                                   <  1:1453060/93‑1 (MQ=255)

                        gttatGATCAACAAAGACACCACCACTATCATCGATGCCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGAtt                                                                   <  1:3489650/93‑1 (MQ=255)

                          tgtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTga                                                                 >  1:872304/1‑93 (MQ=255)

                           gtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaa                                                                >  1:4459538/1‑93 (MQ=255)

                             gATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGTTCAGATCCGTCAGCAGATTgaaga                                                              >  1:3593795/1‑93 (MQ=255)

                             gATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaaga                                                              <  1:448207/93‑1 (MQ=255)

                              aTCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaagaa                                                             >  2:4437989/1‑93 (MQ=255)

                              aTCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaagaa                                                             <  1:3890501/93‑1 (MQ=255)

                                caacaaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGc                                                           <  1:4162512/93‑1 (MQ=255)

                                 aacaaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCa                                                          <  2:3076319/93‑1 (MQ=255)

                                 aacaaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCa                                                          <  2:616031/93‑1 (MQ=255)

                                   caaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAAc                                                        >  1:2741601/1‑93 (MQ=255)

                                    aaaGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACt                                                       <  2:2197642/93‑1 (MQ=255)

                                     aaGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACtt                                                      >  1:3569042/1‑93 (MQ=255)

                                     aaGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACtt                                                      <  1:4098126/93‑1 (MQ=255)

                                      aGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                                                     <  2:43103/93‑1 (MQ=255)

                                      aGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                                                     <  2:3970076/93‑1 (MQ=255)

                                      aGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                                                     <  2:2355583/93‑1 (MQ=255)

                                       gACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCt                                                    >  2:1374974/1‑93 (MQ=255)

                                       gACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCt                                                    >  1:4683989/1‑93 (MQ=255)

                                       gACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCt                                                    >  1:3460595/1‑93 (MQ=255)

                                        acacCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTg                                                   <  2:4749381/93‑1 (MQ=255)

                                          accaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                                                 <  2:2026732/93‑1 (MQ=255)

                                          accaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                                                 >  2:3565734/1‑93 (MQ=255)

                                           ccaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACt                                                <  1:1013332/93‑1 (MQ=255)

                                           ccaccacTATCATCGATGGCGTGGGAGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACt                                                >  2:1806945/1‑93 (MQ=255)

                                            caccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTa                                               <  2:400110/93‑1 (MQ=255)

                                            caccacTATCATCGATGGCGTGGGATAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCCGATTCCAGAAGAAACTTCTGACTa                                               >  2:3336901/1‑93 (MQ=255)

                                             accacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTAc                                              >  2:4045850/1‑93 (MQ=255)

                                               cacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGa                                            >  2:4109030/1‑93 (MQ=255)

                                                attATCATCGATGGCGGGGGCGAAGAAGCTGCAAGCCAGGGCCGCGTTGCTCAGAGCCTTCAGCAGATTGAAGAAGCAACTTCTGACTACGAc                                           >  2:3277395/3‑93 (MQ=255)

                                                 cTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTTAGCAGATTGAAGAAGCAACTTCTGACTACGAcc                                          >  1:4576601/1‑93 (MQ=255)

                                                 cTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGAcc                                          >  1:4278637/1‑93 (MQ=255)

                                                  tATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCg                                         >  1:1504578/1‑93 (MQ=255)

                                                   atgatcGATGCAGTGGGTGAAGAAGCTGCAATCCAGGGCGGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGt                                        <  1:1728773/90‑1 (MQ=255)

                                                   atcatcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGt                                        <  1:2421437/93‑1 (MQ=255)

                                                    tcatgGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTg                                       <  2:664960/93‑1 (MQ=255)

                                                    tcatcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTg                                       <  2:3034418/93‑1 (MQ=255)

                                                       tcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaa                                    >  1:1154965/1‑93 (MQ=255)

                                                         gATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaaaa                                  <  2:3442129/93‑1 (MQ=255)

                                                         gATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaaaa                                  >  2:3847944/1‑93 (MQ=255)

                                                         gATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaaaa                                  >  2:52498/1‑93 (MQ=255)

                                                          aTGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAAc                                 >  1:1020687/1‑93 (MQ=255)

                                                          aTGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAAc                                 >  1:1082134/1‑93 (MQ=255)

                                                          aTGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAAc                                 >  1:988108/1‑93 (MQ=255)

                                                           tGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACt                                <  2:216734/93‑1 (MQ=255)

                                                             gCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGc                              <  1:1770375/93‑1 (MQ=255)

                                                               gTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAg                            <  2:3941361/93‑1 (MQ=255)

                                                               gTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAg                            <  2:46888/93‑1 (MQ=255)

                                                                 gggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGa                          <  2:3161867/93‑1 (MQ=255)

                                                                  ggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCTGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGaa                         <  2:2938425/93‑1 (MQ=255)

                                                                  ggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGaa                         <  1:686583/93‑1 (MQ=255)

                                                                   gTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAc                        <  1:2067188/93‑1 (MQ=255)

                                                                    tGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcg                       >  1:3714855/1‑93 (MQ=255)

                                                                    tGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcg                       <  2:1104294/93‑1 (MQ=255)

                                                                       agaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGt                    >  1:586585/1‑93 (MQ=255)

                                                                       agaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGt                    >  1:334991/1‑93 (MQ=255)

                                                                        gaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa                   <  1:3924366/93‑1 (MQ=255)

                                                                        gaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa                   >  1:51860/1‑93 (MQ=255)

                                                                        gaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa                   >  2:2899541/1‑93 (MQ=255)

                                                                          agCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGc                 <  1:3004731/93‑1 (MQ=255)

                                                                          agCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGc                 <  1:3827113/93‑1 (MQ=255)

                                                                          agCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGc                 >  2:4340900/1‑93 (MQ=255)

                                                                           gCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTATGACCGTGAAAAACTGCAGGAACGCGTAGCg                >  2:4636995/1‑93 (MQ=255)

                                                                              gCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGaaa             >  1:1706014/1‑93 (MQ=255)

                                                                              gCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGaaa             >  2:1751948/1‑93 (MQ=255)

                                                                              gCAATCCAGGGCCGCGCTGCTCAGATCCGTCAGCAGATTGAAGAAACAACTTCGGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCCaaa             >  2:1141611/1‑93 (MQ=255)

                                                                               cAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAAc            <  1:3823660/93‑1 (MQ=255)

                                                                                aaTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACt           >  1:93677/1‑93 (MQ=255)

                                                                                aaTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACt           <  2:1943579/93‑1 (MQ=255)

                                                                                aaTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACt           <  1:1918228/93‑1 (MQ=255)

                                                                                 aTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg          <  2:3753929/93‑1 (MQ=255)

                                                                                 aTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg          <  1:2024251/93‑1 (MQ=255)

                                                                                  tCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTgg         >  2:4237823/1‑93 (MQ=255)

                                                                                   ccAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGc        >  2:922961/1‑93 (MQ=255)

                                                                                    cAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGcg       >  2:1532172/1‑92 (MQ=255)

                                                                                    cAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCa       >  1:2861873/1‑93 (MQ=255)

                                                                                    cAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCa       <  1:2224805/93‑1 (MQ=255)

                                                                                      gggCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAgg     <  2:1543904/93‑1 (MQ=255)

                                                                                       ggCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAggc    <  2:1366233/93‑1 (MQ=255)

                                                                                        gCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAggcg   <  1:2139917/93‑1 (MQ=255)

GAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGG  >  NC_000913/4369545‑4369726

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